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Abstract 

The features of electromagnetic processes are considered which connected with 
finite size of space region in which final particles (photon, electron-positron pair) 
are formed. The longitudinal dimension of the region is known as the formation 
length. If some external agent is acting on an electron while traveling this distance 
the emission process can be disrupted. There are different agents: multiple scatter- 
ing of projectile, polarization of a medium, action of external fields, etc. The theory 
of radiation under influence of the multiple scattering, the Landau-Pomeranchuk- 
Migdal (LPM) effect, is presented. The probability of radiation is calculated with 
an accuracy up to "next to leading logarithm" and with the Coulomb corrections 
taken into account. The integral characteristics of bremsstrahlung are given, it is 
shown that the effective radiation length increases due to the LPM effect at high 
energy. The LPM effect for pair creation is also presented. The multiple scattering 
influences also on radiative corrections in a medium (and an external field too) in- 
cluding the anomalous magnetic moment of an electron and the polarization tensor 
as well as coherent scattering of a photon in a Coulomb field. The polarization of 
a medium alters the radiation probability in soft part of spectrum. Specific fea- 
tures of radiation from a target of finite thickness include: the boundary photon 
emission, interference effects for thin target, multi-photon radiation. The theory 
predictions are compared with experimental data obtained at SLAC and CERN 
SPS. For electron-positron colliding beams following items are discussed: the sepa- 
ration of coherent and incoherent mechanisms of radiation, the beam-size effect in 
bremsstrahlung, coherent radiation and mechanisms of electron-positron creation. 
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1 Introduction 



1.1 General outlook 

The electromagnetic interaction is of the local nature, e.g. the interaction Hamiltonian 
in quantum electrodynamics is 

H int {x) = e^(x)^(x)A^x), (1.1) 

i.e. the current in a point x interacts with an electromagnetic field A^(x) in the same 
point. However the process of photon radiation takes place not in one point but in some 
domain of space-time, in which a photon (an emitted wave in the classical language) is 
originating. The longitudinal dimension of this domain is called the formation (coherence) 
length. It is evident that the minimal size of the formation length is the wave length of the 
emitted photon. One can estimate the formation length from the phase factor entering 
the expression for the probability of photon emission. In the quasiclassical operator 
method (see Eq. (jA.l|) . Appendix A) this factor is (in this paper we employ units such 
that h — c — 1) 



exp 



is 



1.2) 



--(ujt-kr(t)) , 

where s(s') is the energy of the initial (final) electron, uo is the photon energy, s' = s — u, 
k is the wave vector of the photon, r(t) is the particle location on a classical trajectory, 
t is the time. The main contribution into the integral over time in the general expression 
for the probability of radiation (|A.1J) gives the region 



~(a;*-kr(f))~l, (1.3) 



s 



Taking into account that 



2 2 4:7ian e 



£(u))um, r(t) = vt, S(u) = 1 — — |, u>'Z = '"^'" e ; nv = t>cosf9, (1.4) 

u 2 m 

where S(u) is the dielectric constant, uq is the plasma frequency, n e is the electron density, 
we have from Eq.JOJ) for •§ <C 1 

2ss' 

t~tf = h( u ) = / 7^r> u p = ^o7- (1-5) 

urn 2 l + 7 2 tf 2 + -| 

In this paper we consider the processes at high energies, i.e. the Lorentz factor 7 = 
s/m 3> 1. Then from Eq. (|1.5J) follows: 

1. ultrarelativistic particle radiates into the narrow cone with the vertex angle 1? < 1/7 
along the momentum of the initial particle, the contribution of larger angles is 
suppressed because of shortening of the formation length; 
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2. the effect of the polarization of a medium described by the dielectric constant 8(uo) 
manifest itself for soft photons only when 00 < uj b) since on the Earth lj < 100 eV 
we have u b < eu /m < 2 • 10~ 4 e; 

3. when d < I/7, uo^/uo <C 1 the formation length is of the order 

^W = ^- (1-6) 

urn 2 

In the classical case uo <C e one has lj = 2^ 2 /uo. So for ultrarelativistic particles the 
formation length extends substantially For example for a e = 25 GeV electron emitting 
aw = 100 MeV photon, lf = 10/im, i.e.~ 10 5 interatomic distances. 

The formation length is important in many electromagnetic processes including cre- 
ation of electron-positron pair, magnetic bremsstrahlung, transition radiation, Cerenkov 
radiation. There is also a set of applications involving other forces. 

There is a number of physical interpretation of the formation length. 

• The momentum transfer q to a nucleus in the process of the photon emission by the 
particle with the momentum p 

q = p-p'-k, (1.7) 

takes the minimal value when it is longitudinal q\\ = q m in = p — p' — k, and in this 
case q m i n = 1/Z/o = (^A/o), i-e. from the uncertainty principle follows that if the 
minimal momentum transfer is small, then the formation length is large. 

• For creation of the electron-positron pair the formation length is the distance re- 
quired for the final state particles to separate enough that they act as independent 
particles. 

• This is the distance over which the amplitudes of several interactions can add co- 
herently to the total probability. 

• The transverse size of the region, where the radiation is formed g± is determined 
by the minimal transverse momentum transfer q™ m : g± = \jq™ n = (h/q™ m ). 
Although usually qj_ ^> q» and this dimension is much smaller, we consider below 
the situation when the constraint of q± is the main effect. 

The formation length concept is important because within it the amplitudes can add 
coherently if a charged particle is moving freely. In the opposite case, when some pertur- 
bation is acting on the particle within the formation length, breakdown of the coherence 
happens and radiation will be suppressed. 
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1.2 Processes in a medium 



Landau and Pomeranchuk were the first to show that if the formation length of the 
bremsstrahlung becomes comparable to the distance over which the multiple scattering 
becomes important, the bremsstrahlung will be suppressed pQ. They considered radiation 
of soft photons. Migdal 0, [Sj developed a quantitative theory of this phenomenon. Now 
the common name is the Landau- Pomeranchuk -Migdal (LPM) effect. 

Let us estimate a disturbance of the emission process due to a multiple scattering. As 
it is known, the mean angle of the multiple scattering at some length lf is 

= = e a = mJ^=21.2MeV, (1.8) 
v £ V 2L rad V ol 

where a = e 2 = 1/137, L ra d is the radiation length. Since we are interesting in influence 
of the multiple scattering on the radiation process we put here the formation length (jl.fij) . 
One can expect that when d s > I/7 this influence will be substantial. From this inequality 
we have 

e' uj m^L. 



>-, e LP ='^p±, (1.9) 

here e^p is the characteristic energy scale, for which the multiple scattering will influence 
the radiation process for the whole spectrum. It was introduced in and denoted by 
Eq = Elp- For tungsten we have Elp = 2.7 TeV, and similar values for all the heavy 
elements. For the light elements the energy elp is much larger. When e <C Elp the 
multiple scattering will influence emission of soft photons only with energy 

E 2 

cj<u lp = , (1.10) 

elp 

e.g. for particles with energy e = 25 GeV one has ulp ~ 230 MeV. Of course, we estimate 
here the criterion only. For description of the effect one has to calculate the probability 
of the bremsstrahlung taking into account multiple scattering. 

Side by side with the multiple scattering of an emitting electron one has to take into 
account also the influence of the medium on the radiated electromagnetic field. As it was 
shown above Eq. (jl.5|) . this effect leads also to suppression of the soft photon emission 
(the dielectric suppression, the Ter-Mikaelian effect, see in 0). 

A clear qualitative analysis of different mechanisms of suppression is presented in early 
papers 0,0- More simple derivation of the Migdal's results is given in [Hj. 

The next step in a quantitative theory of the Landau-Pomeranchuk-Migdal (LPM) 
effect was made in This theory is based on the quasiclassical operator method in QED 
developed by the authors (see [7j, [S], [ID], [II])- One of the basic equations (obtained 
with the use of kinetic equations describing the motion of an electron in the medium in 
the presence of an external field) is the Schrodinger equation in an external field with 
imaginary potential (Eq.(3.3) of jU]; see (jA.22|) in Appendix A). 

New activity with the theory of the LPM effect is connected with a very successful 
series of experiments ^2] - QHj performed at SLAC recently. In these experiments the 
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cross section of the bremsstrahlung of soft photons with energy from 200 keV to 500 MeV 
from electrons with energy 8 GeV and 25 GeV was measured for a variety in materials with 
an accuracy of the order of a few percent. Both the LPM and the dielectric suppression 
was observed and investigated. These experiments were the challenge for the theory since 
in all the published before papers calculations are performed to logarithmic accuracy 
which is not enough for a description of the new experiments. The contribution of the 
Coulomb corrections (at least for heavy elements) is larger then experimental errors and 
these corrections should be taken into account. In papers [Hj-[2S] the problem was 
investigated using different approaches. 

Very recently the LPM was studied at CERN [16 j using iridium target and electrons 
with energy 287 GeV, 207 GeV and 149 GeV. For iridium according to Eqs. (jl.9|) and 
(ll.lOj) Elp=2.2 TeV, u LP ~ 37 GeV for the highest energy and the LPM suppression of 
bremsstrahlung was observed in much wider interval of photon energies than at SLAG 

The crossing process for the bremsstrahlung is the pair creation by a photon. The cre- 
ated particles in a medium can undergo the multiple scattering. It should be emphasized 
that for the bremsstrahlung the formation length increases strongly if u <C e. Just 
because of this the LPM effect was investigated at SLAC at a relatively low energy. For 
the pair creation the formation length 

i r = 2 -^ (i.n) 

m z u) 

attains maximum at e = uj/2 and this maximum is l P) max = (u>/2m)\ c . Because of this 
even for heavy elements the effect of multiple scattering on pair creation becomes notice- 
able at photon energies uj > 10 TeV [21]. Starting from these energies one has to take into 
account the influence of a medium on the pair creation and on the bremsstrahlung hard 
part of the spectrum in electromagnetic showers being created by the cosmic ray particles 
of the ultrahigh energies. These effects can be quite significant in the electromagnetic 
calorimeters operating in the detectors on the colliders in TeV range. 

The nonlinear effects of QED are due to the interaction of a photon with electron- 
positron field. These processes are the photon-photon scattering, the coherent photon 
scattering, the photon splitting into two photons, and the coalescence of two photons into 
photon in the Coulomb field. 

The process of creation of the electron-positron pair by a photon connected with the 
coherent photon scattering |26] . In the quasiclassical approximation the amplitude M of 
the coherent photon scattering is described by diagram where the electron-positron pair 
is created by the initial photon with 4-momentum k± (u, ki) and then annihilate into the 
final photon with 4-momentum k 2 . For high energy photon uj ^> m this process occurs 
over a rather long distance, known as the time of life of the virtual state 

'' = £■ (L12) 

where q c > m is the characteristic transverse momentum of the process. When the virtual 
electron (or positron) is moving in a medium it scatters on atoms. The mean square of 
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momentum transfer to the electron from a medium on the distance If is 

2 4 
Hmin ^ 



q 2 s =4nZ 2 a 2 n a L(qc)lf, L(q c ) = In Js_, g 2 ^ = a" 2 + A 2 + ?L, (1.13) 



where Z is the charge of nucleus, n a is the number density of atoms in the medium, A is 
the photon momentum transfer (A = |k2 — ki| <C q c ), cl is the screening radius of atom. 

The coherent photon scattering amplitude M can be obtained from general formulas 
for probabilities of electromagnetic processes in the frame of the quasiclassical operator 
method (see e.g. [ID])- It can be estimated as 

M~— 2| = q 2 . 1.14 

27ilfn a q£ 7in a 

We use the normalization condition ImM = ua p for the case A = 0, where a p is the total 
cross section of pair creation by a photon. 

In the case of small momentum transfer q s = ^Jq^ Cmwe have in the region of small 
A < m 



2Z 2 a 3 u . m 2 



q l = m\ M —In T . (1.15) 

At an ultrahigh energy it is possible that q s 3> m. In this case the characteristic 
momentum transfer is q c , the combination q 2 /m 2 enters in the formation length Lf fll.5|) . 
The self-consistency condition is 



9c 



where L(g c ) is defined in Eq. fll,13|) . So using Eq. (|1.14|) one gets for the estimate of the 
coherent photon scattering amplitude M (the influence of the multiple scattering manifests 
itself at the high photon energies such that m 2 a/uj 1) 



2Z 2 a 3 cu / Ai 



M ~^^f n a -2 + A 2 > K = p™Z^n a ^>A 2 . (1.17) 

1.3 Radiation in the presence of an external field 

The authors have shown j27j that an external electromagnetic field can also modify the 
bremsstrahlung process. We consider the effect of an external field on the photon emission 
in electron-electron(positron) collision. This effect is associated as well with a reduction 
of the formation length of a photon (either real or virtual) due to relatively large turn of 
particle velocity over this length, and the corresponding increase in the vertex angle of 
the radiation cone. If a photon of energy Q is emitted by an electron (positron) at an 
angle to its velocity, the formation length of such a photon will be according to Eq. (J1.5)) 
(we neglect the polarization of a medium) 

2(e - Q) 
em 2 euti&l 



*/(«) = = — (1-18) 
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where uq = Q/(e — Q), q c defined in Eq. fll.l6|) . The characteristic emission angle in weak 
field is $ c = I/7, and we can neglect the influence of the external field if 

. „ ^ N eHlAVL) 1 

vZ/(fi = ^«~. 1.19 

€ 7 

Substituting Eq. (jl.l8}) into Eq. (jl.l9j) . we have the criterion for a field to be weak 

7 !- = -± < 1. (1.20) 

Here we introduced the important parameter x, which characterizes the quantum effects 
in an external field, when x < 1 we are in the classical domain and with \ > 1 we are 
already well inside the quantum domain 

X = ^? X=\X\, F = E ± + (vxH), (1.21) 

-TO 

where E(H) is an electric (magnetic) field, Eq. (jl.21j) contains the transverse with respect 
to particle velocity v component of fields, F is one of 

Eo = ™L. = (^Lpj = 1.32 ■ 10 16 V/cm, or H = ^ = fe^jpj = 4.41 • 10 13 Oe. (1.22) 

In strong fields, where ^/u> 1, with characteristic radiation angles $ 3> I/7, the effective 
emission angle is determined by a self-consistency argument: the deviation angle of 
the particle in the field over the formation length must not exceed i? c , i.e. 



1 (2L) 1/3 = ( 2eH \ 1/3 

7 \u n J \e 2 uaJ 



7 VMo/ \e 2 unj 



muu \2xJ \4une 2 H 2 

It can be seen from Eq. (jl.2Hj) that when ^/m > 1, neither the characteristic emission 
angle nor the photon formation length depends on the mass of the radiating particle. 
A parameter characterizing the effect of an external field on the radiation process was 
derived in 0. 

It follows from Eq. (jl.23|) that, for x/ u n ^ 1, the characteristic photon formation 
length is reduced by factor (x/w^) 2 / 3 , and the emission angle is increased by factor 
(x/un) 1 ^ 3 - Since the relevant parameter x/ u n depends on the frequency £1 (via u n ), 
this effect is manifested earlier for soft photons. Because of this a specific situation arises 
in bremsstrahlung in a collision of an electron and a positron(an electron) in colliding- 
beam experiment. The point is that the external factors act differently on the radiating 
particle and on the recoil particle. For the radiating particle the criterion of influence of 
external factors is the same in scattering by a nucleus as in collision of particles. For the 
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recoil particles the effect turns out to be enhanced by 7 2 times (in the laboratory frame). 
This is due to the fact that the main contribution to the bremsstrahlung cross section 
give the virtual photons emitted by recoil particle with very low energy 

(1*) 

where u is the frequency of the real photon. The formation length of the virtual photon 
is 

h = *,(*) = (1-25) 



for example, for e = 2 GeV and u — 10 MeV the formation length is enormous: l v 
4 • 10 3 cm. This is the reason why the field has a significant effect on virtual photon 
emission in e~e + collisions even for moderate field and at relatively low particle energies, 

when the parameter — = — is small while the value — — ~ 47 2 — is large. In 

oj u qo u 

that event, the picture of emission of a bremsstrahlung photon (radiation vertex) does 
not change directly, but a significant change takes place in the virtual photon spectrum 
at momentum transfers \q\ < g m j n (47 2 x/-u) 1 / 3 , increasing the lower bound of the effective 
momentum transfer and resulting in a corresponding decrease of the cross section. The 
bremsstrahlung cross section under these conditions was derived to logarithmic accuracy 
in |2Zj using the equivalent photon method. In |2Hj the problem was solved to relativistic 
accuracy (terms 0(1/7) were neglected). 

It follows from Eq. (|1.2Uj) . and was demonstrated in jU] that in the case x/ u > 1 there 
is a change of the radiation vertex describing absorption of a virtual photon and emission 
of a real one. Since for \/u ^> 1 the formation length of a real photon falls off as 
(m/x) 2 / 3 (see Eq. (jl.23|0 . the bremsstrahlung cross section falls off in just the same extent. 



2 The effect of a medium in an infinitely thick target 
2.1 The LPM effect for radiation 

We consider first the case where the formation length is much shorter than the thickness 
of a target 1(1 f <C /). The derivation of basic formulas for probability of radiation under 
influence of multiple scattering is given in Appendix A. The spectral distribution of the 
probability of radiation per unit time in this case is given by Eq. (|A.20|) and functions 
(Pn((p ,<p) satisfy the Eq. (|A.22|) . Introducing new variables 




where the quantities a and b are defined in Eqs. ()A.14|) . (|A.15|) . we obtain for the spectral 
distribution of the probability of radiation per unit time 

00 

-j^ = ^Re J dte~ u [Rm(0, t) + i? 2 p^(0, t)} , (2.2) 
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u 2 e e' 

where Ri = — -, R2 = — H , and the functions ip^ = f^(fo, V 5 ) now satisfy the equa- 
tion 



H<p m H = p 2 -iV(g) 1 p = -tV g , V(g) = -Qg 2 (\n 1 2 i9 2 1 



dt 

g 2 _„ ^ „ 2nn n Z 2 a 2 ee' 
4 ; " m 4 uj 



+ ln^- + 2C-l), Q= a 4 , C = 0.577216... (2.3) 



with the initial conditions ip (g,0) = 5(g), ip(g,0) = p5(g). Here e,e',uj are defined 
in Eq. (jl.2|) and Z, n a are defined in Eq. (|1.13jl . The functions ip and ip in Eq. (j2.2|) are 
rescaled according with the initial conditions (factors I/7 2 and I/7 3 , correspondingly). 
Note, that it is implied that in formula (|2.2j) subtraction at V = is made. 

The potential V(g) Eq. (j2.3|) corresponds to consideration of scattering in the Born 
approximation. The difference of exact (as a function of Za) potential V(g) and taken 
in the Born approximation is computed in Appendix A of ^7]. The potential V(g) with 
the Coulomb corrections taken into account is 



V(g) = -Qf> 2 (ln 7 2 tf 2 + ln^- + 2C-l + 2/(Zc 

= -Qg 2 (\n 7 2 tf 2 + l n ^ + 2C) , (2.4) 

where $2 = $1 exp(/ — 1/2), the function / = f(Za) is 

00 1 

/(o = Re mi + ^ - m = e e n{n 2 +ey ( 2 - 5 ) 

where i{>(£) is the logarithmic derivative of the gamma function. 

In above formulas g is space of the impact parameters measured in the Compton wave- 
lengths A c , which is conjugate to space of the transverse momentum transfers measured 
in the electron mass m. 

An operator form of a solution of Eq. (J2.3)) is 

<A)(f? ; t) = exp(—iHt)(p (g, 0) =< g\ exp(—iHt)\0 >, H = p 2 — iV(g), 

cp(g, t) = exp(—iHt)p(f (g, 0) =< g\ exp(— iHt)p\0 >, (2.6) 

where we introduce the following Dirac state vectors: | g > is the state vector of coordinate 
g, and < £>|0 >= 5(g). Substituting Eq. ()2.6|) into Eq. (j2.2|) and taking integral over t we 
obtain for the spectral distribution of the probability of radiation 

^ = ^ImT, T=<0| J R 1 (G- 1 -G - 1 ) + J R 2 p(G- 1 -G - 1 )p|0>, (2.7) 

where 

G = p 2 + 1-^, G = p 2 + 1. (2.8) 
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Here and below we consider an expression < 1 1 > as a limit: lim x — > 0, 

lim x' — > of < x|...|x' >. 

Let us estimate the effective impact parameters q c which give the main contribution 
into the radiation probability. Since the characteristic values of g c can be found straight- 
forwardly by calculation of Eq. (|2.7|) . we the estimate characteristic angles d c connected 
with q c by an equality g c = l/(7$ c ). The mean square scattering angle of a particle on 
the formation length of a photon If Eq. ()1.5|) has the form 

®s = —^—nalf In ^ = — In (2.9) 

where ( = 1 + 7 2 i9 2 , we neglect here the polarization of a medium. When $1 <C I/7 2 
the contribution in the probability of radiation gives a region where C, ~ l(i9 c = I/7), in 
this case £ c = 1. When d s ^> I/7 the characteristic angle of radiation is determined by 
self-consistency arguments (compare with Eq. (|1.23jl ): 

Vs~V c - 72 in ^ , Q to ^ 1, 4G^ C to 1. 

It should be noted that when the characteristic impact parameter q c becomes smaller than 
the radius of nucleus R n , the potential V(g) acquires an oscillator form (see Appendix B, 
Eq.(B.3) in [T?]) 

V{q) = Qg 2 (hi || - 0.041^ , a s = OMclbZ- 1 / 3 , (2.11) 

here a s is the screening radius, ag is the Bohr radius. 

Allowing for the estimates Eq. (|2.1()j l we present the potential V(g) Eq. (|2.3jl in the 
following form 

V(q) = V c (q) + v(q), V c (g)=qg 2 , q = QL C , L c = L(g c )=\n ' 



^1)=]^ = !^, (2^ fcg), (2.12) 

here a s2 = a s exp(— / + 1/2). So, one can to redefine the parameters a s and #1 to include 
the Coulomb corrections. The inclusion of the Coulomb corrections (f(Za) and -1) into 
In $2 diminishes effectively the correction v(g) to the potential V c (g). In accordance with 
such division of the potential we present the propagators in Eq. (|2.7|) as 

G- 1 - G 1 = G- 1 - G- 1 + G- 1 - G^ (2.13) 

where 

G c = p 2 + 1 - iV c , G = p 2 + 1 - iV e - iv 

This representation of the propagator G~ l permits one to expand it over the "perturba- 
tion" v. Indeed, with an increase of q the relative value of the perturbation is diminished 



11 



V 1 

(— ~ — ) since the effective impact parameters diminishes and, correspondingly, the 

c c 

value of logarithm L c in Eq. ()2.12|) increases. The maximal value of L c is determined by 

the size of a nucleus R„ 



In J| ~ 2 In = 2L 1; L x = 2(ln(183Z- 1 /3) _ f ( Za )), 



[2.U) 



The value Li is the important parameter of the radiation theory. 

The matrix elements of the operator G~ x can be calculated explicitly. The exponential 
parametrization of the propagator is 



(X) 

G- 1 = ij dte- u exp(-iH c t), H c = p 2 - 



iqg 



(2.15) 



The matrix elements of the operator exp(— iH c t) has the form (details see in [T7] ) 



< ei|exp(-zW c t)|e 2 >= K c (g 1 ,g 2 ,t), 



K c {g 1 ,Q 2 ,t) 



Ani sinh i^t 



exp < — 



^i + 62) cothz/t — 



sinh vt 



Q1Q2 



, (2.16) 



where v = l^fiq (see Eq. ()2.12|) ). 

Substituting formulas (|2.15jl and (j2.1(i|l in the expression for the spectral distribution 
of the probability of radiation ()2.7j) we have 



a. 



duo 2%Y 



$(V) —v \ die 



■Im $(V) 



Ri 



- - ) - ivR 2 



sinh z z 



1 



R\ I In p — ip I p 



+ R 2 lip (p) — hip + 



sinh z z 2 
1 \ 



2p 



(2.17) 



where z — vt, p — i/(2v), let us remind that ip(x) is the logarithmic derivative of 
the gamma function (see Eq. (J2.5j) ). If we substitute in Eq. (j2.17j) t — > texp(— m/A) the 
probability will be transformed into the form containing the real functions only. Then it 
can be written as 

<nV ' n ' H ViG( S ) +2R 2( f>( S )], 



where 



duo 



G(s) = 12s 5 



6s" 



12 7 2 

Ti 7 sin(sx) 

/ exp(— sx) dx 

4 J smh x 



o 

J exp(— sx) sin(sx) cothxrfx — — 
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(2.18) 



(2.19) 



are the functions introduced by Migdal (see Eqs.(46), (47) in [2]), note that parameter 
s = — r in Eq. (j2.19)) is two times larger than that used by Migdal. If in Eqs. (j2.18|) . 



V2\ 



(12.19)1 one omits the Coulomb correction, then the probability 1)2.18)1 coincides formally 
with the probability calculated by Migdal (see Eq.(49) in [2]). 
Expanding the expression G^ 1 — G^ 1 over powers of v we have 



G 



-i 



G 



-i 



G-'MG- 1 + G^^G-'MG- 1 + ... 



[2.20) 



In accordance with Eqs. (J2.13)l and (J2.20)) we present the probability of radiation in 
the form 

dW dW c dW 1 dW 2 , nM . 

+ ^ + ^ + ... (2.21) 



duo 



duo 



duo 



duo 



1 



At Q > 1 the expansion (J2.2U)) is a series over powers of — . It is important that variation 

L 

< >f the parameter q c by a factor order of 1 has an influence on the dropped terms in ()2.2())1 

dW c . . 

only. The probability of radiation — — is defined by Eq. ()2.17j) . 



The term 



dW 1 
duo 



duo 

in Eq. ()2.21j) corresponds to the term linear in v in Eq. (J2.20j) . The 



explicit formula for the first correction to the probability of radiation ^7j is 



dW l 
duo 

£>i0o) 



a 



47T7 2 L, 



-ImF(u); -1m F{u) = D 1 {u )R 1 + -D 2 {u )R 2 ; s 



dze 



sinh z 



TV 



d(z) sin sz + -^g{ z ) cos sz 



D 2 (vo) = J 



dze 



sinh z 



x 



d(z) 
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(sin sz + cos sz) + —g[z) (cos sz — sin sz) 



d(z) = (In 1/0$ (1 — u ) — In sinh z — C)g(z) — 2G(z) cosh^, 



[2.22) 



where 



z 

g(z) = 2; cosh 2; — sinhz, G(z) = J (1 — y cothy)dy 



-z---^r-z\n(l- e~ 2z ) + -Li 2 (e~ 2z 



2 12 

here Li 2 (x) is the Euler dilogarithm; and 



(2.23) 



9 ,19 \ 87m a Z 2 a 2 se' 
v 2 = \u\ 2 =Aq = AQL(q c ) = a —_ L(g c ), 



m A uo 



[2.24) 



The functions G(s) and <p(s) are given in Fig. 1 and the functions Di(i> Q ) and D 2 (i> Q ) are 
shown in Fig. 2. 
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As it was said above (see Eqs. (j2.10j) . (|2.14j) ). g c — 1 at 



= u f = 4QLx < 1. (2.25) 

The logarithmic functions L c = L(g c ) and Li are defined in Eqs. (j2.12|) and ()2.14|) . If the 
parameter \v\ > 1, the value of g c is defined from the Eq. ()2.10|) . where $i — > up to 
g c = R n /X c . So, we have two representation of \u\ depending on g c : at g c — 1 it is \v\ — v\ 
and at g c < 1 it is \v\ — Vq. The mentioned parameters can be presented in the following 
form 

2 -2 2 | ,2 2 A . „/ 2 ^ 1 - X CJ 

y Li i J e e X £ 

(SixZ 2 a 2 n a \lLi) , L c ~ L : M + ^^(^i " 1) ] (2.26) 



m 



It should be noted that in the logarithmic approximation the parameter g c entering into 
the parameter v is defined up to the factor ~ 1. However, we calculated the next term 
of the decomposition over v (q) (an accuracy up to the "next to leading logarithm") and 
this permitted to obtain the result which is independent of the parameter g c (in terms 
oc 1/L). Our definition of the parameter g c minimizes corrections to Eq. (|2.17|) practically 
for all values of the parameter g c . 

The approximate solution of Eq. (j2.10J) given in Eq. (j2.26|) has quite good numerical 
accuracy:it is ~ 2 % at v\ = 100 and ~ 4.5 % at v% — 1000. The LPM effect manifests 
itself when 

Vl {x c ) = 1, x c = ^j = — ^— . (2.27) 

So, the characteristic energy e e is the energy starting from which the multiple scatter- 
ing distorts the whole spectrum of radiation including its hard part. If the radiation 
length L ra d is taken within the logarithmic approximation the value e e coincides with Elp 
Eq. (jl.9|) . The difference between these values is 

1 <2% 



e e 9^! 

The formulas derived in and written down above are valid for any energy. 

In Fig. 3 the spectral radiation intensity in gold (e e = 2.5 TeV) is shown for different 
energies of the initial electron. In the case when e <C e e (e = 25 GeV and e = 250 GeV) 
the LPM suppression is seen in the soft part of the spectrum only for x < x c ~ ej ' e e 1 
while in the region e > e e (e = 2.5 TeV and e = 25 TeV) where x c ~ 1 the LPM effect 
is significant for any x. For relatively low energies e = 25 GeV and e = 8 GeV used in 
famous SLAC experiment J7j, ^3] we have analyzed the soft part of spectrum, including 
all the accompanying effects: the boundary photon emission, the multiphoton radiation 
and influence of the polarization of the medium (see below). The perfect agreement of the 
theory and data was achieved in the whole interval of measured photon energies (200 keV< 
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uj <500 MeV), see below and the corresponding figures in JTJ , JE] 5 • It should be 
pointed out that both the correction term with F(u) and the Coulomb corrections have 
to be taken into account for this agreement. 

When a scattering is weak [v\ <C 1), the main contribution in Eq. (J2.22j) gives a region 
where x<l. Then 



-Im F{u) = -Im v 2 (R 2 - R t ) , L c -> L x , 

<l>(v)~^-{R 1 + 2R 2 ), (H<1) 
o 



(2.28) 



Combining the results obtained in Eq. (j2.28|) we obtain the spectral distribution of the 
probability of radiation in the case when scattering is weak <C 1) 



dW dW, 



cL; do; 
a 2Q 

~ 2vt7 2 T" 
4Z 2 a 3 n Q 
3m 2 w 



+ 



a. 



du 



2vt7 2 
1 



Im 



1 



2Lx 



F{v) 



2R 2 L x 



to 



ln(l83Z" 1/3 



1 

6 



1 

f(Za 



2 i 1 + iP) ( ln ( 183Z " /3 ) + T2 ~ f{Za] 



(2.29) 



where L\ is defined in Eq. (l2.14j) . This expression coincide with the known Bethe-Maximon 
formula for the probability of bremsstrahlung from high-energy electrons in the case of 
complete screening (if one neglects the contribution of atomic electrons) written down 

within power accuracy (omitted terms are of the order of powers of — ) with the Coulomb 

corrections, see e.g. Eq. (18.30) in or Eq.(3.83) in [37]. 
At v Q > 1 the function Im F(u) Eq. ([2~!?2j) has the form 



7 



-Im F(v) 
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In2-C 



7T 



4 



R2 



Under the same conditions (z/ 1) the function Im $(f) Eq. (j2.17|) is 

Im Mv) = i Rl + - R2 _ « R2 . 



(2.30) 



:2.3V 



So, in the region where the LPM effect is strong the probability Eq. (j2.17j) can written as 



dW r aRo Z 2 a 2 ee'n 



1/2 



du 



7TLO 



(2.32) 



This means that in this limit the emission probability is proportional to the square root 
of the density. This fact was pointed out by Migdal [2j, the expression ()2.32j) coincides 
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formally with the probability calculated by Migdal (see Eq.(52) in 2\) with the same 
remarks which were made after Eq. (j2.19|) . 

dWi 

Thus, at i/q ^> 1 the relative contribution of the first correction — - — is defined by 

du 



dWi 1 
dW c ~ 2L(g c ) 



In 2 



c + i 



0.451 

Hoc)'- 



(2.33) 



where L(g c ) = In 



'•«2 



\2 n 2 
A C Q C 



In this expression the value r with the accuracy up to terms 



~ 1/L 2 . doesn't depend on the photon energy:L c ~ L\ + ln(e/e e )/2. Hence we can find 
the correction to the total probability at e ^> e e . The maximal value of the correction is 
attained at e ~ 10e e , it is ~ 6% for heavy elements. 

In the above analysis we did not consider an inelastic scattering of a projectile on 
atomic electrons. The potential V e (g) connected with this process can be found from 
formula (j2.4j) by substitution Z 2 — > Z, $i — ► i9 e = 0.153$i (an analysis of an inelastic 
scattering on atomic electrons as well as the parameter d e can be found in [HI]). The 
summary potential including both an elastic and an inelastic scattering is 



V(q) + V e (g) = -Q(l + I) e 2 [ln 7 2 ^ + In ^ + 2C + ' 



-Q efe 2 (W$l f + ln°- + 2C) 



z+i r*f 



V 



[2.U) 



where 



Q ef = Q(l + — ), & ef = exp 



—=(Zf(aZ) - 1.88) - - 



2.2 Integral characteristics of bremsstrahlung 

Integrating Eq. (j2.17j) over x = uj/s we obtain the total intensity of radiation in the 
logarithmic approximation 



'-if 



rad 



2^Im 



i 



dx / 



X 



[2(l-x)+x 2 



+ 



/ (^( p +l ) ~^(p + \)) +2 j xdx (if) (jp + 1) - hip) 



(2.35) 



where 



V 



9V 
2 ' 



x 



x 



g = exp it 
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'Ll£e 
L r € 



L® ad is the radiation length in the logarithmic approximation. The relative energy loss 
of electron per unit time in terms of the Bethe-Maximon radiation length L® ad : —L° rad 
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in gold is given in Fig. 4 (curve 1), it reduces by 10% (15% and 25%) at e ~ 700 GeV 
(e ~ 1.4 TeV and e ~ 3.8 TeV) respectively, and it cuts in half at to ~ 26 TeV. This 
increase of effective radiation length can be important in electromagnetic calorimeters 
operating in detectors on colliders in TeV range as well as in analysis of high-energy 
cosmic rays. The relative energy loss of electron per unit time in terms of the exact 
Bethe-Maximon radiation length L ra d in iridium and lead calculated taking into account 
the correction terms Eq. (j2.22J) is shown in Fig. 19 (see Sec. 3. 10.2). 

The spectral distribution of bremsstrahlung intensity and the spectral distribution 
over energy of created electron (positron) as well as the reduction of energy loss and the 
photon conversion cross section was calculated by Klein ^3] using the Migdal [2] formulas. 
As will be explained below (see Sec. 3.7) we use more accurate procedure of fine tuning 
and because of this our calculation in logarithmic approximation differs from Migdal one. 
We calculated also the correction term and include the Coulomb corrections. For this 
reason the results shown here in Figs. 3-4 are more precise than given in |15j . 

In the case e <C e e the correction term to the total intensity of radiation (or to the 
radiation length) can be written as [21] 

I am 2 / 1 4tt e \ T -\ ( ^ e 

1 = 4^ \ 1 + 9T 1 ~ TE7 e ) ~ rad v ~ TE7 t 

1 1 ^'), * = 2gqya (2 . 36) 



Lrad L° rad \ 9L l J ' LP rad m 2 
In the opposite case e ^> e e we have [2T] 

1 - 2.37a/ — + 4.57— + — In 0.3455 

V e e ALi \ e e ) 



(2.37) 



£L ra( i 2 V e 

Although the coefficients in the last expression are rather large at two first terms of the 



decomposition over \je e je this formula has the accuracy of the order of 10% at e ~ 10£ e . 
Analogous expression for the integral probability of radiation has the form 



llirZ 2 a s n a fe 



w = — = — 

2y/2m 2 



1 - 1.23 J— + 1.645— + — In — + 2.53 

£ e 4Li V e e 



'2.38) 



Ratio of the main terms of Eqs. ()2.37J) and (|2.38j) gives the mean energy of radiated 
photon 



2.3 The LPM effect for pair creation 

The probability of the pair creation by a photon can be obtained from the probability of 
the bremsstrahlung with help of the substitution law: 

u 2 du> — > e 2 de, to —u, e — > —e, (2.40) 
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where u is the photon energy, e is the energy of the particle. Making this substitution in 
Eq. ()2.7|) we obtain the spectral distribution of the pair creation probability [2T] (over the 
energy of the created electron e) 

dW p _ 2auf lm , |si _ ^ + _ p|Q v ^ (2 41) 



ee' 

where 

si = l, 3 2 = - ' ° , e' = w-e, (2.42) 



e 2 + e' 2 



w 2 



the other notations are given in Eos. ([27%]) . (jQj) , (1273)) . (1275)1 . 

Just as in the analysis of the LPM effect for radiation we present the potential Eq. ()2.4j) 
in the form V(g) = V c (g) +v(g) and will proceed as we did in Eqs. (j2.12|) - ()2.16)) . Substi- 
tuting the expression ()2.16)) in the formula for the spectral distribution of the pair creation 
probability (|2.41j) we have 



am 2 



(l\V;, ,.„. 

— = :1m %(u) 



de 2nee' 

oo 

= v j dte~ u 
o 



si -IVSi 



sinh z z) * V sinh 2 z z 2 



= s x (inp - $ [p + + s 2 U{p) - Inp + , (2.43) 

where we used the same notations as in Eq. (|2.17)l . This formula gives the spectral distri- 
bution of the pair creation probability in the logarithmic approximation which was used 
also by Migdal [3]. It should be noted that the parameter g c entering into the parameter v 
(see Eqs. (l2.16|) . (|2.12|) and (j27TUJ)) is defined up to the factor ~ 1, what is inherent in the 
logarithmic approximation. However, below we will calculate the next term of the decom- 
position over v(g) (an accuracy up to the "next to leading logarithm") and this permits 
to obtain the result which is independent of the parameter g c . It will be shown that the 
definition of the parameter g c minimizes corrections to Eq. ()2.43j) practically for all values 
of the parameter g c . It should be emphasized also that here the Coulomb corrections are 
included into the parameter v in contrast to j3J. 

Substituting the expansion (|2.20)) in Eq. (j2.41j) we obtain the decomposition of the 
probability of the pair creation: 

dW p dW; dWl dW 2 

—r 1 - = —r^ + —r^ + —r^ + ... 2.44 

de de de de 

dW v dW} 

The probability of pair creation — — - is defined by Eq. (j2. 43)1 .The term — in Eq. (j2.44j) 

de Mn-rrm. 

corresponds to the first term (linear in v) in Eq. ()2.20j) . The explicit expression for the 



first correction to the pair creation probability is [21] 

dWl am 2 „ „ , s „ , s °r dze 



tw„ am* „ . . ^ . . r dze " „ , . „ , . 
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fx (z) = ( In q c + In — — In sinh z — C ) g(z) — 2 cosh zG(z) 



h(z) 



sinh z 



t = tx+t 2 



(2.45) 



where the functions and £7(2) are defined in Eq. ()2.23|) . L c is defined in Eq. fl2.12|) . 
Use of the last representation of function G(z) simplifies the numerical calculation. As in 
radiation case the parameter v has two representations: Eqs. (j2.24|) and (J2.25|) . 

When the scattering of created particles is weak {yx <C 1), the main contribution in 
Eqs. fl2.43|) and ()2.45|) gives the region where z <C 1 and one can expand the functions 
of z (compare with Eq. (|2.28jl ). As result we obtain the spectral distribution of the pair 
creation probability for \v\ <C 1: 



dW n 



de de 
4Z 2 a 3 n, 



dW c p dwl 



am 



3m 2 uj 



de 2nee' 
In ( 183Z- 1/3 ^ 



Im 



2L C 



+2 



e 2 + e' 2 



f(Za 

(in (l83Z~ 1/3 ) - f(Za 



31uf 
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H 


16z/f\ 
21 j 


i 

+ I2 



(2.46) 



where Li is defined in Eq. (l2.12j) . Integrating Eq. (|2.46p over e we obtain 



28Z 2 a i n a 
9m 2 



(in^Z- 1 / 3 ) - f(Za] 



3312 J 1 " 
2401^2 



1 

42 



where 



m 



(2TTZ 2 a 2 n a X 3 r Lx 



{2A7) 



'2.48) 



Note that u e is four times larger than e e Eq. (j2.26|) . in gold it is u e = 10.5 TeV. This is 
just the value of photon energy starting with the LPM effect becomes essential for the 
pair creation process in heavy elements. If one omits here the terms cx v\ and oc (u/uj e ) 2 
these expressions coincide with the known Bethe-Maximon formula for the probability of 
pair creation by a high-energy photon in the case of complete screening (if one neglects 
the contribution of atomic electrons) written down within power accuracy (omitted terms 

are of the order of powers of — ) with the Coulomb corrections, see e.g. Eqs.(19.4) and 

(19.17) in 0. 

The pair creation spectral probability dW/dx vs x = e/uo in gold is shown in Fig. 5 for 
different energies. It is seen that for u = 2.5 TeV which below u e the difference with the 
Bethe-Maximon probability is rather small. When u > u e there is significant difference 
with the Bethe-Maximon spectrum increasing with u growth. In Fig. 5 are shown the 
curves (thin lines 2,3,4) obtained in logarithmic approximation dW^/de Eq. (j2.43|) . the 
first correction to the spectral probability dW^/de Eq. (j2.45|) . curves c2,c3,c4 and the 
sum of these two contributions: curves T1,T2,T3,T4. It should be noted that for our 
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definition of the parameter g c Eq. (j2.1U|) the corrections are not exceed 6% of the main 
term. The corrections are maximal for z/ ~ 3. 

Integrating Eq. ()2.43|) over y = e/uj we obtain the total probability of pair creation in 
the logarithmic approximation 



9 uj £ 

1m 

14 uj 



dy 



+ {l-2y + 2y< 



y(i -y) 

if) (p) — In p + 



\np 

1 " 



ip(p+- 



2 P/ 



(2.49) 



where 



bs 



1/(1-2/) 



7T 

(>xp |l- 



L fi a; 



W 7 ^ is the Bethe-Maximon probability of pair photoproduction in the logarithmic ap- 
proximation. The total probability of pair creation W£ in gold is given in Fig. 4 (curve 
2), it reduced by 10% at uj ~ 9 TeV and it cuts in half at uj ~ 130 TeV. 
At uq ^> 1 the function F p (v) (see Eq. (|2.45j l) has the form 



-Im F p (y) 
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vSl - s 2 ) + 



^0 / 

v/2 V 



In2-C + 



7T 



«2, 



(2.50) 



Under the same conditions (uq ^> 1) the function Im Q p (v) Eq. ()2.43|) is 

Im $» = - s 2 ) + ^=s 2 -> i| S2 . (2.51) 

So, in the region where the LPM effect is strong the probability Eq. (j2.43j) can written as 



aw; 

de 



a(e 2 + e' 2 ) 



ijj- 



Z 2 a z n a 



1/2 



Tiuoee 



(2.52) 



This means that in this limit the spectral probability of pair creation is proportional to 
the square root of the density just as the emission probability Eq. (j2.32| ). 

dW 1 

At v o 3> 1 the relative contribution of the first correction — — ^- is defined by 



dwl 

dW; 



1 



2L r 
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In2-C + - 
4 



de 
0.451 



(2.53) 



In this expression the value r with the accuracy up to terms ~ 1/L 2 . doesn't depend on 
the electron energy: L c ~ Li + \n(u/u e )/2. Hence we can find the correction to the total 
probability at uj ^> uj e . The maximal value of the correction is attained at uj ~ 10uj £ , it is 
~ 6% for heavy elements. In this region i/ > 1 or w > w e the total probability of pair 
creation in terms of the Bethe-Maximon total probability can be presented as 



W BM 



2.id- 

V uj 



[ZJ7 uj e 1 

1 - 0.836,/— - 0.548— + — 
I w uj 4Li 



UJ 



In — + 0.274 



UJp 



(2.54) 
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2.4 Anomalous magnetic moment of the electron in a medium 

The contributions of higher orders of the perturbation theory over the interaction with 
an electromagnetic field give the electromagnetic radiative corrections to the electron 
mass and lead to appearance of the anomalous magnetic moment (AMM) of the electron 
[2H|. It is known that under influence of an external electromagnetic field these effects in 
particular the AMM of electron are changed essentially |2|, J5J. Here we consider how 
the mentioned effects modify in a medium *.TQ • 

The total probability of radiation W is connected with imaginary part of the radiative 
correction to the electron mass according to 

mAm = eAe, -2ImA£ = W, Im Am = — —W. (2.55) 

2m 

Since T in Eq. (j2.7|) is the analytic function of the potential V we have that 



Am = —am 




(2.56) 



o 



One can derive this formula considering the self-energy diagram and the corresponding 
amplitude of forward scattering of electron (see Eqs.(12.18) and (12.39) in [H])- 

In the presence of a homogeneous external field the Hamiltonian 7i (j2.3|) . (j2.6J) ac- 
quires the linear over coordinate term. One can find the explicit form of this term using 
the Eqs. (jA.18|) and ()A.22|) of Appendix A and carrying out the scale transformation of 
variables according to Eq. (j2.1|) : 

AH = - X Q, (2.57) 

u 

where x is the known parameter characterizing quantum effects in the radiation process 
(see Ea. ljPT|) ). u = u/ef. 

For polarized initial electron (see (jA.4|) in Appendix A) have for the addition to T in 
Eq. (j2.7|) depending on the spin vector £: 

T —>T + T^, T^ = ^(<0|(G- 1 p-pG- 1 )|0> xv)C- (2.58) 
For the radiative correction to the electron mass we have respectively 

Am -> AM = Am + Am^. (2.59) 

It should be noted that expressions for T ()2.7)1 and ()2.58|) have universal form while the 
specific of the particle motion is contained in the propagator G^ 1 through the effective 
potential 

V f (q) = -iV c (q) + -xq, H F = p 2 + V F , Gp = Hp + 1, (2.60) 
u 

where we restrict ourselves to the main term in the decomposition over v, see Eqs. (j2.12|) . 
(|2.2(J|) . This means that result has the logarithmic accuracy over the scattering (but not 
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over an external field). With regard for an external field the parameter g c in Eq. (j2.24|) is 
defined by a set of equations: 



g c =l for 4^+4QL 1 <l 



u 



X 2 o 



4— Q c + \v{Qc) 



a 



1 for 4^r + 4QL 1 >l, 



if 



(2.61; 



where z/ = 2i/2g, g = QL C , L±, Q and L c are defined in Eqs. (j2.25|) and (|2.24|) . Substituting 
into the depending on the spin term in (J2.58|) we find 



OKG^p-pG^lO 

<p = <p(x, M) = ex P 



X 

2nu 



exp(-it) 



o 



cosh 2 f dtl 



^ H ^1-ltanh^ 
vt 2 



V 2 U 2 



Thus, in the used approximation {G = Gp) we have for Eq. (j2.58ft 



(2.62) 



T, 



C = 2 



itu e 



exp(-it) 



V 3 



cosh 2 f 



dt(Cxv) 



(2.63) 



where (Cx v ) — (C' {x x v ))- Within relativistic accuracy (up to terms of higher order over 
I/7) the combination of vectors entering in Eq. ()2.63|) can be written as (see Eq.(12.19) 
in [H] and [33]) 



(Cxv) = 2^o 



CH 



R 



m 



(2.64) 



where H^j = 7(Hj_ + E x v) is the magnetic field in the electron rest frame. Here we use 
that x = X s 1 where s is the unit vector in the direction of acceleration (this vector is used 
ini). 

In the electron rest frame one can consider the value Re Am^ depending on the 
electron spin as the energy of interaction of the AMM of electron with the magnetic field 
H R (see Eq. (12.23) in and [33]) 



Re Am, 



yen. 



;2.65) 



Taking into account Eas. (l2~K6l . (I23K1) . (J23HJ), (l2~E3l . (l2~H4l . (gHH) we obtain the 
following general expression for the AMM of electron moving in a medium in the presence 
of an external electromagnetic field |31j : 



E. 



a 



Im 



00 00 
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2 vt 
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(2.66) 
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In the absence of scattering [v — > 0) the expression ()2.66|) gets over into the formula 
for the AMM of electron in external field (see Eq. (12.24) in f\ and P]^ 



fi a 
— = 1m 

(JLq IX 



du 



exp 



-it 1 



3m 2 



dt 



(2.67) 



In the weak external field (x <C I,*/ 9 — 1) we obtain the formula for the AMM of 
electron under influence of multiple scattering 



— — Im 
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[1 + uf 
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dt 
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(2.68) 



where 



ReJ, J = 2i 



du 



l + u) 



exp (—it) 



cosh 



2 i/t 



dt 



1/ 





^4 



;i + M ) 3 



2? 



ft 



(2.69) 



where ^(x) is defined in Eq. (l2.5jl . 

The dependence of the AMM of electron on its energy e in gold is shown in Fig. 6. It 
is seen that at energy e = 500 GeV value of AMM is 0.85 part of standard quantity (SQ) 
of AMM (r = 1) , at energy e — 1 TeV it is 0.77 part of SQ and at energy e = 5.5 TeV 
it is 0.5 part of SQ. Actually in all heavy elements the behavior of AMM of electron will 
be quite similar, i.e. the scale of energy where the AMM of electron deviates from SQ is 
of order of TeV. 

In the case of weak effect of multiple scattering when £<e e , where e e is defined in 
Eq.flOHD we obtain for AMM of electron 
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At very high energy e ^> e e the effect of multiple scattering becomes strong. In this 
case the asymptotic expression for AMM of electron is 
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(2.71) 



Redefinition of the characteristic energy e e e c is connected with enlargement of the 
radiation cone (in comparison with I/7) or in another terms with increasing of the char- 
acteristic transverse momentum transfers due to the multiple scattering. 
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Let us discuss one of possibilities of experimental observation of influence of medium 
on the AM M of electron. At a very high energy where the observation becomes feasible the 
rotation angle <p of the spin vector £ in the transverse to the particle velocity v magnetic 
field H depends only on the value of AMM and doesn't depend on the particle energy: 

( m u'\ eH a H I . . 

<P = - + - — 1 - ^Fr. ( 2 - 72 ) 



where I is the path of electron in the field, H = m 2 /e = 4.41 • 10 13 Oe (see Eq. fll.22|) ). 
The dependence of r on energy e is found above and shown in Fig. 6. 

Since at radiation of hard photons in a medium the picture is quite complicated:energy 
losses, cascade processes, spin flip and depolarization, it is desirable to measure the par- 
ticles which don't radiate photons on the path I. The number of such particles N is 
determined by the total probability of radiation in a medium found in [2*T] . Eqs.(3.12)- 
(3.14) (see Eas. (l2~3S|) and (J2~33|) in Sec.2.2 above): 

iV = iV exp(-^)), ^e) = W{e)l=^^- ] k(e) = W(s)L° rad , 

(O-^f^, = (2.73) 

47r e e m Ho 

where ip$Q is the rotation angle for standard value of AMM in QED (r = 1), N is the 
number of initial electrons. With energy increase the function k(e) decreases (see Eq.(43) 
in EH) 

k(e<^e e ) ~ - fin— + 1.96V He = e e ) ~ 3.56, k(e > e e ) ~ ^=\t^- (2.74) 
3 V e J 4V2 V e 

The crucial part of the experiment is an accuracy of measurement of electron polarization 
before target and after target. If one supposes that spin rotation angle can be measured 
with accuracy noticeably better than 1/10 then we can put that (1 — r)(fsQ = 1/10. In 
the gold we find for the energy e = e e = 2.61 TeV and the magnetic field H = 4 • 10 5 Oe 
that 1 — r = 0.371, the path of electron in the target is I ~ 1 cm and number of electron 
traversing the target without energy loss is iV ~ 3.2 ■ 10 -5 iVo. This estimates show that 
the measurement of the effect found in this paper will be feasible in the not very distant 
future. 



2.5 Propagation of high-energy photon in a medium in presence 
of an external field 

As known, the propagation of electromagnetic wave in a medium is defined by its dielectric 
tensor £ ik {oj). For relatively low frequency to (e.g. visible light) the dielectric tensor is 
defined by atomic phenomena. When the frequency of wave is much higher the atomic 
frequencies, the dielectric tensor has a form 

AtTTI P 2 

£ ik {u)=8 ik £{u), £H = 1--|, u 2 = e —, (2.75) 
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where luq is the plasma frequency, in any medium Uq < 100 eV. So that for to ^> luo 
an influence of atomic phenomena on propagation of electromagnetic wave in a medium 
becomes small. 

At very high energy the nonlinear effects of QED enter into game. One of them is the 
polarization of the vacuum by a photon. In the presence of an external electromagnetic 
field the polarization of vacuum was considered first in the pioneer papers [33]. In the 
strong field this effect can be essential for propagation of high-energy photons [36] , 

To evaluate the polarization tensor one has to consider the amplitude of photon scat- 
tering which included the polarization operator (221 • As above we use the quasiclassical 
operator method [7j, jH], [ID]- In this method the mentioned amplitude is described by 
diagram where the virtual electron-positron pair is first created by the initial photon 
with 4-momentum k(u, k) and polarization e! and then annihilate into final photon with 
4- momentum k and polarization e 2 . This corresponds to use of the non-covariant pertur- 
bation theory where at high energies (uj ^> m) the contribution of this diagram survives 
only. For this energy of photon this process occurs in a rather long time (or at a rather 
long distance) known as the lifetime of the virtual state 

h = (2-76) 

where q c > m is the characteristic transverse momentum of the process, the system 
h = c = 1 is used. When the virtual electron (or positron) is moving in a medium it 
scatters by atoms and changes the velocity under influence of external electromagnetic 
field. The mean square of momentum transfer to the electron from a medium and an 
external field on the distance // is 

Qf = (l 2 s+ Oik Qa = 4:7rZ 2 a 2 n a Ll fl L = L(q 2 c ) = ln(#?a 2 ), 

q F = eFl f , F = Ei + vxH (2.77) 

where a = e 2 = 1/137, Z is the charge of nucleus, n a is the number density of atoms in 
the medium, a is the screening radius of atom, Ej_ is the electric field strength transverse 
to the velocity of particle v ~ n = k/a>, H is the magnetic field strength (see Eq. (jl.21j) ). 

In the case of small momentum transfer qj = y^gj m the influence of a medium and 
an external field is weak, in this case q c = m. At high energy it is possible that q c > m. In 
this case the characteristic value of the momentum transfer (giving the main contribution 
into the spectral probability) is defined by the value of qf. The self-consistency condition 
is 



UJ 



9 9 2nu}Z 2 a 2 n a L (qA m 6 K 2 9 

q 2 c = qj = 2 — + -TT- >™ 2 , k = — eF 

3 qj Aqj m 6 

F 7 = E — n(nE) + n x H, (2.78) 

here k = \k\ is known parameter characterizing the pair production process in a homo- 
geneous external field F 7 . With q c increase the lifetime of the virtual state Eq. (j2.76|) 

decreases. 
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We will use the following normalization condition for the amplitude under considera- 
tion 

M = 2uAu. (2.79) 

The amplitude M is the contraction of the tensor ej ej,^' (e^ and e^) are the polar- 
ization vectors of the initial and the final photons) and the polarization tensor Mjk- We 
select the basic vectors as 

F 

- - e 2 = n x ei. (2.80) 



F 



Since the tensor Mj k is invariant under the space inversion then in the selected basic 
vectors it has the diagonal form 

M jk = \ M M H + M 22) + (eyeifc - e 2j e 2k )(M u - M 22 )} (2.81) 

In absence of external field it is convenient to describe the process of photon scattering 
using the helicity polarization vector e x (A = ±1) connected with momentum transfer A 
(see Eq. (j2.97|) in the next subsection). In presence of external field and for A = we 
choose the polarization vectors in the following way: 

e x = i (ei + i\e 2 ) , (e x e* x ) = 1, (e A e* A ) = 0, e A x n = i\e x . (2.82) 

In terms of helicity amplitudes M ++ and M + _ the tensor Mjk and the corresponding 
dielectric tensor Sjk has a form 

Mjk = 5 jk k 2 ++ + (e^eife - e 2j e 2k )k 2 + _, 

S jk = 6 jk -\M jk , k 2 ++ = M ++ , k 2 + _=M + _ (2.83) 

00 

The polarization tensor k 2 k is diagonal in the basic vectors ei and e 2 Eq. (j2.8()j) . The 
corresponding mass squared are 

k 1 = kn = k ++ + k^ — , k 2 = k 22 = k ++ — k^ — (2.84) 

The probability of pair creation by a photon with polarization e is 



(e) = — Im 



ee x fkl + {ee 2 fkf\ = --Im \k 2 ++ + &k 2 + _] , (2.85) 



where £3 is the Stokes' parameter. For unpolarized photon one has 

= — Im k 2 (2.86) 
v U) 

The method of calculation of the functions k\ + and k 2 + _ which define the tensor Mjk 
is similar to used in previous subsection (see details in [H2])- The general expressions for 
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photon masses squared under simultaneous influence of multiple scattering in a medium 
and an external electromagnetic field has a form 
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For k\_ we found respectively 
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where k is defined in Eq. fl2.78|) . Q is defined in Eq. fl2.3j) . the parameter g c is defined by 
the set of equations (compare with Eq. (j2.61|) ): 



Qc =\ for A(k' 2 + QL x ) < 1; 

4^ (V 2 £ 2 + QL C ) = 1 for 4(k' 2 + QLi) > 1 

In the absence of external field («/ = 0, ip p = 1) we have 
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where p = i/2u, ip(p) as above is the logarithmic derivative of the gamma function. 
Subsisting the result obtained into formula (|2.12|) we have the probability of pair creation 
which agrees with Eq. (j2.43|) . 

In the absence of multiple scattering [u — > 0) we get 
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where 

(f F = exp ( J (2-93) 

For this case the expressions for k\ 2 ()2.10|) after substitution results of Eq. ()2.92|) agree 
with masses squared of photon calculated in an external electromagnetic field (see jH] and 
references therein). 

In the case when the both characteristic parameters are small {y\ = AQL\ C 1, k < 
1), the main terms of decomposition of the functions k 2 are the sum of independent 
corrections to the photon mass squared both on account of the multiple scattering and 
an external field 
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UJ <ti UJ e , UJF = flT^T) UJ <C cu^, (2.94) 

F 

where a; e is defined in Eq. fl2.48J) and H in Eq. (jl.22j) . The correction oal/L\ follows from 
the first term of decomposition of the type Eq. (j2.20|) . Let us remind that in gold the value 
uj e is uj e = 10.5 TeV, this is the typical value for the heavy elements. 

In Fig. 7 the functions Re k\ + (curve 2) and Im k\ + (curve 1) are given for the case 
when the influence of a medium is taken into account only (Eq. (|2.91jl ). The both curves 
are normalized to the asymptotics given by Eq. ()2.94|) in the limit uj f — > oo. 

It should be noted that beginning with some photon energy uj = uj^ the radiative 
correction to the value Re Sjk in the absence of a field (k = 0) becomes larger than uJq/uj 2 
(see Eq. (|2.75jl ). The estimate of uJb is following 

uJb ~ xf^^r^uo- (2.95) 

V 71 A r 



For gold one obtains UJb ~ 40 GeV. 

A propagation of high-energy photons in oriented single crystals is one of interesting 
applications of the result obtained in this subsection. In this case we have both the dense 
matter with strong effect of multiple scattering and high fields of crystal axes or planes. 
As we saw, the LPM effect is most pronounced in the heavy elements. The same is valid 
for the process under consideration. Let the high-energy photon incident on crystal. The 
angle of incidence is small and such that the distance from axis g (or the distance from 
plane x) can be considered as a constant on the formation length of process (the constant 
field approximation is applicable, see Sections 12,15 in [TO]). 
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For orientation of a crystal along an axis the ratio of density of atoms in the vicinity 
of axis n(g) to the mean density n a is 



t ( n \ n ^ exp(-g 2 /2^) 

U{B) = — = 2,u\dn a ' (2 - 96) 

where u\ is the amplitude of thermal vibrations of atoms, d is the mean distance between 
atoms which form the axis. This ratio is maximal at g = 0. For numerical estimates we use 
for defmiteness the tungsten single crystal. For the axis < 111 > in W the ratio £ax(0)=370 
at the room temperature (T = 293 K) and £ ax (0)=1020 at T = 77 K. The effect of 
multiple scattering becomes strong at characteristic photon energy u) e {n a ) ~ 11 TeV and 
this value is inversely proportional to the density. So we have that ui e {g = 0) ~ 30 GeV 
at T = 293 K and uo e (g = 0) ~ 11 GeV at T = 77 K. It should be noted that within 
logarithmic accuracy we neglect by relatively small variation of L\ due to substitution the 
screening radius a 2 by the value 2u\. 

It is useful to compare these estimates with known threshold energies uj t at which the 
probability of pair creation in the field of axis is equal to the probability of the Bethe- 
Heitler mechanism, see Table 12.1 in ^H]- For photon energy u > u t the process of pair 
creation in the field of axis dominates. In W crystal for < 111 > axis uj t = 22 GeV at 
T=293 K and u t = 13 GeV at T=77 K. It is seen that for these energies the ratio u>/u) e 
which characterize the strength the LPM effect is of the order of unity. At u ~ Ut the 
maximal value of the parameter k(q) which determines the probability of pair creation by 
a photon in a field is also of the order of unity (at g ~ u±, see the mentioned Table). 
So we reach the conclusion that at some energy (for axial orientation of crystal) all 
the discussed effects are essential simultaneously. The analysis in this situation will be 
published elsewhere. For example, to calculate the influence of the field of axis on the 
polarization tensor one have to average the general formula (J2.87J) over all values of g 
(this is integration over d 2 g with the weight n±, where n± = n a d is the density of axis in 
the perpendicular to them plane). 



2.6 An influence of multiple scattering on coherent scattering 
of photon 

Investigation of the coherent photon scattering (the Delbruck scattering) has rather long 
history which can be found in review (see e.g. [HI]). There is a special interest to study the 
process for heavy elements (the cross section oc Z 4 ). However, in this case the contribu- 
tions of higher orders of Za into the amplitude of photon scattering are very important. 
This means that one needs the theory which is exact with respect to the parameter Za. 
The amplitudes of the coherent photon scattering valid for any Za for high photon energy 
oj ^> m and small scattering angle (or small momentum transfer A) were calculated by 
several group of authors (see jHH], [SHI, jlO], |H] |12])- In the last paper the process of the 
coherent photon scattering was considered in frame of the quasiclassical operator method 
which appears to be very adequate for consideration of this problem. 
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In the quasiclassical approximation the amplitude M of the coherent photon scattering 
is described by diagram where the electron-positron pair is created by the initial photon 
with 4-momentum k\ (u, kj.) and then annihilate into the final photon with 4- momentum 
k 2 (u, k 2 ), so that the photon momentum transfer is A = |k 2 — ki|. For high energy 
photon lj ^> m this process occurs over a rather long distance (see Eq. (|1.12(l ). 

It is convenient to describe the process of photon scattering in terms of helicity am- 
plitudes. We choose the polarization vectors with helicity A 

eA = —7= (ei + i\e 2 ) , e 1 = v=— e 2 = n x i/, n = — , 

A = ±l, e A e* = l, e x e*_ x = 0, e A xn = iAe A , (2.97) 

There are two independent helicity amplitudes: 

M++ = M__, M+_ = M_+, (2.98) 

where the first subscript is the helicity of the initial photon and the second is the helicity 
of the final photon. When the initial photons are unpolarized the differential cross section 
of scattering summed over final photons polarization contains the combination 

2[|M ++ | 2 + |M+_| 2 ]. (2.99) 

The calculation of the amplitudes of photon scattering on the separate atom appears to 
be most simple in the case when the screening radius of atom a s is much smaller than 
the formation length of process (a s <C X c u/2m = l p>max Eq. ()l.lip ). The imaginary part 
of helicity amplitudes in this case can be written in the form [43J 

Im M\\i = — dx dy Hxx'fxx 1 , (2.100) 

m J J 
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H ++ = 1 — 2x(l — x) + 4x(l — x)y(l — y), = x(l — x)y 

f ++ = \n(ma s ) - ~^^= In (s + VTT^) - f(Za) + 

2sy 1 + s z v ' 4z 

U- = l-^==ln(s + VT+^), s=^, (2.101) 
syl + s 2 v ' 2 

where the function /(£) is defined in Eq. (|2.5jl . The important property of Eq. (|2.100j) is 
that the dependence on the screening radius a originates in it from the Born approxima- 
tion. In this approximation in the case of arbitrary screening the radius a s enters only in 
the combination 
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\ + ql gn= n 9m x , q m = ^-- (2.102) 
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Because of this we can extend Eq. fl2.101|) on the case of arbitrary screening making the 
substitution 

- y/4 + *7 2 = lef, s=^- (2.103) 

U s ^Hef 

In the case a s ^> uo/m? (the screening radius is very large, or in other words we consider 
the photon scattering in the Coulomb field) we have to substitute in Eq. (|2. 101(1 

1 A Auj . . , 

a s -> — , s -> s c = — = —^x{l - x)y. (2.104) 
?|| 2g|| m 2 

In the case of the complete screening (a s <C uj/m 2 ) the functions /aa' are independent 
of x and y and the corresponding integrals are 

ii ii 

J dx J dy fi ++ = -, J dx J dy = — . (2.105) 



For the scattering amplitudes we have 



_ 28Z 2 a 3 c^ 2Z 2 a 3 a;^ 

Im M ++ = /++, Im M+_ = 2 /+_, 

Re M++ = 0, Re M + _ = 0. (2.106) 



The real part of amplitudes are calculated using dispersion relations. The photon scatter- 
ing amplitude in this case for arbitrary value of parameter a s m 2 /2u> for A = was found 
recently in |H] and for arbitrary A in [42J. 

When a photon is propagating in a medium it also dissociates with probability oc a 
into an electron-positron pair. The virtual electron and positron interact with a medium 
and can scatter on atoms. In this scattering the electron and positron interaction with the 
Coulomb field in the course of the coherent scattering of photon is involved also. There is 
a direct analogue with the LPM effect: the influence of the multiple scattering on process 
of the bremsstrahlung and pair creation by a photon in a medium at high energy. However 
there is the difference: in the LPM effect the particles of electron-positron pair created by 
a photon are on the mass shell while in the process of the coherent scattering of photon 
this particles are off the mass shell, but in the high energy region (this is the only region 
where the influence of the multiple scattering is pronounced) the shift from the mass 
shell is relatively small. To include this scattering into consideration the amplitude of the 
coherent scattering of photon should be averaged over all possible trajectories of electron 
and positron. This operation can be performed with the aid of the distribution function 
averaged over the atomic positions of scatterer in the medium in the same manner as it 
is done in Appendix A [43J. 

The photon scattering amplitudes can be written in the form 

M++ = M c ++ + MS, M+_ = M%_ + M [ +l, (2.107) 
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where Ml r+ is the main (logarithmic) term and M+\_ is the first correction to the scattering 
amplitude. The main terms are |43| 
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(2.108) 



where a s2 is defined in Eq. (l2.12|) . The structure of the amplitude Ml_ + is quite similar 
to the structure of the pair creation probability Eq. (j2.43|) under influence of the multiple 
scattering (the fact that the amplitude is expressed in terms of the probability is the com- 
mon feature of the coherent photon scattering description). This result is not occasional 
and connected with origin of coherent scattering of photon as a process which is going 
through dissociation of a photon into electron-positron pair. However the expression for 
the basic parameter v s differs from corresponding expression in Eq. (j2.43|) . 
The first corrections to the amplitudes are defined by 
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here the functions f x (z), f 2 (z) and g(z) are defined in Eq. (j2.45| ). £(s, a) is the generalized 
Riemann zeta function. Use of the given representations of functions F 3 (i/ S ) and G(z) in 
fi{z) simplifies the numerical calculation. The functions F\{z) and F 2 (z) in Eq. ()2.108|) 
encountered in the radiation theory. 

In the region of the weak effect of scattering (\v 8 \ <C 1, g c = 1) we obtain 
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here 



1 = 2 
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(2.111) 



a; e is denned in Eq. (j2.48|) and Li is denned in Eq. ()2.12|) . The characteristic energy uj e 
encountered in analysis of influence of the multiple scattering on the probability of pair 
photoproduction [2T], in gold a; e =10.5 TeV. The amplitudes Eq. (j2.110j) coincide with the 
formulas in absence of multiple scattering if we neglect the terms oc uj/uo e and (u/uj e ) 2 . 
The terms oc w/w e define the appearing real part of the scattering amplitudes while the 
terms oc (u/uj e ) 2 are the corrections to the imaginary part. 

In the region where the effect of scattering is strong (\u a \ ^> 1) we obtain 
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(2.112) 



the real and imaginary parts of the amplitudes are equal (if we neglect the term oc 1/L S 3 
in M++). Moreover the amplitudes Eq. (j2.112j) don't depend on the electron mass m. In 
place of it we have the value A a . 

The influence of a medium on the process of the coherent photon scattering illustrated 
in Fig. 8 a and b, where Im M ++ and Re M ++ as well as Im M + _ and Re M + _ are given 
as a function of photon energy u in gold. This influence is due to the multiple scattering 
of electron and positron of the virtual pair on the formation length of the process (see 
Eqs.flTHJ) and l03|) ) 
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(2.113) 



In the region A 2 < + m 2 this formation length is independent of A and it value is 
coincide practically with the formation length of pair creation by a photon l c considered in 
previous subsection. There is some difference connected with the logarithmic dependence 
of v 2 value on A 2 : 
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where q. 
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A 2 + a~ 2 , 



is the screening radius of atom ()1.13j) . This 



defines the weak (logarithmic) dependence of ImM ++ on A in the region A < yq 2 + rn 2 . 
This can be seen in Fig. 8a where the curves 1 and 3 represent behavior of Im M ++ for 
A=0.4435 m and A=0.0387 m respectively. For lower value of A the minimal momentum 
transfer q min (jl,13|) diminishes thereby the interval of contributing the multiple scattering 
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angles increases, so the multiple scattering affects the photon scattering amplitude at a 
lower energy (and smaller formation length). Because of this the curve 3 is shifted to the 
left respect the curve 1. Note that the curve 3 (A" 1 ~ a s2 ) is very similar to the curve 
2 in Fig. 4 which represents the behavior of the probability of pair photoproduction in 
gold vs photon energy. 

The new property of influence of a medium is the appearance of the real part of the 
coherent photon scattering amplitudes at high energy uj. In the region w<w e the value 
of ReM is small accordingly Eq. ()2.110|) . In the asymptotic region uj ^> uj e we have - 
Re M = Im M according to Eq. (|2.112|) . This property is seen clearly in Figs. 8 a and b. 
So the value of -Re M is small at low and very high energies of photon. At intermediate 
energies the value of -Re M have the maximum at uj ~220 TeV for A =0.4435 m and 
at uj ~80 TeV for A =0.0387 m. In Fig. 8b the same curves are shown for amplitude 
M + _. These curves are very similar to curves in Fig. 8a. The curves in both figures 
are normalized to imaginary part of the corresponding amplitude in the absence of the 
multiple scattering. The ratio these imaginary parts r is very smallat is r = 0.04435 for 
A=0.4435 m and r = 0.003018 for A=0.0387 m. 



2.7 The polarization of a medium and the bremsstrahlung pro- 
cess 

When one considers bremsstrahlung of enough soft photons uj < co>o7 (see Eq. ([1.4jl ). one 
has to take into account the effect of a polarization of the medium. This effect diminishes 
the formation length (see Eq. (jl.5j) ) as well as the probability of radiation (see ^j, the 
qualitative discussion may be found in [6 ). For analysis we use the general expression for 
the probability of radiation, see Eq. ljA.lj) . Appendix A. The factor in front of exponent 
in this expression (see Eq. (jA.2jO contains two terms A and B, the term A is not changed 
and the term B contains combination 

v ~tf + n— , «o = — , ujp = uj 'y, (2.115) 

UJ UJ 

and its dependence on k q (term of the order I/7 2 ) may be neglected also. So, the depen- 
dence on ujq manifests itself in the exponent of Eq. (jA.l|) . and respectively in the exponent 
of (lA~20l) only: 

uj 2 \ e' 
" 1 an = a, k = 1 + k 2 , Kq = — Kq. (2.116) 




2euJ ' u ' u e 

where the term with uj 2 . arises in the case of virtual photon: k 2 = k^k^ 7^ 0. Performing 
the substitution a — > a in Eq. (j2.1|) we obtain for the potential (see Eqs. ()2.1j) ()2.3|) . f!2.4|) . 

V(q) V(~q) = Qq 2 U Lj=\ - 2C\ = V c (q) + v(q), q=\q\ = qVZ, 
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- - 2 / -\ J , - , A Q T I ~\ I ''' 



V c (q) = QQ , g = QL(Qc), Q = L(g c ) = In 



vm = -lt{2C + ln^. (2.117) 

The substitution Eq. ()2.116|) in the expression for the probability of radiation Eq. (j2.2j) 
gives 

Rx -> i? 2 -> R 2 k = R 2 . (2.118) 

The value of the parameter g c in (|2.117jl is determined by equation (compare with 
Eqs.flrmj), (l2~2oT) ) 

Ag c A QL{g c ) = 1, for 4QL(1) > 1. (2.119) 
In the opposite case g c = 1 and this is possible in two intervals of the photon energy uo: 

1. for k,q <C 1 when the multiple scattering and effects of the polarization of a medium 
are weak; 

2. for k 3> 1 when effects of the polarization of a medium become stronger then effects 
of the multiple scattering (y Q < k). 

In an intermediate region we substitute g 2 — > g 2 c K in Eq. ([2.119j) . After it we obtain the 
equation for g c which coincides with Eq. (J2.10J) . see also Eq. (j2.25J) : 

\ = u 2 (g c ), u 2 (g c ) = 4QL(g c ). (2.120) 

Be 

Thus, for g c < 1 we have 



z7 = V4QL(<? C ) = — = — = -, L(g c ) = L(g c ), (2.12i; 



while for z/q < 1 we have 



,2 



v / lWD = 1 \/(?ln|^). (2.122) 



c 



The spectral distribution of the probability of radiation Eq. (j2.29|) with allowance for the 
polarization of a medium have the form 



dW a T 
-lm 



2L(g c ) 



(2.123) 



duo 2iv~f 2 
where 

i? 2 ) = $(i?i, 4), R 2 ) = F(R U R 2 ), 

In the case vq -C Kq> the characteristic momentum transfer in the used units (£ c ) are 
defined by value K^{g 2 = 1), one can use asymptotic expansion Eq. ([2.28) 1 and we have for 
the spectral distribution of the probability of radiation 

% + — -/(£«) =- L p + — -/(Za) , (2.124) 



dio 3 m?ujKQ \ p 12 / 37T m7 2 V p 12 
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where f(Za) is defined in f!2.5|) . L p = In (l83Z _1 / 3 Ko) . The results obtained agree with 
given in |4] where calculations are fulfilled within a logarithmic accuracy and without 
Coulomb corrections. It is seen that the dependence of spectral distribution on photon 
energy (udu) differs essentially from the Bethe-Maximon one (dw/u), and the probability 
is independent on the density n. 



3 Impact of a medium for a target of finite thickness 
3.1 General consideration 

In the case under consideration probability of radiation is defined not only by the relative 
time t = t 2 — t\ as in Sec. 2.1. The used radiation theory is formulated in terms of 
two times (see Eqs. f)A.l|) - (jA.3|) of Appendix A). Proceeding from this formulation we 
can obtain more general expression for the radiation probability. With allowance for 
the polarization of a medium we have for the spectral distribution of the probability of 
radiation (compare with Eq. (J2.2)) ) 

d 4a °° * 2 / *2 \ 

— = — Re / dt 2 I dtiexp —i / ji(t)dt 
dui to J J \ J 

— OO — OO \ t\ J 

x [n<p o {0, t 2 , h) - ir 2 V<p(0, t 2 , h)}, (3.1) 

where 

e' to 2 e' 2 

fj,(t) = 1 + KQg(t), k'q = -Kq, r x = —, r 2 = l + —, (3.2) 

here the function g(t) describes the density of atoms on the projectile trajectory, n is 
defined in Eq. fl2.115|) . The functions ip^(g,t 2 ,ti), (tp^ = (ip ,<p)) satisfy the equation 
(12. 3|) . but now the potential V depends on time 

i^ = H(t)<p„ H(t) = p 2 -iV(Q)g(t); 

<Po(Q,ti,ti) = 5(g), <p(Q,t 1 ,t 1 ) = p6(Q). (3.3) 

So, the functions and V(g) in Eqs. (|3.2|) and (|3.3p depend now on mean density of 
atoms on the projectile trajectory while the function g(t) describes the modulation of the 
density as a function of time 

m = ^ 

n 

Using an operator form of a solution of Eq. (j3.3|) (compare with Eq. (|2.6jl ) we can present 
the probability Eq. (|3.1j) in the form 

OO *2 / *2 \ 

-Re J dt 2 J dti exp I — i J ji{t)dt 



dw 4a, 
duj uj 



x (0\ ri S(h, h) + r 2P S(t 2 , tx)p\0) , S{t 2 , h) = T exp 



*2 

-% J H(t)dt 



ti 



(3.4) 
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where the symbol T means the chronological product. Note, that in Eqs. (j3.1|) and (|3.4|) 
it is implied that subtraction is made at V = 0, /x(£) = 1. 



3.2 Boundary effects for a thick target 

For the homogeneous target of finite thickness / the radiation process in a medium depends 
on interrelation between / and formation length If Eq. (jl.5j) . In the case when I > lj we 
have the thick target where radiation on the boundary should be incorporated. In the 
case when I < we have the thin target where the mechanism of radiation is changed 
essentially and in the case when I ~ If we have intermediate thickness. Integrals over 
time in Eq. ([3.4)1 we present as integrals over four domains: 

1. ti < 0, < t 2 < T; 

2. < t x < T, < t 2 < T; 

3. < h < T, t 2 > T; 

4. t x < 0,t 2 > T; 
where 

= tf(_t) + #(t - T) + Kd{t)d{T -t), T = -L = K = 1 + nl (3.5) 

i fo Lee 

In two more domains t± i2 < and t± j2 > T an electron is moving entirely free and there 
is no radiation. We consider in this subsection the case, when the thickness of a target / 
is much larger than formation length If Eq. (jl.5|) or (uq + k)T ^> 1. In this case domain 
4. doesn't contribute. The contributions of other domains are 

oo 

I\ ~ J dt\ J dt 2 exp (i(ti — Kt 2 )) exp (—iHt 2 ) exp (iH tx) - 



1 1 



H + kH 

h = J dt 2 J dtaexp (-i(H + K){t 2 - t x )) ~T J dr exp (-i(H + «)r) 



7 T 1 II 

- / rrfrexp(-2(i/ + fi;)r) = -z— — - + — — — -, 7 3 ~ - , (3.6) 

j H + k (H + kY Hn + 1 H + k, 



-oo 
T t 2 



T I , 1 I 

u 

where i^o — P 2 - The term in I 2 : —iT/(H + k) describes the probability of radiation 
considered in previous subsections. All other terms define the probability of radiation of 
boundary photons 1 . So, making mentioned subtraction we have for the spectral distri- 
bution of the probability of radiation of boundary photons 

dwh 4a 

— - = —Re (0 nM + r 2 pMp , 

M = {jTZ ^] ( 3 - 7 ) 

\H + k p 2 + IJ 



1 Radiation of boundary photons in an inhomogeneous electromagnetic field was considered in 
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In the case when both the LPM effect and effect of polarization of a medium are weak 
one can decompose combination in M in Eq. (j3.7|) 

1 1 1 , r , 9s 1 , , 

' V -^o)-^-T- (3-8) 



h + k P 2 + i p 2 + r u/ p 2 + i 

Using this decomposition one can find an estimate for probability of radiation 

dwh a 



dhJ 7T 



-c x v\ + c 2 Kq) , (3.9) 



where c\ and c 2 are some positive coefficients. So, in the case under consideration the 
probability of boundary photon radiation is negligible small. 

In the case when v\ < k 2 > 1 one can omit the potential V in Eq. (j3.8j) . so that the 
effect of polarization of a medium is essential, then 



dwh 4a 



duo c<j(27r) ; 



(Q\nM + r 2 pMp|0) dfp 



a f ( , 1 2 i > 
— ri H In K) +r 2 

TIUJ IV K K — 1 



K ~ 1 



In/? - 2 



(3.10) 



This result is the quantum generalization of the transition radiation probability [58J. 

Let us estimate the probability of radiation Eq. ()3.7j) in the case when both the LPM 
effect and effect of polarization of a medium are strong (i/q + k ^> 1). In this case we can 
neglect the term 1/(H + k) in the term oc T\ in M Eq. ()3.7|) . In the term oc r 2 one can 
put V(g) ~ UqQ 2 ~ ^o/P 2 so tnat 

P 2 (l^- ? Vl) 2d2 ^-^o + «). (3.11) 

Substituting these results into Eq. (|3.7jl we have up to logarithmic accuracy 

dwh a r , . . 

[ra+r 2 ln^ + «)]. (3.12) 



In the case ^Ck this result agrees with Eq. (j3.10|) . The case vq > k will be considered 
below (see Eas. fl3~3oT) and (CTTTIl ). 

3.3 A thin target 

We consider now a situation when the formation length of radiation is much larger than 
the thickness I of a target jTTj 

l«h= l f, C = l + 7 2 ^ 2 , T=-L« 1 -, (3.13) 
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where If, Ifo are defined in Eqs.(JT3J), ()1.6|) . In the case kT <C 1 the radiated photon 
is propagating in the vacuum and one can neglect the polarization of a medium. The 
spectral distribution of the probability of radiation from a thin target is 



dw 



th 



duo 



a 



7C 2 LU 



d 2 g nK 2 (g) + r 2 K 2 (g) (1 - exp(-V(g)T)) 



(3.14) 



where V(g) is defined in Eqs. ()2.3|) . ()2.4|) . K n is the modified Bessel function. If V(g = 
1)T <C 1 one can expand the exponent (the contribution of the region g ^> 1 is ex- 
ponentially damped because in this region Kq^(q) cx exp(— g)). In the first order over 
VT using the explicit expression for the potential Eq. (j2.4)l one obtains in this case the 
Bethe-Maximon formula with the Coulomb corrections Eq. ()2.29j) . 

In the opposite case when the multiple scattering of a particle traversing a target is 
strong (V(g = 1)T ^> 1, the mean square of multiple scattering angle ft 2 S> 1/7 2 )- We 
present the function V(g)T (see Eqs.JZSJ, (J23J) and pT4jl as 



k(? 1 



mf 
L t ln Qt, 



Ag 2 . ln = 
Qt 



Vln^= Aq 2 ln^-ln^ 
Q V Qt Qt / 



L t = ln 



Xt 
Q? 



In 



4a 



s2 



V c Q 2 t 



2C, 



(3.15) 



where Qt is the lower boundary of values contributing into the integral over g. Substituting 
this expression into ()3.14j) we have the integral 



(JXJ s 

2n J gdgKf(g) < 1 — exp 



-kg 2 | 1 ~ ln 

V U Qt, 



ttJ. 



(3.16) 



In this integral we expand the exponent in the integrand over 1/L t keeping the first term 
of the expansion. It is convenient to substitute z = kg 2 . We find 



exp(—z)zdz, 



1 °° 

J 2 = J K 2 (J^jexp(-z)lnzzdz, J = J x + J 2 . (3.17) 
* o 

Expanding the modified Bessel functions K n (x) at i C 1 and taking the integrals in the 
last expression we have 



J = Ji + J 2 
, nZ 2 a 2 



I + - ) (\nAk-C) 



— - I — 

2k + Lt 



-nJL t 



(3.18) 
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In the term with K$ in Eq. (|3.14|) the region g ~ 1 contributes. So we have 

OO OO 

J 3 = 2 J Kl(g) (1 - exp(-yr)) Qd Q ~2j K 2 {g)gdg = 1. (3.19) 



Substituting found J and J 3 into Eq. (j3.14|) we obtain for the spectral distribution of the 
probability of radiation in a thin target at conditions of the strong multiple scattering 

dwth a ( i i\ a on\ 

— — = — (n + r 2 J) . (3.20) 

du nuj 

The logarithmic term in this formula is well known in theory of the collinear photons 
radiation at scattering of a radiating particle on angle much larger than characteristic 
angles of radiation ~ I/7. It is described with logarithmic accuracy in a quasi- real 
electron approximation (see jSU], Appendix B2). 

The formula ()3.14|) presents the probability of radiation in the case when the formation 
length If I. It is known, see e.g. [Hj, that in this case a process of scattering of a 
particle is independent of a radiation process and a differential probability of radiation at 
scattering with the momentum transfer q can be presented in the form 

dW^ = dw s (q)dw r (q, k), (3.21) 

where dw s (q) is the differential probability of scattering with the momentum transfer 
q which depends on properties of a target. The function dw r (q, k) is the probability 
of radiation of a photon with a momentum k when an emitting electron acquires the 
momentum transfer q. This probability has a universal form which is independent of 
properties of a target. For an electron traversing an amorphous medium this fact is 
reflected in Eq. fl3.14|) . Indeed, passing on to a momentum space we have 

adu 



dw r {q, k) = I d 2 g \nK%(g) + r 2 Kl(g)\ (1 - exp(-iqg)) 
adu 



1TL) 



(3.22) 



where the functions F\(z) and ^(z) are defined in Eqs. (j2.109J) and ()2.108|) . Remind that 
q is measured in electron mass. The probability of radiation in this form was found in 
For a differential probability of scattering (here we consider the multiple scattering) 
there is a known formula (cp Eas. (|A~T3l . (IA~T8l) . Q and (JOJ)) 

dw s (q) = F s (q)d 2 q, F s (q) = J d 2 gexp (-iqg) exp (-V s (g)l) , 

V s (g) = n J d 2 q (1 - expHqg)) <r(q), (3.23) 
where <r(q) is the cross section of single scattering. 
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The formula ()3.23|) is the exact solution of the kinetic equation for the pure scattering 
problem (without radiation). Let us consider scattering of relativistic electrons in the 

case when scattering angles $ s <C 1 and the condition $ s 3> i?i = — is fulfilled (the 

a s e ___ 

angle $1 is defined in Appendix A, Eq. flA.13jl ). For combination V s (g)l in Eq. (j3.23j) one 
can use Eq. ()3.15|) . Let us remind that the impact parameters g measured in the Compton 
wavelengths A c , which is conjugate to space of the transverse momentum transfers q 
measured in the electron mass m. 



Vs (g)l^V(g)T = ^--^\n*- 2 , 



Qt L t g\ g 2 t 

1 

Lt = h + L t , — x = AL t 
Qt 

^2„2 i„2\ 



l t = \n(Axt) = \n(47rZ 2 a 2 n a la 2 s ) + 1-2C- 2f(Za). (3.24) 
Substituting Eq. ()3.24|) into Eq. ()3.23|) and passing to the variable v = qgt/2 we obtain 

dw s (v) = f(u)udu, 

/( u 2 u 2 u 2 \ 
uJ {uu)exp [ -— + — In — J du. (3.25) 



Solution of scattering problem in angular variables is given in the end of Appendix A (see 
Eq.flESU)). 

On adopted assumptions = A\t 3> 1, It ^> 1, Lt 3> 1) one can expands in 

Eq. (|3.25j) over powers of 1/L t : 



oo 1 «> 

f( u ) = S~i / M < / o(^)exp 

n=0 n 





.r n, (-) + ^/ (1) W + ^/ (2) W--- (•>.^» 
^ L t 




It is instructive to compare these results with the classical paper Bethe on multiple 
scattering jH]. The obtained formula Eq. (j3.26|) is consistent with Eq.(25) of [H], However 
there is difference in evaluation of the Coulomb corrections entering into the expression for 
the function l t Eq. (j3.24j) (see e b in Eq.(22) in 01]). In the last formula in 03] the Coulomb 
correction is contained in the form — ln(l + 3.3AZ 2 a 2 ) (it should be noted that a) Bethe 
was referring to the numerical coefficient at {Za) 2 calculated by Moliere 45j as "only 
approximate"; b) the Coulomb corrections to the cross sections of bremsstrahlung and 
pair creation by photon which contain the famous function f(Za) defined in Eq. (|2.5|) were 
found after publication of |44j ) while in Eq. ()3.24j) enters different correction: —2f(Za). To 
compare these results numerically let us consider case Za -C 1. In this case the correction 
in @3! is -3.3AZ 2 a 2 while for Eq.flHjS) we have the correction -2((3)Z 2 a 2 = -2A0AZ 2 a 2 . 
Besides, for accounting the inelastic scattering on electrons the use of Eq. (j2.34|) for the 
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potential V s (g) gives more accurate results than the simple substitution Z 2 — > Z(Z + 1) 
done in 

The main term of the expansion in Eq. ()3.26|) has the Gaussian form 



f i0 \v) = 2e~ v ' 



(3.27) 



The next terms contains the corrections to the Gaussian form. They were analyzed in 
detail in jUj . 

Using the formula (|3.22j) one can easily obtain to within logarithmic accuracy expres- 
sions (|3.2(Jj) . (j3.12|) . Both a radiation of boundary photons and a radiation in a thin target 
may be considered as a radiation of collinear photons (see e.g. [20] ) hi the case when an 
emitting particle deviates at large angle {$ s ^> I/7, q ^> 1). Using Eq. (j3.22|) at 2; > 1 we 
find 



dw r (q) 



adu r 



1TLU 



J d 2 qdw r (q)F s (c[ 



r i + r 2 (hi q 2 — 1 
adu r 



TILd 



+ r 2 (in q 2 



- 1 



(3.28) 



For a thin target value of q 2 is defined by mean square of multiple scattering angle on a 
thickness of a target I, and for boundary photons is the same but on the formation length 
If. However, if we one intends to perform computation beyond a logarithmic accuracy, 
the method given in this subsection has advantage since there is no necessity to calculate 
F s (q). Radiation from a thin target using the path integral method was considered in 
|E3j and [S3]. 

3.4 A target of intermediate thickness 

In the case when a target has intermediate thickness (I ~ If) the separation of contribu- 
tions on photon emission inside target and boundary photon emission becomes senseless. 
We consider this case ^B] neglecting by the correction term v(g). We present the spectral 
probability of radiation as 

00 T 

jM = J M = Jdhjd 


T 

j| m ) = J dt(T - t)e~ iftt 


00 00 
jj m ) = J dt x J dt 2 e~ 




+ r 2 Re 4 2) ~ 



{ti + tz)- m -N? 



■;{t L +t 2 +KT) 



( tl + t 2 + Ty m -N™ 
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v ' 



-it± 



/ it2 _ iKt2 \ f dt (T t) / it _ 
V / J t- iO) m V 



(t - z'O) 7 



oo oo 



+ dh dt 2 



-i(h+t 2 ) 







'(ti+t 2 + T) r 



(■ 



e lT - e~ lKT 



(3.29) 



where T is defined in (|3.2|) . the functions are 



iVi 



V 



sinh + vt\ cosh ' 

7/ 



No 



7/ 



sinh 7v(t — zO) 



'1 + v 2 t\t-2) sinh 7/T + v (t\ + t 2 ) cosh t/T ' 



1 + 2 
v = — —v , 



V2 



(3.30) 



here k and t/ are defined in Eq. (j2.116|) and in Eq. ()2.26|) . 

We confine ourselves to the case eCe e when the LPM effect manifests itself for soft 
emitted photons (uj <C e). So we can neglect the terms with r x in ()3.29|) . 

Let us consider first the case when LPM effect is weak iv\ <C 1). We assume here 
that condition vi{uj p ) > 1 (definition of u p see in Eq. (ll.5|) ) is fulfilled, that is in the region 
where v\ <C 1 one has u ^> u p and effects of the polarization of a medium are negligible. 
This is true for high energies (e > 10 GeV). Then for thickness T <ti \jv\ the transverse 
shift of the projectile due to the multiple scattering in a target as a whole have no influence 
on coherent effects defined by the phase = ul(l — nv) in the factor exp(— Indeed, 
for the projectile traversing a target in the case v{T C 1 an increment of the phase <fi is 
small 

A? 



v\T 2 < 1 



r 



(3.31; 



The angle of multiple scattering $ s is small also comparing with an characteristic angle of 
radiation I/7 (7 2 ^ = v\T <C 1). So, in the case v\ -C 1, v{F <C 1 the radiation originates 
on separate atoms of a target and an interference on target boundaries is defined by the 
value ioliX — v) = T. At T 1 this interference is weak, while at T 3> 1 there is a damping 
of the interference terms due to integration over photon emission angles. Expanding over 
ui in Eq. (j3.29|) we obtain (k = 1): 



j( 2 )(T)=Re £4 2) 0T) 
k=i 



{x - If 



1 + 3T 




•')/ / 7T^— sm{xT)dx 

3 T ' 

si(T) - 2Tci(T) + - sin T - - cos T 



For case T< 1 



J(2) (T) 



i_^ T + 6T 2 ( ln L + 1 

2 V T 



c 



and for case T 1 



J(2) (T) 



6- 



,cosT 



(3.32) 

(3.33) 
(3.34) 
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Thus, in the case ui -C 1, v{F <C 1 the probability of radiation is defined by Bethe- 
Maximon formula both for T < 1 and for T ^> 1. However, for T ~ 1 the interference 
on the target boundaries is essential. If we present as above as product of z/ 2 T '/3 
(Bethe-Maximon formula) and some interference factor, then this factor attains 0.53 at 
T = 0.32 (minimum of the interference factor) and 1.33 at T = 1.84 (maximum of the 
interference factor). 

When the parameter v{F is large [v\ <C 1, v{T ^> 1) the radiation is formed inside 
a target and the interference terms are damped exponentially. In this case formulae 
derived for thick target in Sec. 2.1 are applicable. Taking into account the contribution of 
boundary photons (see Sec. 3. 2) we have 

v\T ( 1 16i^ N 



J( 2 )( 



oo 



1 



2v\ 
~2p 



(3.35) 



3 V" ' 6Li 21 

where L\ and V\ are defined in Eqs. ()2.12|) and f)2.25j) . 

We consider now the case when the LPM effect is strong (z/ ^> 1) and the parameter 
T 1 while the value which characterize the thickness of a target z/ T ~ 1. We find for 
kT < 1 (r 2 ~ 2) 

*E _ 2a ( 7 (2) , j{2)\ T {2) U - D 2 

du 

J(2) 



J, 



In I 

T 



Re ^ ln(z/sinhz/T) - 1 - C 



7TT 



2i 



i/ tanh z/T 

vT 



VK 



ln(z/ tanh z/T) +1-C- — 



+ ikT In 



cosh z/T 



z/T 



tanh z/T tanh z/T 




sinh 2t sinh t 



v 



71 

( ' x l> l 1 ^ 



U . 



(3.36) 



For a relatively thick target (z/qT ^> 1) we have from Eq. (j3.36j) 




In 2 + 



V2 



7r 

K h In Z/n 

24 




(3.37) 



Here the terms without T are the contribution of boundary photons (see Eq.(4.14) in [T7j ) 
while the term cx T gives the probability of radiation inside target (with correction ~ k/vq 
but without corrections ~ l/£)- The relative value of the last corrections at z/ ^> 1 is 
given by Ea.(|2~33|). 

In the limiting case when a target is very thin and z/ T <C 1 but when z/qT > 1 we 
have from Eq. (j3.36|) 

2 A H^n + 



J(2) 



i 



C 



2 + 5, 



8 



(z/ T) 4 , 2(z/ T) 



+ ^T(lnz/ 2 T-C)-f(z/ T) 2 . 



180 



45 



(3.38) 
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The terms without 5 in this expression coincide with Eq. (j3.20J) (up to terms oc C/L t ). 

In the photon energy region where z/ T <C 1 the contribution of the terms (fc=l,2,3) 
is very small (~ 5) and decreases with photon energy reduction (oc to), so that in the 
spectral distribution of radiation only the terms J^, jjf"^ contribute. We obtain for the 
function in the case when (1 + vq)T <C 1 and when the parameter VqT, which char- 
acterizes the mean square angle of the multiple scattering in a target as a whole, has an 
arbitrary value 



Jfe- T 



7 U * J 7 

r rdti r dt 2 

J dxx JHJl[ exp 



1 + 



x 

oo 



{h+h)) l-exp 



4 J dxK\ (2y/xj [l - exp (-aci/gr)] = 2 J dggK 2 (g) [l - exp (-kg 2 ^ 



4k 



(3.39) 



where Ki(g) is the modified Bessel function. Formula ()3.39|) corresponds at k = 1 to 
result for a thin target obtained above (see Eq. (|3.14j l) without terms oc 1/L. Since the 
dependence on the parameter k is contained in Eq. (|3.39|) as a common phase multiplier 

(2) 

expy—inT), one can write more accurate expression for J\ (with terms oc 1/L) using the 
found results (see Eq. (|3.15jl ): 



J, 



(2) 



2e 



-inT 



dggKl{g) [1 - exp {-V(g)T)} . 



(3.40) 



For the case UqT ^> 1 it has the form 



JkT t(2) 
e j 4 



2k 
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(3.41) 



In the case when parameter k is not very high one has to use an exact expression Eq. (j3.14j) . 
For k <C 1 one can expand the exponent in the integrand of Eq. 1)3.40)1 . Then we find 



JkT t(2) 
c J 4 



1 + 



6Li 



AnZ 2 a 2 nl 



At kT <C 1 the spectral distribution of probability is 



dw 
doj 



2a 
3tyuj 



v{T 1 + 



6Li 



(3.42) 



(3.43) 



This is the Bethe-Maximon formula for not very hard photons (terms oc ^— \ are omit- 
ted). 
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When a photon energy decreases, the parameter k increases as well as the combination 
kT oc 1/uo, while the value {vqT) 2 decreases oc to. Just this value defines an accuracy of 
Eq. (j3.4U|) . Using Eq. ()3.29|) at T 1, kT > 1 we find for the probability of transition 

(2) 

radiation (J5 ) the following expression 



duo 



+ ^sin(KT) 




dwtr 2a I /. 2 

rv 



In k - ci(kT) + cos(kT) (In T + C) 
+ kTsi(kT) - 4 sin 2 — \ . (3.44) 



In the limiting case kT 3> 1 the probability Eq. ()3.44|) turns into standard probability of 
the transition radiation (j3.10|) with oscillating additions 



dw tr 2a 



J tr + cos(kT) (InT + C + 1) + I sin(KT) 



duO 7T 

J tr = (l + ^— In ft - 2. (3.45) 

Note, that there is a qualitative difference in a behaviors of interference terms in Eqs. ()3.34|) 
and (|3.45|) . In the former an amplitude of oscillation with uo increase decreases as 1/uo 2 
while in the latter the corresponding amplitude weakly (logarithmically) increases with uo 
decrease. 

From the above analysis follows that in the case when uqT <C 1 (v® ^> 1) the spectral 
distribution of probability of radiation with the polarization of a medium taken into 
account has the form 

aJ= a7 + cos(kT) ^7' (3 ' 46) 

where dw t h/duo is the spectral distribution of probability of radiation in a thin target 
without regard for the polarization of a medium. In the case 4k = v^T 3> 1 the probability 
dwth/duo is defined by Eqs. (|3.18|) . (|3.2U|) and for the case k 1 it is defined by Eq. (j3.35| ). 
More accurate representation of the probability of radiation dw t h/duo may be obtained 
using Eq. (j3.40J) . It follows from Eqs. (j3.45|) and ()3.46|) that if we make allowance for 
multiple scattering at kT ^> 1 this results in decreasing of oscillations of the transition 
radiation probability by magnitude of the bremsstrahlung probability in a thin target. 

The radiation of the boundary photons with regard for the multiple scattering was 
considered in jHI] (for uo e), the effect of polarization of a medium was included into 
consideration of the problem in [52^ and j^S] • I n these papers the probability of radiation of 
boundary photons (under condition of applicability of Eqs. ([3.37|) and ()3.38jl ) was analyzed 
also to within the logarithmic accuracy (see Eq.(20) in and Eq.(15) in jSHl)- This 
accuracy is insufficient for parameters connected with experiment ^1]-PE]- For example, 
for e = 25 GeV and in heavy elements the value h> equates k for z/ ~ 20. One can see 
from Eqs. (j3.37|) and ()3.38|0 that in this case Inuo is nearly completely compensated by 
constant terms. Our results, which are consistent with obtained {52], are more accurate 
and presented in more convenient for application form and the Coulomb corrections are 
included. 
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3.5 Multiphoton effects in energy loss spectra 



It should be noted that in the experiments [12]- [HI the summary energy of all photons 
radiated by a single electron is measured. This means that besides mentioned above 
effects there is an additional " calorimetric" effect due to the multiple photon radiation. 
This effect is especially important in relatively thick used targets. Since the energy loss 
spectrum of an electron is actually measured, which is not coincide in this case with the 
spectrum of photons radiated in a single interaction, one have to consider the distribution 
function of electrons over energy after passage of a target ^H]- As it is known, this 
distribution function is the solution of the corresponding kinetic equation. However, the 
problem can be simplified essentially if one is interested in the soft part of the energy loss 
spectrum. Just this situation was in the experiments smce the measurements 

were performed in the region of the photon energies of one to five orders of magnitude 
lower than the electron energy. 

We consider first the spectral distribution of the energy loss. The probability of the 
successive radiation of n soft photons with energies 

u>i, u>2, ■ ■ -0J n by a particle with energy e {oJk <C e, k = 1, 2 . . . n) on the length I in the 
energy intervals duj\duj2 ■ ■ ■ duj n is given by expression 

l ll ln-1 

dw(uJi,uj2,...itJn) = aJ dW{u\)dl\ J dW{u 2 )dl2 ■ ■ ■ J dW(u n )dl n 



= — dw(u>i)dw(uj2) ■ ■ ■ dw(uj n ), (3.47) 
n\ 

where A is the normalization constant, dW(uj)/duj is the differential probability of the 
photon radiation per unit length, dw(uj)/duj = IdW (uj) / duo is the differential probability 
of the photon radiation per length I. If the probability dw/duj doesn't depend on the 
particle energy e then integrating Eq. (|H.47jl over all photon energies we obtain 

Wn = — w n , w = / —du. (3.48) 
n\ J duj 

The value A is defined by the condition that probability of all the possible events with 
radiation of any number of photons or without photon radiation is equal to unit. 

OO OO 

V w n = — w n = Aexpw = I, A = exp(-w). (3.49) 

n=0 n=0 

Using Eqs. (|3.47j) . (|3.49|) we can write the expression for the differential distribution of the 
energy loss in the form 



1 de ^ 1 , . f dw dw dw . ( J\ \ , 

-— = > — exp(-w) / — — 6 > u> k - uj \dui . . . dw r , 

uj duj £r x n\ J dui duo 2 du n \^[ j 



(3.50) 



Here on the right-hand side we have the sum of the probabilities of radiation of n photons 
with summary energy uj. 
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Using standard parametrization of 5-function 
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we obtain 



de oj 
du) 2tt 



oo oo / ^ \ n 

J exp (isoj — w)^2 —lj — — exp(— isu>i)du>i j rfs 



i 

oo 



2^ 

1 

TT 
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Re/ exp (zx) exp -/^ 



1 — exp 



CO 



duj\ \ dx 



(3-51) 



The formula (|H.51|) was derived by Landau [001 (see also PQ) as solution of the kinetic 
equation under assumption that energy loss are much smaller than particle's energy (the 
paper [SU] was devoted to the ionization energy loss). Let us notice the following. The 
energy loss are defined by the hard part of the radiation spectrum. In the soft part of 
the energy loss spectrum Eq. (j3.5ip the probability of radiation of one hard photon only 
is taken into account accurately. To calculate the probability of the emission of two and 
more hard photons one has to take into account step by step the recoil in previous acts of 
the photon emission. The probability of radiation of two and more hard photons is of the 
order (l/L ra d) 2 and so on. Thus, the formula (j3.51J) is applicable for the thin targets and 
has accuracy l/L ra d. If we want to improve accuracy of Eq. (J3.51|) and for the case of thick 
targets I > L ra d one has to consider radiation of an arbitrary number of hard photons. 
This problem is solved in Appendix of ^21 for the case when hard part of the radiation 
spectrum is described by the Bethe-Maximon formula. In this case the Eq. (13.51)1 acquires 
the additional factor. As a result we extend this formula on the case thick targets. 

We will analyze first the interval of photon energies where the Bethe-Maximon formula 
is valid. We write it in the form (within the logarithmic accuracy) 



10- 



dw 
du 



I 



J rad 



(3.52) 



where I is the thickness of the target, L ra d is the radiation length, 

4/ 



3L 



rad 



UJ 

x x = -. 

e 



The integral in the curly brackets in Eq. 1)3.51)1 is 

J ~ ft (inx + /i + i~V u = ln— -\ + C 
\ 1) uj 8 



(3.53) 
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Substituting this result into Eq. (|3.51j ) we have 

de _ exp(-/fy) 

where T(z) is the Euler gamma function. If we consider radiation of the one soft photon, 
we have from Eq. (j3.51j) de/duj = (3. Thus, the formula (J3.54JI gives additional "reduction 
factor" fsH which characterizes the distortion of the soft Bethe-Maximon spectrum due 
to multiple photons radiation. The derivation of Eq. ()3.54j) is based on expression (j3.51|) 
in which the statistical independence of photon radiation acts is assumed (the Poisson 
distribution). This statistical independence is broken when two or more hard photons 
(energies of which are of the order of the energy of the initial electron) are emitted. 
Since the probability of hard photon radiation is of the order of /3, Eq. ([3.54)1 is applicable, 
strictly speaking, at (3 <C 1 and it has the accuracy up to factor (l + 0(j3 2 )). However, the 
contribution of arbitrary number of hard photons radiated in the Beth-Maximon region 
is exactly calculated. This contribution is contained in the factor g((3) (see Appendix of 
[T5], Eqs.(A.15), (A. 16)). So we obtain the reduction factor valid for (3 > 1 

fBB = = (-Y (i + (i + CV /4 • (3.55) 



r(l + /3) \e) 

In Fig. 9 the function /bh(u;) is given for electron with the energy e = 25 GeV. One can 
see that for uj = 100 MeV the difference between the value fsH^) and unit is 5-6 times 
greater than value (3. Let us discuss this circumstance. The emission of accompanying 
photons with energy much less or of the order of uj (we consider in this figure the situation 
connected with SLAC experiment ^2]-|II] where uj <C e, /3 < 1) changes the spectral 
distribution on quantity order of f3. However, if one photon with energy uj r > uj is emitted, 
at least, then photon with energy uj is not registered at all in the corresponding channel of 
the calorimeter. Since mean number of photons with energy larger than uj is determined 
by the expression (see Eqs. ([3.47|) - ([3.51|l ) 

00 f dw 

V nw n = w UJ = —du, (3.56) 
n=o I duJ 

the probability of the event when no photon with energy uj r > uj is radiated is defined 
by exp (—w u ). This is just the main factor in the expression for the reduction factor f BH 
(13.55)1 . In the case, when radiation is described by the Bethe-Maximon formula the value 
w w increases as a logarithm with uj decrease (w u — (3 In e/uj) and for large ratio e/uj the 
value Wu is much larger than f3. Thus, amplification of the effect is connected with a large 
interval of the integration (uj — e) at evaluation of the radiation probability. 

The Bethe-Maximon formula becomes inapplicable for the photon energies uj < uj c , (uj c 
is defined in Eq. (BTTTjl ). where LPM effect starts to manifest itself (see Eqs. (l2~T7)l . 
Calculating the integral in Eq. ()3.51j) we find for the distribution of the spectral energy 
loss 

SH^&™' /-M = ff OT(l + ^jDexp(- Wc ), (3.57) 
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where Jlpm is the reduction factor in the photon energy range where the LPM effect is 
essential, 



UK. 



C 2 ~ 1.959. 



(3.58) 



In this expression the terms ~ 1/L (see Sec. 2.1) are not taken into account. In Fig. 10 
the function /lpa/(^) is shown for u <100 MeV (e = 25 GeV, uj c = 228 MeV). It is seen in 
Fig. 10 that the reduction factor /lpm changes appreciably in the region of high photon 
energies solely. This is due to the fact that here the total probability of photon radiation 
is finite (in contract to the Bethe-Maximon formula) and the integral which defines this 
probability converges at uj — > 0. 

In the above analysis we neglected an influence of the polarization of a medium on 
the bremsstrahlung. This is correct if k <C vq (see Eqs. (j2.115|) . (J2.1 16)1 . (J2. 26)1 . Sec. 2.6) 
where 



^'7 



K 



UJ, 



1 + -| = 1 + Kq, uj p = 7^0, = — • (3.59) 

UJ UJ 

We assume, for definiteness, that uj v <C uj c . This is true in any case for a dense matter if 
electron energy e > 10 GeV . In the opposite case when k > ^ > 1 there is an additional 
suppression of the bremsstrahlung (see Sec. 2.6): 



dw 



31 



dui AnL r 



ad 



K 



(3.60) 



and this contribution into reduction factor (J3.57|) can be neglected. 

The main contribution into effect considered for the photon energies such that uj uj p 
gives the transition radiation Eqs. ()3.10|) . ([3.45)1 . Let us note, that for the transition radi- 
ation the given above derivation of the formula (|3.51j) is directly inapplicable. However, 
for application of formula (|3.51j) it is enough that acts of radiation of soft photons are sta- 
tistically independent. Taking integrals in Eq. 1)3.51)1 and taking into account that rj 1 
(but it may be that rja ~ 1) we obtain 



de g(j3) 



du 



71 



exp(— w c ) exp 



-7] \a - 



3vr 2 



x sin [2TT7] (a + - 2rja))] T(l - 2rja), 



where if)(z) = dlnT(z)/dz, g(/3) is defined in Eq. ([3.55)1 . 

r) = — , <7 = ]n-£ + (7-l. 

7T UJ 



(3.61) 



(3.62) 



In the case i]o 1 we have 

2r] ( ln ~uj ~ / 9 ^ ex P( _Wc )^' 
f tr = exp 



de 
du 




(3.63) 
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The case of a thin target, where the photon formation length 

l f = — (K + u )>l, (3.64) 

I is the target thickness, should be analyzed separately. We consider situation when effects 
of the polarization of a medium are weak (kq < u ) In this case the scattering takes place 
on a target as a whole during the radiation process and the spectral probability dw/duj 
in Eqs. (l3.48|) - (j3.51|) depends on the momentum transfer q = 7$ (see Eq. fl3.22j) ) 



dw r 2adu_ F Jg\ 



On the final step one have to average de r /du over q with the distribution function F s (q) 
Ea. (l3~23l . 

Since the expression for the spectral distribution ()3.65j) has the same form as in the 
Bethe-Maximon case, one can use Eq. 1)3.54)1 for calculation of the reduction factor. As a 
result we have for the contribution of the thin target region 



fth = jim {pliq)exv 



x 2a (g (I) 



7T (F 2 U)) 



ln^ ; (3.66) 



where u;^ is the boundary energy starting from which the target becomes thin. The 
lower limit of the applicability of Eq. ()3.66j) is determined from the condition If — L The 
expression (J3.66)) depends essentially on the mean square of the momentum transfer q 2 
associated with the target thickness. In the case q 2 <C 1 when an influence of the multiple 
scattering is weak {ft 2 <C I/7 2 ) 



AZ 2 a 2 nl A c 
(q 2 + a 2 ) m 2 a s2 



F s(l) = , „ , on2 -^ « s = — , (3-67) 



where a S 2 is the effective screening radius Eq. (j2.14|) . When the multiple scattering is 
strong (q 2 S> 1, ft 2 ^> I/7 2 ), the reduction factor can be calculated using the Gaussian 
distribution for F s (q) 



1 / q 2 \ 2 47rZ 2 a 2 , 7 g 2 rfg 2 4vrZ 2 « 2 „ g 2 + 1 
— exp — - , q s = — nl / — - — ~ — nZ In — . 



In the first case (q 2 <C 1) one has (F 2 2 ) <C (F2) and one can neglect the correction to unit 
in Eq. (j3.66J) . In the opposite case (q 2 ^> 1) the main contribution into entering mean 
values of F 2 gives values g > 1 where 

F 2 (q) ~lng 2 -l + ^. (3.69) 

In this case the effect under consideration may be noticeable. 
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From the above analysis one can find an approximate expression for the reduction 
factor in the general form. One can present the integral over u)\ in Eq. (j3.51|) as 



f dw \ ( . uj\\ 1 , f dw , f aw 

/ - — exp —ix — — 1 auJi ~ — / - — duoi + / - — 
J duj\ V \ uj J J J aui J auJi 



exp 



-ix — ) — 1 

UJ 



duj!. (3.70) 



The second integral on the right-hand side of Eq. (J3.7U|) (as well as terms omitted in the 
first integral) has an order udw/duj. For the case (3 <C 1 this value is small 



dw 41 
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ln^-1) =2r)(a-C), rj 
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71 



(uj uj p ) . 



(3-71) 



Let us note that a contribution of the transition radiation can be enlarged n times if one 
makes a target as a collection of n plates conserving the total thickness / provided that 
definite conditions are fulfilled for the plate thicknesses and gaps between plates. In that 
case 77 — > nr\ and above formulae are valid if nr\ < 1. 

Expanding the exponential in Eq. (j3.51|) with the second integral of Eq. (j3.70|) and 
integrating over x we have 
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27 



exp(ix)dx 
dw 



duj\ 
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UJi 



dw 
dui 



exp 



-ix- 



UJi 



UJ 



duj\ 



UJ 



duj\ = uj 



dw 
duj 



(3.72) 



Thus, the spectral distribution of the energy loss and reduction factor / have the following 
general form 



de dw 

T~ = J = ex P 

duj duj 



dw 
dui 



duj\ 



(3.73) 



3.6 Radiation from structured target 

The radiation from several plates for the relatively hard part of the spectrum in which 
the bremsstrahlung in the condition of the strong LPM effect dominates was investigated 
recently in jHI]. A rather curious interference pattern in the spectrum of the radiation 
was found which depends on a number (and a thickness) of plates and the distance 
between plates (the polarization of a medium was neglected). In the present subsection 
the probability of radiation in a radiator consisting of iV plates is considered |20j . The 
transition radiation dominates in the soft part of the considered spectrum, while the 
bremsstrahlung under influence of the strong LPM effect dominates in the hard part. 
The intermediate region of the photon energies where contributions of the both mentioned 
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mechanisms are of the same order is of evident interest. We analyze this region in detail. 
In this region effects of the polarization of the medium are essential. 

Here we consider the case when the target consists of N identical plates of thickness 
h with the equal gaps l 2 between them. The case h <C If will be analyzed where If 
is the characteristic formation length of radiation in absence of a matter (uo = 0), see 
Eqs.flUD-dH) 



/ 



fo 



h 



1 + 7 2 $2 1 + p 



i 



fa 



k 
h 



;i + p^Tt < 1, 7\ + T 2 = T, (3.74) 



where d c is the characteristic angle of radiation. 

One can split the spectral distribution of the probability of radiation into two parts 



dw dwbr dw tr 
duo duj duj 



(3.75) 



where dwbr/duo is the spectral distribution of the probability of bremsstrahlung with al- 
lowance for the multiple scattering and polarization of a medium and dwt r /duj is the 
probability of the transition radiation for the iV-plate target. For dw^ /duo one has: 



dw 



(AT) 
br 



duo 



a 

71UO 



N-l 

Re ^ / dr 2 / dr\ exp [— i (r x + r 2 + kXi + (n — 1)kT)\ 



n\ ,ri2=0 



X 



n {G n -G n {0))+r 2 (G 2 n -G 2 n (0) 



(3.76) 



where 



G n 1 — Pn (d n — d. 



■n-i) - 1 +i(d n - d n -i) (n + r 2 ) - (3 n 1 t 1 t 2 ; 

N > n = n 2 — n\ + 1 > 1; rii = 0, < T\ < oo; n 2 = N — 1, < r 2 < oo; 

n x >l, < n < T; n 2 < N - 2, < r 2 < T, 
2 li 



1 + K , K — K\ 



h + h 



(3.77) 



Note for the subtraction procedure one has that when q — > 0, the function b, (3 n — > oo. 
The function d n can be presented in the form 



sinh nrj / ; r 

d n = — r-; , sinhr] = 2JiqT 1 T (1 + iqTiT). 



sinh?] 

For the case r\ 1 one has 

t? ~ 2 iqTiT = 2T^iq 

In above formulae 

T = T x + T 2 , 6 



h 



h + h 



' 0n=^, G n (0) = G n (g = 0) 



(3.78) 



(3.79) 



(3.80) 
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For one plate (N = 1) we have ri\ = n 2 = 0, n = 1 and 



cr 1 



i(ri + t 2 ) - -nr 2 , Gi ^0) = i(n + r 2 ) 



(3-81) 



In the integral ()3.76|) we rotate the integration contours over n, r 2 on the angle —tt/2 and 
substitute variables r lj2 — > —ix 1)2 . Then we carry out change of variables x = Xi+x 2 , x 2 = 
zx. We have 



dw 
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1 r 



du 
ab 
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TXU) 



cos kT\ J exp(— x)dx J 
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n l- 



+ r 2 1 



g 2 (x,z) 



f/2 



= — cos(kTi) / exp(— bx) 
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T 

g(x,z) 



riF x 



I )+r 2 F 2 (,/ T 



l + -z(l-z), 



(3.82) 



where the functions Fi(x) and -F 2 (a;) are defined in (J2.1U9)) and (J2.1U8)) . 

The special case is when a target consists of two plates (N = 2). Since the formation 
length is enough long for the soft photons (u <C e) only, we consider the term with r 2 in 
Eq. ()3.76p as far as r\ = u 2 /e 2 <C 1. For the case (N = 2) the sum in Eq. ()3.76|) consists 
of three terms: l)n\ = n 2 = 0; 2)n 1 = n 2 = 1; 3)ni = 0,n 2 = 1. For the two first n = 1 
and we have from Eq. (|3.81|) 



dwfX dwfX ar 2 



du 



du 



Re 



1X0J 



exp(-mTi) y g?t 2 / exp(-z(ri + r 2 )) 



x 



1 



(ri + r 2 ) 2 (n + r 2 + zrir 2 /6) 2 
For the third term n = 2, d 2 = 2 cosh?] = 2(1 + 2iqTiT) (see Eq. (j3.78j) ) and we have 



(3.83) 
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(3.84) 



so that 



dWfo 3 _ ar 2 

dw TXUJ 



Re 



CXJ CXJ 

exp(— «(kT + kTi)) J dx\ J exp(— (x\ + x 2 )) 



{iT + xi + x 2 ) 2 \^ iT + Xl + X2 + 2 XiX2 + ix [ Xl + X2 + 



dx 2 



(3.85) 



Here we rotate the integration contours over n, r 2 on the angle —tt/2 and substitute 
variables T\ 2 — > — zxi 2 . 
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In the case of the weak multiple scattering (6 ^> 1) neglecting the effect of the polar- 
ization of a medium (k = 1) and expanding the integrand in Eqs. (j3.83|) and ()3.85|) over 
1/6 we have for the probability of radiation 



,(2) 



,(2) 
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Note that the main term of the decomposition in Eq. (J3.86|) is the Bethe-Maximon prob- 
ability of radiation from two plates which is independent of the distance between plates. 
This means that in the case considered we have independent radiation from each plate 
without interference in the main order over 1/6. The interference effects appear only in 
the next orders over 1/6. 

In the case of the strong multiple scattering we have 6 <C 1, p 2 c 3> 1. When the 
formation length is much longer than the target thickness as a whole (T 6) we find for 
the probability of radiation 
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In the opposite case & < T < 1 neglecting the effect of the polarization of a medium 
(kTi <C 1) we have 



dw 



(2) 
br 



du 



ar 2 



TYU 



1 i 26) jhip + l-C) -2 



(3- 



The case of the strong multiple scattering (6< 1) for the photon energies where the 
value T > 1 is of the special interest. In this case we can neglect the polarization of a 
medium k — K = 1, and disregard the terms oc kT\ in the exponent of the expressions 
(I3.83J) and ()3.85|) since Ti <C 1. We obtain as a result 
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where si(z) is the integral sine and ci(z) is the integral cosine. At T ^> 1 we have 
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Let us turn to the case of large N (N > 3). If the formation length of the bremsstrahlung 
is shorter than the distance between plates (T > 1) the interference of the radiation from 
neighboring plates takes place. Using the probability of radiation from two plates ()3.86|) 
we obtain in the case of weak multiple scattering (b ^> 1) 



dwl!' Nar 2 
duj 3nujb 



-4 + 4^> 



(3.91; 



where for T ^> 1 the function G(T) is defined in Eq. (|3.8fij) . In the case of strong multiple 
scattering (6< 1) and large T we have (compare with Eqs. (j3.89J) and ()3.90p ) 

The behavior of the spectral distribution uj—— can be discussed using as an example 

duj 

the case of two plates with the thickness l\ and the distance between plates l 2 > h which 
was analyzed in detail above. For plates with the thickness l\ > 0.2%L ra d and in the 
energy interval uj > uj p , in which the effects of the polarization of a medium can be 
discarded, the condition (|3.74|) is fulfilled only for enough high energy e so that uj p <C uj c 
(see Eqs. (ll.5|) . (|2.27j) ). We study the situation when the LPM suppression of the intensity 
of radiation takes place for relatively soft energies of photons: to < uj c <C e. 

For the hard photons uj c <C uj < e the formation length Ifo Eq. ()1.6|) is much shorter 
than the plate thickness l\ (Xi ^> 1), the radiation intensity is the incoherent sum of 
radiation from two plates and it is independent of the distance between plates. In this 
interval the Bethe-Maximon formula is valid. 

For uj < uj c the LPM effect turns on, but when uj = uj c the thickness of plate is still 
larger than the formation length Iq 

U m / m 27T /1 

—4 = T x {uj c ) = T c ~ — -J- > 1, (3.93) 

so that the formation of radiation takes place mainly inside each of plates. With uj 
decreasing we get over to the region where the formation length l c > h, but effects of the 
polarization of a medium are still weak (uj > uj p ). Within this interval (for uj < uj t h) the 
main condition 1)3.74)1 is fulfilled. To estimate the value ujth we have to take into account 
the characteristic radiation angles (p 2 c in Eq. (J3.74j) ). connected with mean square angle of 
the multiple scattering. Using the definition of the parameter b = l/(4gTx) in Eq. ()3.80J) ) 
we have 
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where the parameter gt, is defined by a set of equations (we rearranged terms in (|3.15j) ): 
g b = 1, L b = Li for AQL{Ti = — - — < 1, q = QLr, 

AQL{ Qh )T lQ l = —-^qI ( 1 - ^4 J = 1, for 4QL 1 T 1 > 1, (3.95) 
a L rad \ Li ) 

here L\ = ln(a S 2/A c ) 2 - One can show (see discussion after Eq.(3.6) in [IBj) that ui t h — 4a>6- 
Naturally, uj t h < ^c/T c , and when uj = uj c /T c one has T\ — 1, l\ — l . It is seen from 
Eq. (j3.94j) that when the value l\ decreases, the region of applicability of results of this 
subsection grows. 

So, when uj < uj t h the formation length is longer than the thickness of the plate 
li and the coherent effects depending on the distance between plates Z 2 turn on. For 
the description of these effects for T = {\ + l 2 /l\)Ti > 1 one can use Eq. (j3.89| ). For 
T > 7r ^> 1 one can use the asymptotic expansion (J3.9U|) and it is seen that at T = 7r 
the spectral curve has minimum. Let us note an accuracy of formulas is better when uj 
decreases, and the description is more accurate for T 3> T\ (/ 2 ^> h). 

With further decreasing of the photon energy u the value T diminishes and the spectral 
curve grows until T ~ 1. When T < 1 the spectral curve decreases oc InT according with 
Eq. (|3.88|) . So, the spectral curve has maximum for T ~ 1. The mentioned decreasing 
continues until the photon energy a> for which (1 + 2/b)T ~ 1. For smaller uj the thickness 
of the target is shorter than the formation length. 

The results of numerical calculations are given in Fig. 11. The formulas ()3.83|) . ()3.85|) 
were used. The spectral curves of energy loss were obtained for the case of two gold plates 
with the thickness l\ = 11.5 fim with different gaps I2 between plates. The initial energy 
of electrons is 25 GeV. The characteristic parameters for this case are: 

oj c ~ 240 MeV, T c ~ 2.9, r 1 ~ 3.3, u th ~ 80 MeV, 

w p ~3.9MeV, c^-SOkeV, — = k = ±-±-1 = 3, 5, 7, 9, 11. (3.96) 

At uj > 80 MeV the radiation process occurs independently from each plate according 
with theory of the LPM effect given above. The interference pattern appears at uj < 
80 MeV where the formation length is longer than the thickness of one plate and the 
radiation process depends on the distance between plates T. According to Eqs. (|3.90j) . 
(I3.91|) the curves 1-5 have minimums at uj ~ nuth/k (T = n) which are outside of Fig. 
11 and will be discussed below. In accord with the above analysis the spectral curves 
in Fig. 11 have the maximums at photon energies u ~ uJth/k (T — 1). These values 
(in MeV) are uj ~ 27,16,11,9,7 for curves 1,2,3,4,5 respectively. At further decrease 
of uj (T) the spectral curves diminish according to Eq.( (j3.88)l ) and attain the minimum 
at uJmin = u t h/{k{l + 2/6)) (T(l + 2/6) = 1). The corresponding values (in MeV) are 
uj ~ 3.5,2,1.4,1.2 for curves 1,2,3,4. The least value of uj m i n ~ 1 MeV has the curve 
5. However, one has to take into account that at uj < 1.5 MeV (kq = p 2 c ~ 2/6) the 
contribution of the transition radiation becomes significant. Starting from uj < 0.6 MeV 
the contribution of the transition radiation dominates. 
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The approach presented in this subsection is applicable in the interval of photon en- 
ergies where effects of the polarization of a medium are essential. It includes also the 
soft part of the LPM effect. For the case given in Fig. 11 our results are given up to 
Umax ~ 20 MeV. On the other hand, in [UT] the hard part of the LPM effect spectrum was 
analyzed where one can neglect the effects of the polarization of a medium (u > 5 MeV 
for the mentioned case). Although here (as in [20]) in and in [HT] the different methods 
are used, the results obtained in overlapping regions are in a quite reasonable agreement 
among themselves. 

It is interesting to discuss behavior of the spectral curves obtained in [UT] from the 
point of view of our results. We consider the low value l\ (T c is not very large) and a 
situation when corrections to the value b in Eq. ()3.94|) which neglected in [61] are less 
than 20%. This leads to the difference in results less than 10%. We concentrate on the 
case of gold target with the total thickness Nli = 0.7% L ra( i. The case N — 1 where 
T c = 5.8, b~ l = 7.3, tu c ~ 240 Mev, u th = Au c /(T C (1 + T c )) ~ 24 MeV was considered 
in detail in [UT] (see Sec. 2). The curves in figures in [UT] are normalized on the Bethe- 
Heitler probability of radiation (no the Coulomb corrections), i.e. they measured in units 
cer2T c / (3ir) . In the region where our results are applicable uj < uo t h — 24 MeV in Fig. 2 
of [HI] (G = 0) one can see plateau the ordinate of which is 10% less than calculated 
according with Eq. (|3.82|) . The case of two plates (T c ~ b^ 1 ~ 3, uo c ~ 240 MeV, u th ~ 
w c /T c = 80 MeV) is given in Fig. 3 of [fJT]. The lengths of the gaps are the same as 
in our Fig. 11, except k = 9. The positions and ordinates of the minimums and the 
maximums in the characteristic points (u ~ nuth/k for minimums and uj ~ ooth/k (T = 1) 
for maximums, see above) as well as behavior of the spectral curves is described quite 
satisfactory by our formulas (see e.g. asymptotic Eqs. ([3.87|) - ([3.90p ). 

3.7 Comparison of the results obtained in different papers 

As it was said above, Migdal have calculated the probability of the bremsstrahlung and 
pair creation in the logarithmic approximation 2j.This result, now considered as classic 
one, was confirmed in the most of subsequent calculation. 

We developed the method of calculation [T7|[TBl[IHll2n] which goes beyond the logarith- 
mic approximation. The formula ()2.17|) gives the spectral distribution of the probability 
of radiation which coincides formally with the probability calculated by Migdal (Eq.(49) 
in |2J). However, here Coulomb corrections are included into parameter v in contrast to 
j2], As it was noted, the parameter g c entering into the parameter v (see Eqs. (|2.16j) . 
(I2.12|) and (|2.1U|) ) is defined up to the factor ~ 1, what is inherent in the logarithmic 
approximation. However, we calculated also the next term of the decomposition over 
v(g) (an accuracy up to the "next to leading logarithm") and this permits to obtain the 
result which is independent of the parameter g c . It was shown that the definition of the 
parameter g c minimizes corrections to Eqs. (j2.17|) and (|2.43J) practically for all values of 
the parameter g c . It should be emphasized also that here the Coulomb corrections are 
included into the parameter v in contrast to j^j. Moreover, the procedure of fine tuning 
presented above (see [T7] ) permits one (after adding the term dWi/du (I2.22J1 ) to obtain 
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more accurate expression for the probability of radiation, which has power accuracy in the 
limit of weak multiple scattering, see ([2.29)1 . while the probability calculated by Migdal 
[2] has logarithmic accuracy only in the mentioned limit. The interrelation between the 
main term dW c /(ko and the first correction dWi/du depends essentially on fine tuning 
procedure which was chosen, as it was just said, in a such way that minimize the value of 
the first correction. For example, one can see that the maximal value of dWi/du is of the 
order 7 % of the main term for all four spectral curves shown in the Fig. 3. This means 
that in our case there are more rigid requirements on fine tuning procedure than used by 
Migdal. 

In the approach used by Blankenbecler and Drell (BD) [22| the scattering process 
was considered in the second approximation of the eikonal formalism and an analysis 
of scattering in medium was based on the model of a random medium. The radiation 
was calculated for all possible paths and averaged. One advantage of this approach is 
that it treats naturally finite thickness targets. There are a few short-comings in the BD 
approach. The Gaussian distribution of scattered electrons was obtained by Migdal [2] as 
a solution of the Fokker-Planck equation, while in BD approach this is an input. Actually 
there are deviations from the Gaussian distribution, see end of Sec. 3. 3. The calculation of 
the probability dW/ dw is quite unwieldy in BD approach and general formula for dW/ duo 
(Eq.(5.32) in [22])is written down in an implicit form only, and the functions r](z2,Zi,l) 
and A (,22, z\, I) in the mentioned formula have to be calculated in the each case separately. 
This means that comparison of BD probability with that calculated by Migdal can be 
done in limiting cases only. In the case of weak multiple scattering the probability in [22] 
coincides with the Bethe-Heitler probability taken with logarithmic accuracy (just as in 
Migdal paper [2J). However, in the limit of the strong LPM effect, where the probability 
of radiation is given by Eqs. (|2.17|) and (|2.31j) above and which coincides with obtained by 
Migdal Eq.(52) in j2] (if one omits the Coulomb corrections included into the parameter 
z/ ), the probability calculated by BD (the term with R 2 ) differs from (j2.31|) by factor 
y^37r/8=1.0854. The origin of this discrepancy was not analyzed in |22j . 

In the paper of R.Baier, Dokshitzer, Mueller, Peigne, and Schiff (BDMPS) [2H] the 
multiple scattering of high-energy electrons off a large number of scatterers (presented as 
a screened Coulomb potential) is analyzed within eikonal formalism. Radiation of soft 
photons only is included using the classical current approach, although this approach 
gives no significant simplification for the problem under consideration. All the results 
are obtained within logarithmic accuracy. No general formula of the type of Eq. (|2.18j) is 
written down. In the "Bethe-Heitler limit" the result of [22] does not completely match 
the standard Bethe-Heitler formula. Since the method used consists in direct summation 
over scatterers it naturally included the consideration of a target of finite thickness. For 
thin target (T < lf ) the result of [25 j is similar to logarithmic term in Eq. ()3.18J) where 
r 2 = 2 for soft photons. The explicit formula for radiation spectrum is given in [23] in 
the limit of strong LPM effect where it is similar to Eq. 1)2.32)1 if uj <C e and the Coulomb 
corrections are omitted. 

Note that in Migdal papers [2], and in our papers the kinetic equation is the basic 
element of the used methods while in the both papers |25^ and [22] the consideration 
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does not included the kinetic equation. Due to this feature the calculations in the both 
papers [22] and [22| are more cumbersome and no general formulas for radiation spectra 
are obtained. 

A path integral method is applied successfully for analysis of multiple scattering (see 
e.g. [63j). Zakharov (221, [211 is used this method to study of the LPM effect. Using 
a transverse Green function based on a path integral he succeeded in reduction of the 
problem to the two dimensional Schrodinger equation, i.e. he re-derived equation a la 
Eq.(J22J), which was obtained earlier in 0. The probability of radiation was calculated 
(see Eq.(19) in (231) within logarithmic accuracy (as in Migdal papers [2j and [3J) although 
it is written down in a more cumbersome form. The procedure of fine tuning so important 
for Migdal calculation [2] is oversimplified in [2B]. For earlier attempt to do study of the 
LPM effect in a such way see jHSl- The formalism used in j2H| has the same region 
of validity as in Migdal' s papers [2], [3] and in our papers 0, [T7j and [TH] in spite of 
claims about "rigorous treatment of the LPM effect". The developed formalism allows 
naturally to consider the LPM effect in finite thickness targets [25- As an example the 
calculation of spectral intensity of radiation is given for electron energy e = 25 GeV and 
gold target with thickness I = 0.7%L ra d. This calculation was carried out numerically 
with the complete potential (see Eq.(j2HU)) outside of scope of logarithmic approximation. 
The Coulomb corrections were discarded as well as a multiple photon radiation. After 
arbitrary diminishing of calculated value by 7% it was found very good agreement with 
a portion of data (in the region of the LPM drop up to beginning of plateau). Actually 
this diminishing should be attributed to Coulomb corrections found in ^7] (for details see 
[T%] and below in Sec. 3.10). More late results of Zakharov see in [64] . 



3.8 Qualitative behavior of the spectral intensity of radiation 

We consider the spectral intensity of radiation for the energy of the initial electrons when 
the LPM suppression of the intensity of radiation takes place for relatively soft energies 
of photons: u < u c <C e: 



16nZ 2 



a 



u iu c 



UJr 



l 2 n a In 



a S 2 
A ' 



(3.97) 



see Eos.([2~3|). (j2H21), (I2~2lj) - (J2~2"TJ) . This situation corresponds to the 
experimental conditions [121 EH EU EES] • 

The ratio of the thickness of the target I and the formation length of radiation 1/(uj) 
(ll.5j) is an important characteristics of the process. In Eq. (jl.5j) both the multiple scatter- 
ing and the polarization of a medium are taken into account. This ratio may be written 

as 
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where we put that vq ~ J — . Below we assume that uj c ^> uj p which is true under the 

experimental conditions. 

If /3(uj c ) ^> 1 (T c ^> 1) then at u > uj c a target is thick and one has the LPM 
suppression for to < u c . There are two opportunities depending on the minimal value of 
the parameter (3: 

^i = ^(^) 1/3 , P m ^2T c (^y /3 . (3.99) 

If j3 m <C 1 then for photon energies uj > u\ it will be u>b such that 

(3(u b ) = l, u b ~^ (3.100) 

and for uj < u>b the thickness of a target becomes smaller than the formation length of 
radiation so that for uj <C ujj, the spectral distribution of the radiation intensity is described 
by Eqs. (|3.36j) . (|3.40|) . Under these conditions for 4k = v\T = T c ^> 1 the spectral curve 
has a plateau 

dl 2aJ 

— = = const 3.101 

dU 7T 

in accordance with Eq. (j3.20|) . Under conditions kT <C 1, w < UJ}, the spectral intensity of 
radiation is independent of photon energy w. It should be noted that due to smallness of 
the coefficients in expression for 5 Eq. (j3.38|) . such behaviors of the spectral curve begins 
at u < uJth = 4co>b ~ 4u; c /T c 2 . It will continue until photon energies where one has to take 
into account the polarization of a medium and connected with it a contribution of the 
transition radiation. 

At (3 m >la target remains thick for all photon energies and radiation is described by 
formulae of Sections 2.2 and 2.6. In this case at uj <^ lu c (u ^> 1) and uj ^> (ujp/ujc) 1 ^ 3 uj p 
(u ^> k) the spectral intensity of radiation formed inside a target is given by Eqs. (|2.17|) 
and (|2.22|) and the contribution of the boundary photons is given by Eq. (|3.12| ). 

It is important to include also the contribution of boundary photons. Since the con- 
tribution into the spectral intensity of radiation from a passage of the electron inside the 
target (oc T) is diminishing and a contribution of the boundary photons is increasing 
with uj decrease, the spectral curve has a minimum at u — u m . The value of uj m may be 
estimated from equations (see (|2.31|) . ()3.36|) - ([3.37|) )) 

~T — ^ + In f H 7T A =0, —= ~ 1 + 



duj \ y/2 24 u Q J ' y/2 4:V2u ' 

x = ^. (3.102) 




UJ 



When the value of T c is high enough, the solution of Eq. (j3.102|) does not satisfy the con- 
dition z/q 3> k and in this case the equation ()3.102|) ceases to be valid. For determination 
of u) m in this case we use the behavior of the spectral intensity of radiation at k ^> u . 
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In this case the contribution into radiation from inside passage of the target is described 
by Eq. (|2.124|) while the radiation of the boundary photons reduces to the transition ra- 
diation and its contribution is given by Eq. (j3.10|) . Leaving the dominant terms (vqT is to 
independent) we have 

'tfir , \ v l T T c IT , 

+ hiK =0, — — = 1, K m = — , Um ~J— oj p . (3.103) 

J OK 3 \ J-c 

Since the value 7r 2 /12 ~ 0.8 is of the order of unity, the solution of (|3.102j) at n m ^> 
vq differs only slightly from uj m . Because of this, if the condition 2T c (uj p /uj c ) 2 ^ ^> 1 is 
fulfilled, the position of the minimum is defined by Eq. (j3.102|) . The formulas (|3. 102(1 and 
(I3.103() show the position of a minimum of the spectral intensity curve which is the sum 
of contributions of radiation inside a target and boundary radiation. Due to the LPM 
effect the intensity of radiation inside a target diminishes with uj decreasing. However at 
low uj the contribution of the boundary radiation is sharply increasing and this results 
in appearance of a minimum of spectral curve. Eq. ()3. 102(1 gives its position at small T c 
while Eq. fl3.103() gives minimum position at large T c . 

If /3(u c ) = 2T C <C 1 then at u = uj c a target is thin and the Bethe-Maximon spectrum 

of radiation which is valid at uj ^> uj c ( — ^— ^ = const ) will be also valid at uj < uj c in 

\ aw J 

accordance with Eq. (|3.40(l since 4k = v^T — T c <C 1. This behavior of the spectral curve 
will continue with uj decrease until photon energies where a contribution of the transition 
radiation become essential. In this case the spectral distribution of radiation has the form 
flUH) for all uj 

dw dw tr , dlU^M /„ 1nyl x 

+ cos(kT) — - — , (3.104) 



duj duj duj 
Since for soft photons (uj e) 
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(3.105) 



and T c /3 C la contribution of the transition radiation become visible already at kT 1. 
For uj > uj c (T c <^ 1) the probability of radiation is defined by Eqs. ()3.32() - (j3.34() . In 
this case a considerable distinction from Bethe-Maximon formula will be in the region 

uj ~ uj c /T c . 

As well known, for soft photons (uj e) the Bethe-Maximon formula for the spectral 
intensity of radiation doesn't depend on a photon energy as well. So, the ratio of these 
spectral intensities (see Eqs. (|3.40() and (|3.43(0 is an important characteristics of the 
phenomenon under consideration: 







For T c ^ 1 one has R = 1 and for T c ^> 1 one has using expression 1)3.18(1 

R ^¥ c ( 1 + ^k) [lnAk+1 ~ c] ~ 2 + T; (3 - 107) 
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For estimates one can put with a good accuracy 4k ~ T c and L t — L\. 

There is, in principle, an opportunity to measure the electron energy (in region of high 
energies) using the LPM effect. For this one can measure the spectral curve on a target 
with thickness a few percent of L ra( i and compare the result with the theory prediction. 

Existence of the plateau of the spectral curve in a region of photon energies where a 
target is thin was first found in jJ2] within Migdal approach (quantum theory). Recently 
this item was discussed in jSH] (in classical theory), [22] and [21]. 

3.9 Experimental investigation of the LPM effect 

Bremsstrahlung or pair creation suppression can be studied with high energy electron 
or photon beams. Because pair creation suppression requires photons with u > u e (see 
Eq. (j2.48p ). beyond the reach of current accelerators, pair creation has been studied only 
with cosmic rays, with consequently very limited statistics. The best suppression studies 
have used electron beams at accelerators. Besides the LPM effect, these beams have been 
used to study effect of polarization of a medium and boundary photon emission. 

The first tests of LPM suppression came shortly after Migdal's appeared. These were 
cosmic rays experiments with high energy photons (u > 1 TeV) and studied the depth 
of pair conversion in a dense target [HE]- [EH], see also more recent [70J. All of the air 
shower experiments suffered from the poor statistics. Uncertainties in the photon spec- 
trum complicated the analysis. So, these experiments are at best qualitative verification 
of the LPM effect. 

The first accelerator based study of LPM effect was done in Serpukhov [7T] using 
40 GeV electrons. Photons with energy 20 MeV < uj < 70 MeV emitted from carbon, 
aluminum, lead and tungsten targets were detected in an sodium iodide calorimeter. It 
were several limitations due to substantial background. The results agreed with Migdal's 
predictions. 

A precision measurement of LPM suppression, and study of effect of polarization of 
a medium was performed by the E-146 collaboration at SLAC using 25 GeV and 8 GeV 
electrons. The experimental setup was conceptually similar to the Serpukhov experiment, 
but heavily optimized to minimize background. The target used are listed in Table 1. 
Photons with energy 200 keV < u < 500 MeV were detected in a BGO calorimeter. The 
logarithmic binning of photon spectrum was used with 25 bins per decade. A very large 
data set was collected [12] - [IS]. A detailed, high statistics Monte Carlo simulation was 
used by E-146. Its main purpose was to understand multi-photon emission (see Sec. 3. 5). 

Recently new study of LPM effect at higher energies of electrons (e =149, 207 and 
287 GeV), where the effect has influence upon much wider part of spectrum comparing 
with e =25 GeV, was performed in the H2 beam line of the CERN SPS jTHj. The 
experimental setup was also conceptually similar to the Serpukhov experiment. Photons 
with energy 2 GeV < uo < e were detected in a lead glass calorimeter. The logarithmic 
binning of photon spectrum was used with 25 bins per decade as at SLAC. 
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3.10 Discussion of theory and experiment 
3.10.1 SLAC E-146 

Here we consider the experimental data from the point of view of the above 

analysis. It was shown that the mechanism of radiation depends strongly on the thickness 
of the target. First, we estimate the thickness of used target in terms of the formation 
length. From Eq. ()3.98j) we have that 

T c = > 20 at — — > 2 %. (3.108) 

The minimal value of the ratio of the thickness of a target to the formation length is 
given by Eq. ()3.99j) (£L ~ 2T c (uj p /uj c ) 2 ^ 3 ). For defined value of T c this ratio is least of 
all for the heavy elements. Indeed, the value of u p = depends weakly on nucleus 

charge Z (cu = 30-r 80) eV, while co c — = — oc Z 2 . Furthermore, the ratio uj p /uj c decreases 

with the energy increase. Thus, among all the targets with thickness I > 2 %L ra d the 
minimal value of (3 m is attained for the heavy elements (Au, W, U) at the initial energy 
e = 25 GeV. In this case one has uj c ~ 250 MeV, uj p ~ 4 MeV, (3 m > 2.5. Since the 
parameter T c is energy independent and the ratio oj p /uj c oc 1/e, the minimal value f3 m > 5 
is attained at the initial energy e = 8 GeV for all the targets with thickness I > 2 %L ra d- 
So, all such targets can be considered as thick targets at both energies. 

As an example we calculated the spectrum of the energy losses in the tungsten target 
with the thickness I = 0.088 mm (=2.7 %L ra d) for both initial energies shown in Fig. 
12 (a) and (b). The characteristic parameters for this case are given in Table 2. We 
calculated the main (Migdal type) term Eq. (j2.17J) . the first correction term Eq. ()2.22|) 
taking into account an influence of the polarization of a medium according to Eq. (|2.123|) . 
as well as the Coulomb corrections entering into parameter i> Q Eq. (|2.24j) and value L(g c ) 
Eq. fl2.12j) . The contribution of an inelastic scattering of a projectile on atomic electrons 
was not included into the numerical calculation (it is ~ 1%), although it can be done 
using Eq. (j2.34ft . We calculated also the contribution of boundary photons (see Eq.(4.12) 
in ^7j). Here in the soft part of the spectrum u < uj^ (uj^ ~ 2 MeV for e = 25 GeV) the 
transition radiation term Eq. (J3.10)) dominates, while in the harder part of the boundary 
photon spectrum uo > Ud the terms depending on both the multiple scattering and the 
polarization of a medium give the contribution; for e = 8 GeV one has Ud — 700 keV. 
All the mentioned contribution presented separately in Fig. 12. Under conditions of the 
experiment the multiphoton reduction of the spectral curve is very essential. The curve 
"T" in the Fig. 12 involves the reduction factor which was constructed as the interpolation 
of Eqs. (j3.55|) . (|3.57j) and (|3.63|) with the accuracy up to 1%. The curves "T" in Fig. 12 
(a) and (b) are the final theory prediction in units 2a /tc. 

Experimental data are taken from [14^ and recalculated according with given in men- 
tioned paper procedure: 

' N ^^[-n—\ ■ (3-109) 




Lrad k '' \dlauj , 

exp ' "** \ / exp 
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Since photon energies were histogrammed logarithmically, using 25 bins per decade of 
energy, one has for the coefficient k 



k h = = e — f = 0.096, 

k m = 2 Umax - Umin = = 0.092, 

ki = = 1 - e = 0.88, 



'max 

s = l n ^ = = 0.092, (3.110) 

25 v 



depending on the normalization point within bin. In the review jToj it was indicated that 
k = 0.096 = kh- In more earlier paper [Hj the value = 0.09 was given and we used it in 
our paper [T^J and [IE]. The date recalculated with k = 0.096 are also given in Fig. 12. 
It is seen that there is a perfect agreement of the curves T with data for both energies. 

The gold targets with thickness I = 0.7%L ra d and I = 0.1%L ra d are an exception and 
can't be considered as thick targets. We calculated energy losses spectra in these targets 
for both the initial electron energies e = 25 GeV and e = 8 GeV. The characteristic 
parameters of radiation for these cases are given in Table 3. 

In Fig. 13(a) results of calculations are given for target with a thickness 
/ = 0.7%L rad at £ = 25 GeV. The curves 1,2,3,4 present correspondingly the func- 
tions j} 2) , J 2 (2) , jf\ j| 2) (EH) . At u = 500 MeV the value z^T = J—T c = 8.4 > 1, 



the order 10% according with asymptotic damping factor 1 



the interference terms are exponentially small and one can use formulae for a thick 
target. In this case the parameter v\ = 0.69 and contribution of boundary photons 

(J b = j[ 2) + jf ] + jf >) is small (J 6 ~ -~2i> see E> Eq.(4.16)) and distinction Bethe- 

Maximon formula (J BM = T c /3 = 1.94) from Re = Re (j{ 2) + + jj 2) + jf ] ) is of 

16^ \ 

At uj < ujth — 30 MeV for the case T c 3> 1 and j3 m < 1 the spectral curve turns into 
plateau according with discussion in previous subsection. In this photon energy region 
the parameter vq > 3 and Eq. (j3.36J) for a target with intermediate thickness describes the 
spectral probability of radiation with a good accuracy. With the further photon energy 
decrease one can use limiting formula ()3.38|) where UqT = 4k ~ 7.4. For this case the ratio 
of ordinate of the plateau to ordinate of the Bethe-Maximon intensity (see Eq. (|3.107jl ) is 
R ~ 0.57. Note, that in formulae for j[ 2 ^ j] 2 -* the potential V c (q) Eq. ()2.12|) is used which 
doesn't include corrections ~ 1/L (y(g)). These corrections were given above for both 
thin and thick targets (see Eq. (|2.33|) ). In our case (u ^> 1, u T ^> 1) the expressions 
with corrections ~ 1/L are given in Eqs. (j2.33J) and (|3.41|) for a thick target and a thin 
target respectively. Taking into account behaviors of correction in the region Uq < 1 (see 
curve 2 in Fig. 12(a)) we construct an interpolation factor (taking into account the term 
~ 1/L) with accuracy of order 1%. The summary curve (T) in Fig. 13(a) contains this 
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factor. 

The transition radiation contributes in the region on < oo p (function Re curve 5 

f 2) _____ 

in Fig. 13(a)). When kT <C 1 this curve is described by asymptotic of Re Jg Eq. (j3.36j) . 
The contribution of the transition radiation increases with u decrease and for kT ~ 1 
it describes by Eq. (j3.44j) . The contribution of the multiple scattering diminishes due to 
interference factor cos(kT) in Eq. (j3.46J) (at uj = 0.2 MeV, kT = 1.9). The curve T in 
Fig. 13(a) gives the summary contribution of the multiple scattering (the curve S, where 
factor (1 — Ld/e) is included) and the transition radiation (the curve 5). 

In Fig. 13(b) results of calculations are given for target with a thickness 
/ = 0.7%L ra d at e = 8 GeV. The notations are the same as in Fig. 13(a). In this case the 
characteristic photon energy uj c is one order of magnitude lower than for e = 25 GeV, so 
that at uj = 500 MeV the parameter v\ is small iv\ ~ 1/20). Because of this the right 
part of the curve 5* coincides with a good accuracy with the Bethe-Maximon formula 
(the Coulomb corrections are included). Note, that for this electron energy the effect of 
recoil (factor (1 — uj/e)) is more essential. Strictly speaking, a target with a thickness 
0.7%L ra d at e = 8 GeV is not thin target for any photon energy {(3 m = 1.6). However, for 
bremsstrahlung this target can be considered as a thin one for uj < uj t h = 3 MeV . Since 
the polarization of a medium becomes essential in the same region (u p = 1.25 MeV), the 
interference factor cos(ftT) in Eq. (|3.46|) causes an inflection of the spectral curve S at 
uj ~ 1 MeV. The transition radiation grows from the same photon energy uj and because 
of this the total spectral curve T has a minimum at uj ~ 1 MeV . As far as there is some 
interval of energies between oj p and oj t h (&th — cu p ~ 3 MeV), this minimum is enough 
wide. Moreover, the value of its ordinate coincide with a good accuracy with ordinate of 
the plateau of the spectral curve S in Fig. 18(a) because bremsstrahlung on a thin target 
is independent of electron energy (J3.41)) .It should be mentioned also that for this thickness 
the reduction factor due to multiple photon emission (see Eq. (j3.57|) ) /lpm — 0.94 (so this 
is the effect of the same scale) but this reduction is compensated nearly exactly because 
of more precise definition of data (change of a bin width Ak/k = 0.09 in J2| (used in 
|21) to a bin width Ak/k = 0.096 in [H], [T3], see discussion in Sec.3 of [□]). 

In Fig. 14(a) results of calculations are given for target with a thickness 
0.1%L ra d at e = 25 GeV. The notations are the same as in Fig. 13(a). For this thickness 
T c = 0.96 and one has a thin target starting from uj < uj c . So, we have here very wide 
plateau. The left edge of the plateau is defined by the contribution of transition radiation 
(uj ~ uj p ). Since in this case 4k = v\T = T c ~ 1 (see Eqs. (l3.40j) - (j3.43|) ). one has to 
calculate the ordinate of the plateau using the exact formula for a thin target (J3.40|) . For 
this case the ratio of ordinate of the plateau to ordinate of the Bethe-Maximon intensity 
(see Eq. (j3.107J) ) is R ~ 0.85. The same ordinate has the plateau for electron energy 
e = 8 GeV (Fig. 13(b)). However, a width of the plateau for this electron energy is 
more narrow (1 MeV -j- 20 MeV) due to diminishing of the interval between u p and 
uj c . For uj > uj c the formation length of radiation becomes shorter than target thickness 
(T = T c uj/uj c > 1) and the parameter z/ x decreases. The value v{T = T c ^Ju/u c increases 
with uj growth. A target becomes thick and the spectral curves is described by the Bethe- 
Maximon formula in Fig. 14(b) starting from photon energy uj ~ 100 MeV. In Fig. 
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14(b) the contributions of separate terms into Re are shown as well. Their behavior 
at u) > oj c is described quite satisfactory by Eqs. ()3.32J) - (j3.35|) (see also discussion at their 
derivation). 

We compared our calculations with experimental data ^3]. The curves T in Fig. 13, 
14 give theory prediction (no fitting parameters !) in units 2a /ix. It should be mentioned 
that in our papers [T7j and [THj we put fc=0.09 according with the instruction given in 
paper But the multiphoton effects were not included. However, as was stated above 
one has to use fc=0.096. So, from the one side for the gold target with the thickness 
I = 0.7 % L rad and energy e=25 GeV the reduction factor / ~ 0.94 for photons with 
energy ui < 10 MeV (plateau region). From the other side, use of the coefficient k = 
0.096 in (|3. 109(1 instead of k — 0.09 lowers data upon ~ 6% and this is imitate inclusion of 
the reduction factor. As a result, the excellent agreement of the theory and data noted in 
[TS] is not broken. It is seen that in Fig. 13(a) there is a perfect agreement of the theory 
and data. In Fig. 13(b) there is an overall difference: data is order of 10% higher than 
theory curve. For photon energy uj = 500 MeV the theory coincides with Bethe-Maximon 
formula (with the Coulomb corrections) applicable for this energy. Note that just for this 
case it was similar problem with normalization of data matching with the Migdal Monte 
Carlo simulation (+12.2%, see Table II in jTJj). 

For thickness I = 0.1% L rad there is a qualitative difference between our theory pre- 
diction and Monte Carlo simulation in [T3]. There was a number of experimental uncer- 
tainties associated with this target. Nevertheless, we show data for e = 25 GeV which 
are lying higher than theory curve. 

3.10.2 CERN SPS experiment 

We consider now the data JE] from the point of view of the above analysis Energy 
loss spectra were measured in iridium target of 0.128 mm thickness (4.36 % L ra( i). Basing 
on Eq. (j3 .108(1 and following analysis we see that this is a thick target. Since the photons 
with the energy u > 2 GeV were measured, no boundary effects were observed. For Ir the 
characteristic electron energy Eq. ((2.26(1 e e =2.27 TeV and the characteristic photon energy 
Ea. ((2~2T() u c (e) for which the LPM effect well manifests itself is u c (287 GeV) =32 GeV. 
The result of calculation for the initial electron energy e=287 GeV is shown in Fig. 
15. The curves 1, 2, 3 are calculated using Eqs. ([2.29() . ((2.22(1 . ((2.17(1 correspondingly. 
It should be mentioned that the prediction of our theory (curve 4) in the hard end of 
spectrum coincides with the Bethe-Maximon curve within the accuracy better than 10~ 3 . 
For used thickness of target the multi-photon effects are very essential. In this case it 
is convenient to calculate the reduction factor / using Eq. (j3.73j) . where the mentioned 
spectrum Eq. ((2.17(1 is used. The result obtained is in a good agreement with Eqs. ((3.54() . 
(13.57(1 . The reduction factor for this energy is presented by curve 3 in Fig. 18. The final 
prediction with the reduction factor taken into account for the used target thickness is 
presented by the curve T. The data are recalculated according with Eqs. ((3. 109(1 . (|3. 110(1 
using the coefficient k m . The result of calculation for the initial electron energy e=207 GeV 
is shown in Fig. 16 and the result of calculation for the initial electron energy e=149 GeV 
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is shown in Fig. 17. The reduction factors for these energies are presented by the curves 
2 and 1 in Fig. 18 correspondingly. From Fig. 18 it follows that if the electron energy 
decreases the reduction factor / diminishes in agreement with Eqs. (j3.54j) . ([3.57)1 . It is 
seen that for energy 287 GeV there is the quite satisfactory agreement of theory with 
data, for energy 207 GeV the agreement is somewhat less satisfactory, and for energy 
149 GeV data are 10-15% below the theory. The difference is the same in hard part of 
spectrum {oj > 40 GeV) where the Bethe-Maximon intensity spectrum is valid. 

The decrease of relative energy losses of electron — in consequence of the LPM effect 

£ 

can be important in electromagnetic calorimeters operating in detectors on colliders in 

TeV range |21j . This item was discussed in Sec. 2. 2 using the main (Migdal) term, —L° rad 

in gold is given in Fig. 4 (curve 1). An analysis with better accuracy is of evident 
interest. The relative energy loss of electron per unit time including the contribution 
of the correction term Eq. ([2.17|) in terms of the Bethe-Maximon radiation length L ra d 
Eq. (j2.36J) for iridium (e e =2.27 TeV) curve 1 and for lead (e e =4.385 TeV) curve 2 is shown 
in Fig. 19. The increase of effective radiation length is 13% in Ir and 7.5% in Pb for the 
electron energy 1 TeV, is 38% in Ir and 24.5% in Pb for the electron energy 5 TeV, and is 
58% in Ir and 38.5% in Pb for the electron energy 10 TeV. For both media the corrections 
attain the maximal value of order 5% at energy e ~ e e . 



4 Effects in colliding electron-positron beams 
4.1 Mechanisms of radiation 

In this section we consider the radiation at head-on collision of high energy electron and 
positron beams. The properties of photon emission process from a particle are immedi- 
ately connected with details of its motion. It is convenient to consider the motion and 
radiation from particles of one beam in the rest frame of other beam (the target beam). 
In this case the target beam is an ensemble of the Coulomb centers. The radiation takes 
place at scattering of a particle from these centers. If the target consists of neutral par- 
ticles forming an amorphous medium, a velocity of particle changes (in a random way) 
only at small impact distances because of screening. In the radiation theory just the 
random collisions are the mechanism which leads to the incoherent radiation. For collid- 
ing beams significant contributions into radiation give the large impact parameters (very 
small momentum transfers) due to the long-range character of the Coulomb forces. As 
a result, in the interaction volume, which is determined also by the formation length If 
Eqs. fll.5|) . (jl.6|) (in the longitudinal direction), it may be large number of target particles. 
Let us note that in the case when the contribution into the radiation is given by impact 
parameters comparable with the transverse size of target, the number of particles in the 
interaction volume is determined by the ratio of the radiation formation length to the 
mean longitudinal distance between particles. 

However, not all cases of momentum transfer should be interpreted as a result of 
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random collisions. One have to exclude the collisions, which are macroscopic certain 
events. For elaboration of such exclusion we present the exact microscopic momentum 
transfer to the target particle in the form: q =< q > +q s . Here < q > is the mean value of 
momentum transfer calculated according to standard macroscopic electrodynamics rules 
with averaging over domains containing many particles. The longitudinal size of these 
domains should be large with respect to longitudinal distances between target particles 
and simultaneously small with respect to the radiation formation length. The motion of 
particle in the averaged potential of target beam, which corresponds to the momentum 
transfer < q >, determines the coherent radiation (the beamstrahlung) . While the term 
q s describes the random collisions which define the process of incoherent radiation (the 
bremsstrahlung). Such random collisions we will call "scattering" since < q s >= 0. 

One of principal characteristics of particle motion defining the properties of coherent 
radiation is the ratio of variation of its transverse momentum to the mass during the 
whole time of passage across the opposite beam T 



Ap ± eE ± a z 2aN c \ c 

o, (4.1) 



m m a x + a y 



where N c is the number of particles in the opposite beam, a x and a y is its transverse 
dimensions (a y < a x ), a z is the longitudinal size of opposite beam. The dispersion of 
particle momentum during time T is small comparing with m. It attains the maximum 
for the coaxial beams: 

<3| = 72 /,A.^I i<<1 , (4 . 2) 

m 2 \ I a x a y 

here f> 2 s is the square of mean angle of multiple scattering, L is the characteristic logarithm 
of scattering problem (L ~ 10). This inequality permits one to use the perturbation theory 
for consideration of bremsstrahlung, and to analyze the beamstrahlung independently 
from the bremsstrahlung 1 . 

Another important characteristics of motion is the relative variation of particle impact 
parameter during time T 

Agi eE ± a 2 z 2aN c \ c a z = ^ ^ 



l(&x + <Jy)(Ji 

here i is x or y. When the disruption parameter Di 1, the collision doesn't change 
the beam configuration and the particle crosses the opposite beam on the fixed impact 
parameter. If in addition the parameter 5 1 (this situation is realized in colliders with 
relatively low energies) then the beamstrahlung process can be calculated using the dipole 
approximation. The main contribution into the beamstrahlung give soft photons with an 



1 Actually more soft condition should be fulfilled: 

(q 2 s (l f )) /m 2 = (q2) /m 2 l f /a z «l 
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energy 

^<^«1. (4.4) 

e a z 

In the opposite case 5 ^> 1 the main part of beamstrahlung is formed when the angle 
of deflection of particle velocity is of the order of characteristic radiation angle I/7 and 
the radiation formation length l m is defined by 



eE\l m 2aN r XJ m a 2 . . 

1, l m = — ; (4.5) 



r^j 



m (a x + <j y )a z "" 5 

and the characteristic photon energy is 

7 2 A 2 
u; ~ Um = — = exm, Xm = 2aN c j- 6 (x m < I)- (4.6) 

tm (fx + 0- y )a z 

Here x is the invariant parameter Eq. (|1.21|) which defines properties of magnetic bremsstrahlunj 
in the constant field approximation (CFA). For applicability of CFA it is necessary that 
relative variation of Ej_ was small on the radiation formation length l m . As far l m is 
shorter than a z in 5 ^> 1 times the characteristic parameter becomes 

A™ = A- = ^ = — , (i = *,y) (4.7) 
a z b ^Oi 

to that extent. The condition D m i <C 1 is fulfilled in all known cases. The mean number of 
photons emitted by a particle during the whole time of passage across the opposite beam 
T is iV 7 ~ a5, it include the electromagnetic interaction constant. Using the estimate 
(14. 6|) we get an estimate of relative energy loss 

~ ab Xm (x™ < 1) (4.8) 

In the case Xm 1 (this condition is satisfied in all existing facilities and proposed collider 
projects) the soft photons with energy to ~ u m = ex m *C e are mainly emitted. For uj 
oj m the emission probability is exponentially suppressed. So, such photons are emitted in 
the bremsstrahlung process only. The boundary photon energy Ub, starting from which 
the bremsstrahlung process dominates, depends on particular parameters of facility. If 
Xm ~ 1/10 the energy is u>t, ~ e. The formation length for u 3> u m is much shorter 
than l m . On this length the particle deflection angle is small comparing with I/7 and 
one can neglect the variation of transverse beam dimensions (see Eq. (|4.7[l ). This means 
that all calculations of bremsstrahlung characteristics can be carried out in adiabatic 
approximation using local beam characteristics a Xty (t),v(t) etc, with subsequent averaging 
of radiation characteristics over time. Note that actually we performed a covariant analysis 
and the characteristic parameters are defined in a laboratory frame. 

As an example we consider the situation when the configuration of beams doesn't 
change during the beam collision (the disruption parameter D 1), and the total particle 
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deflection angle during intersection of whole beam is small comparing with the charac- 
teristic radiation angle I/7 (the dipole case). The target beam in its rest frame is the 
ensemble of classical potentials centers with coordinates r a (x a , z a ) and the transverse 
coordinate of emitting particle is rj_. In the perturbation theory the total matrix element 
of the radiation process can be written as 



M(r ± ) = J2 m ( r -L - x a) exp(iqnz a ) 



0=1 



We represent the combination MjM* in the form 



(4.9) 



x a )m 7 (r_L - x 6 j 



M i Mj = y £m i (T ± 

a=b 

+ Xl m *( r ^ _ ^a)mj(r± - x 6 ) exp(-iq\\ (z a 



(4.10) 



In the expression Eq. (|4.1U|) we have to carry out averaging over position of scattering 
centers. We will proceed under assumption that there are many scattering centers within 
the radiation formation length If = l/q\\ 



N f 



n z lf > 1, 



where for the Gaussian distribution 



n. 



2na. 



exp 



2a 2 



(4.11) 



(4.12) 



here N c and a z are introduced in Eq. (|4.1jl . Note that in the situation under consideration 
Qmax — — x a \max > o u where a t is the characteristic transverse size of target beam. 
Let us select terms with approximately fixed phase q\\(z a — Zb) = <t> a b in the sum with 
a ^ b in Eq. (j4.10|) . If the condition (|4.11|) is fulfilled, there are many terms for which the 
phase variation is small (A0 o j <C 1). For this reason one can average over the transverse 
coordinates (x a , x^) of target particles in Eq. (j4.10|) without touching upon the longitudinal 
coordinates (z a , Zb) 



Mi.MJ 



N c (m i m j ) ± + (rrii) ± (m j ) ± ^ exp(-iq\\ (z a - z b )) 

a^b 



N c ( (m i m j ) ± - (m i ) ± (m j ) ± ) + (m^ (mj)_ 



X]exp(ig||Z a ) 



(4.13) 



where 



(mi) ± = 
(m i m j ) ± 



J rrii(r± - x)w c (x)d 2 x, 

= J mj(r_|_ — x)mj(r± — x.)w c (x.)d 2 x, 



(4.14) 
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here w c (x.) is the probability density of target particle distribution over the transverse 
coordinates normalized to unity. In Eq. ()4.13|) in the sum with a 7^ b we add and subtract 
the terms with a = b. The first term (proportional to N c ) on the right-hand side of 
Eq. (J4.13|) is the incoherent contribution to radiation (the bremsstrahlung) . The second 
term gives the coherent part of radiation. For Gaussian distribution Eq. (j4.12j) performing 
averaging over the longitudinal coordinate z a one has 



4.2 The beam-size effect in the bremsstrahlung 

It was pointed in Sec. 1.3 that the external factors act differently on the radiating particle 
and on the recoil particle. This leads to a specific effect in the bremsstrahlung process at 
electron-electron(positron) collision. The point is that in the rest frame of recoil particle 
the formation length is longer by e/m times than in the laboratory frame. Just this length 
determines the maximal impact parameters at scattering (e.g. the screening radius). The 
are a few factors which could act on the recoil electron. One of them is the presence of 
an external magnetic field in the region of particle collision [27] (see Sec. 1.3). Another 
effect can appear due to smallness of the linear interval I where the collision occurs in 
comparison with l v Eq. (|1.6|) . This was pointed in [4"o] . 

The special experimental study of bremsstrahlung was performed at the electron- 
positron colliding beam facility VEPP-4 of Institute of Nuclear Physics, Novosibirsk [73] . 
The deviation of the bremsstrahlung spectrum from the standard QED spectrum was 
observed at the electron energy e = 1.84 GeV. The effect was attributed to the smallness 
of the transverse size of the colliding beams. In theory the problem was investigated in 
[74*] . where the bremsstrahlung spectrum at the collision of electron-electron (positron) 
beams with the small transverse size was calculated to within the power accuracy (the 
neglected terms are of the order I/7 = m/e). After the problem was analyzed in |75j . 
and later on in [TH] where some of results for the bremsstrahlung found in [71] were 
reproduced. It should be noted that in [71] (as well as in all other papers mentioned 
above) an incomplete expression for the bremsstrahlung spectrum was calculated. One 
has to perform the subtraction associated with the extraction of pure fluctuation process 
as it was shown above. An analogous problem encountered in the analysis of incoherent 
processes in the oriented crystals [EI],[E3] where the necessity of the subtraction procedure 
was indicated. Without the subtraction the results for the incoherent processes in oriented 
crystals would be qualitatively erroneous. 

The correction to photon emission probability due to the small transverse dimensions 
of colliding beam for unpolarized electrons and photon was calculated in [77] basing on 
subtraction procedure as in Eq. (j4.13|) . It is obtained after integration over the azimuth 
angle of the emitted photon 



2 00 



J2 ex v(iq\\z a ) -> n g exp(iq\\z)dz = iV c 2 exp(-gjc^) 




(4.15) 



a 3 e' duj 



U{QF(u,,C)dC, C = l + 7 2 <? 2 , 



(4.16) 
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where •& is the photon emission angle, e' = e — u, 

4(C - 1) . 



U(C) = v - 



C 2 



+ F(o;,C) = F«(a;,C)-F( 2 )(a;,C), 



V 



K ( V g)K 2 ( V g) - Kf(rjg) 



dg 



-d g 



F®(u>,Q = ^/ (J K 1 ( V g)^w c ^-Q)d 2 g ] j w r (x)d 2 x, 



here 



V 



qminC, Qmin = m 3 u/Ae 2 e', = J w r (x + g)w c (x.)d 2 x, 



(4.17) 



(4.1? 



where w c {x) is defined in (14. 14)1 . w r (x) is the same but for the radiating beam, value q min 
is defined in cm. frame of colliding particles. The term F^(cu, Q is the subtraction term. 
The total probability is dw~, = dwo + dwi, where dwo is standard QED probability. The 
analysis in [77] was based on Eqs. ()4.16j) - (j4.17)) . 

We considered in [77j the actual case of the Gaussian beams. The Fourier transform 
was used 



w x 



(2tt) 2 
w r (q) = exp 



<i 2 gexp(— iqx)w(q) 



-t(q 2 A 2 + q 2 A 2 ' 



w c {q) = exp 



/ 2 2 i 2 2 



(4.19) 



where as above the index r relates to the radiating beam and the index c relates to the 
target beam, A y and A x (a y and a x ) are the vertical and horizontal transverse dimensions 
of radiating (target) beam. Substituting Eq. ()4.19j l into Eq. (|4.18j) we find 



■exp[-Q x E x - g y E y \; S a 
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2(al + A 2 ) ' 



2(a 2 + Aiy 



(4.20) 



Using the relation do\ = $ 1 (0)dwi the following expression for the correction to spectrum 
was found in [77] starting from Eq. ()4.16|) 



da? = 



2a s e' duj 
m 2 e uj 
1 



/ (1) H, f(s) = —(v - 8s 2 )erfc(s) + 4e~ s ' + 2Ei(-s 2 ,, 
2s 



2?r 



nT.xT.y J £~ 2 cos 2 + S~ 2 sin </> 



F 2 {z)f(s)sds, 



(4.21) 



s 2 o" 2 



£~ 2 cos 2 <y9 + S~ 2 sin 2 <p ' 



2x 2 



x 



ln(x + VI + x 2 ) - 1, 



where Ei(a;) is the exponential integral function and erfc(x) is the error function. This 
formula is quite convenient for the numerical calculations. 
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The subtraction term (F^ 2 \uj,C) i n Eq. (j4.17j) ) gives for coaxial beams 



da\ 



(2) 



2cr e' duj 



(4.22) 



where 



ds\ \ ds2g 



G(si,s 2 ), 



G(s 1 ,s 2 ) 



1 / aia 2 bib 2 



1/2 



aia 2 &i& 2 



(4.23) 



Here the function (7(g) is: 



9(o) 



2g 2 



7T 

2<z~ 



g 2 ) erfc(g) + -e~ g + Ei(-g 2 



(4.24) 



In Eq. fl4.23j) we introduced the following notations 

1.1 1 



a 



2 A 2 ' 
A = a,\ + a 2 



2AI 



Ol,2 



1 



s lj2 + 2a2 
B = hi + b 2 + b, s = si + s 2 . 



Sl,2 + 2o* 



(4.25) 



In the case of narrow beams one has q m in/ (S x + E y ) <C 1. In this case of coaxial beams 



da~ = dcr + da \ is 



2a 3 e' duj 



2 

v 

3 



21n- 



C + 2 - J_ 



9 



where 



J_ 



dsi J ds 2 G(s 1 ,s 2 ) 



(4.26) 



(4.27) 



The dimensions of beams in the experiment [13] were <r y 



24 /im, a x 



A x = 450 /xm, so this is the case of flat beams. The estimate for this case gives 
J_ ~ (4/3\/3)ira y /a x <C 1. This term is much smaller than other terms in Eq. (|4.26|) . 
This means that for this case the correction to the spectrum calculated in [71] is very 
small. 

The result of calculation and VEPP-4 (INP, Novosibirsk) data are presented in Fig. 20 
where the bremsstrahlung intensity spectrum uda/dut is given in units 2ar 2 versus the 
photon energy in units of initial electron energy (x = oj/e) |77], The upper curve is the 
standard QED spectrum, the three close curves below are calculated using Eqs. (j4.21|) and 
(14.22)1 for the different vertical dimensions of colliding beams (equal for both colliding 
beams a = a z = A z ):a = 20 \ivn (bottom), a = 24 /im (middle), a = 27 jj,m (top) (this 
is just the 1-sigma dispersion for the beams used in the experiment). All the theoretical 
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curves are calculated to within the relativistic accuracy (the discarded terms are of the 
order m/e). It is seen that the effect of the small transverse dimensions is very essential 
in soft part of spectrum (at u/e — 10~ 4 the spectral curve diminishes in 25 %), while for 
uj/e > 1CT 1 the effect becomes negligible. The data measured in |73] are presented as 
circles (experiment in 1980) and as triangles (experiment in 1981) with 6 % systematic 
error as obtained in [7JJ| (while the statistical errors are negligible). This presentation is 
somewhat different from [73] • It is seen that the data points are situated systematically 
below the theory curves but the difference is not exceed the 2-sigma level [73]. It should 
be noted that this is true also in the hard part of spectrum where the beam-size effect is 
very small. 

One more measurement of beam-size effect was performed at the HERA electron- 
proton collider (DESY, Germany) |78j . The electron beam energy was e=27.5 GeV, the 
proton beam energy was e p =820 GeV. The parameters of beam in this experiment were 
(in our notation): o y = A y = (50 58)/zm, a x = A x = (250 -r- 290)/im. Since the 
ratio of the vertical and the horizontal dimensions is not very small, the contribution of 
subtraction term fEq. (J4.22j) ) is essential (more than 10%). The bremsstrahlung intensity 
spectrum uoda/duj in units 2oTq versus the photon energy in the units of initial electron 
energy (x = uj/e) for the HERA experiment is given in Fig. 21. For details of comparison 
of experimental data [73] , [THj with theory see [77] , where we discussed also possible use 
of beam-size effect for linear collider tuning. It should be noted that for linear collider 
the condition of strong beam-size effect a y q min <C 1 {o y <C <j x ) is fulfilled for the whole 
spectrum. This can be seen in Fig. 22, where the lower curve is calculated using Eq. (j4.21J) 
and the subtraction term is very small since o y /a x < 0.01. As far as the narrow beams are 
considered in Fig. 22, the lower curve is consistent also with Eq. (j4.26|) . This curve depends 
on the energy and the transverse sizes of beams. It will be instructive to remind that 
the analysis in [77] (see Eq.(2.8)) and here is valid if Xm/u <C 1 (see Eq. (j4.6j) . u = u/e'). 
The parameter Xm depends also on number of particles N c and the longitudinal beam 
size. So, for low N c Fig. 22 is valid for any x, but for TESLA project (x m =0.13) it holds 
in hard part of spectrum only. In fact, the probability of incoherent radiation becomes 
larger than the probability of coherent radiation only at x > 0.7 where the lower curve in 
Fig. 22 is certainly applicable. 

4.3 Coherent radiation 

The particle interaction at beam-beam collision in linear colliders occurs in an electromag- 
netic field provided by the beams. As a result, l)the phenomena induced by this field turns 
out to be very essential, 2)the cross section of the main QED processes are modified com- 
paring to the case of free particles. These items were considered by V.M.Strahkhovenko 
and authors [79*] . 

The magnetic bremsstrahlung mechanism dominates and its characteristics are deter- 
mined by the value of the quantum parameter xif) Eq. (|1.21|) dependent on the strength 
of the incoming beam field at the moment t (the constant field limit (CFA)). 

In the CFA the spectral probability of radiation from an electron per unit time has 
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the form 0, Eq. (10.25) (see also Eq.(4.24) in QUI or Eq. (90.23) in [TI] ) 

^ 7 " dW^t) = " a J ^(t)dt<ku, 



dt 7V 7 2^72 

&r(t) = + 7) 'M*) - / KyMdy, (4.28) 



-00 

00 



where K v (z) is the Macdonald functions, 2 = 2«/3x(t), it = ui/e'. The corresponding 
probability for polarized electrons and photons is given in 
For the Gaussian beams 

X(<) = Xo(x,y) ex V (-2t 2 /cr 2 z ), (4-29) 

here the function Xo( x iU) depends on transverse coordinates. 

It turns out that for the Gaussian beams the integration of the spectral probability 
over time can be carried out in a general form: 



dw 1 amcFz 



du v^v^ (l + u) 2 

dy 
yyfiny 



l + u 



d 

x J K2/3 (ay) JL- - 2a J K 1/3 (ay) <f\ny dy 



(4.30) 



where a = 2u/3xo- 111 the case when a> 1 the main contribution into the integral ()4.30|) 
gives the region y — 1 + ^, ^< 1. Taking the integrals over £ we obtain 

dw% F \piama z 1 + u + u 2 ( 2u\ . , 

T° 4 7 „(l + M )3 X ° eXP ('^j (4 ' 31) 

For round beams the integration over transverse coordinates is performed with the 
density 

27r^ eXP 2a 2 



Me) = 7T— r ex p(-^l ( 4 - 32 ) 



The parameter Xo(^) we present in the form 

Xo(o) = Xrd^-, x = —, f(x) = - (l -exp(-x 2 /2)) , 

Xrd = 0.720aiV c7 ^, /'(x ) = 0, /„ = f(x ) = 0.451256, (4.33) 

where iV c is the number of electron in the bunch. 

The Laplace integration of Eq. (|4.31|) gives for radiation intensity dl/du = eu/(l + 
u)dW/du 

d ~t H am2 ^{w\^ir$ f ° 2x%2 exp ("£) ' (434) 
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where f{j = f"(x ) = -0.271678. 

Integration of Eq. (|4.30[) over transverse coordinates gives the final result for the ra- 
diation intensity. For the round beams it is shown in Fig.23 for Xrd — 0.13, the curve 
attains the maximum at u ~ 0.02. The right slope of the curve agrees with the asymptotic 
intensity (J4.34)) and the left slope of the curve agrees with the standard classical intensity. 

dlcl = ^! 3 l/6 r(2/3)x 2/3 Ml/ 3 (435) 
du 71 

It will be instructive to compare the spectrum in Fig.23, found by means of integration 
over the transverse coordinates with intensity spectrum which follows from Eq. ()4.30|) 
(multiplied by uj) with averaged over the density Eq. ()4.32|) value xo Eq. (j4.33j) : Xo — 
Xrd ' 0.8135. The last spectrum reproduces the spectrum given in Fig.23, in the interval 
10~ 3 < u/xo < 1 with an accuracy better than 2% (near maximum better than 1%) while 
for u/xo < 10~ 3 one can use the classic intensity ()4.35|) and for u/xo 3> 1 the asymptotics 
()4.34|) is applicable. 

For the flat beams (a x 3> u y ) the parameter Xo{q) takes the form 

Xo = -77-^ = X m e~ v2 [e y erf (w) - ie x erf (iv)} , Xm = 2NcaiXc ; (4.36) 



here v = x/y/2a x , w = y/y/2a y , erf (z) = 2/y/ix J exp(— t 2 )dt, e x and e y are the unit 



vectors along the corresponding axes. The formula ()4.36|) is consistent with given in [HU] • 
In [7^] the term with was missed. Because of this the numerical coefficients in results 
for the flat beams are erroneous. 

To calculate the asymptotics of radiation intensity for the case u 3> Xm one has to 
substitute 



Xo = |Xol = X m e 



1/2 

erf 2 (w) — erf 2 (iv) 



(4.37) 



into Eq. (|4.31j) and take integrals over transverse coordinates x, y with the weight 



n±(x ' y) = 2^ exp ("^! - ^)' (438) 

Integral over x can be taken using the Laplace method, while for integration over y it is 
convenient to introduce the variable 

2 °° 

77 = — / exp(-t 2 )dt, w = -f-. (4.39) 
As a result we obtain for the radiation intensity in the case of flat beams 

du 8,/2(l-2/ 7 r) zXm « 3/2 (l + nY 6XP I, 3 Xm ) ' ^ 
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It is interesting to compare the high-energy end of intensity spectrum at collision of flat 
beams Eq. ([4.40)1 with intensity spectrum of incoherent radiation with regard for smallness 
of the transverse beam sizes considered above. For calculation we use the project TESLA 
parameters jST|: e = 250 GeV, a x = 553 nm, a y = 5 nm, a z = 0.3 mm, N = 2 ■ 
10 10 , Xm=0.13. The result is shown in Fig. 24, where the curves calculated according to 
Eq. (J4.40)) . and Eq. (J4.26j) . It is seen that for x = u/e > 0.7 the incoherent radiation 
dominates. For x > 0.7 the incoherent radiation may be used for a tuning of beams [77]. 

Along with the spectral characteristics of radiation the total number of photons emit- 
ted by an electron is of evident interest as well as the relative energy losses. We discuss 
an actual case of flat beams and the parameter Xm 1- In this case one can use the clas- 
sic expression for intensity (bearing in mind that starting from \ m ~ 1/10 the quantum 
effects become substantial). In classical limit the relative energy losses are 



' Ae\ 2am 2 
e J d 

Using Eqs.flOHJ), flj~37J) , and (jPHD we get 



(~) = ~3T~ / ' X 2 (t,x,y)n ± (x,y)dtdxdy. (4.41) 



Ae\ 2 [W 2 ao 



- 



, ow ■> — r=-jT (4-42) 
£ J d 9 V 3 7A C 



For mean number of emitted by an electron photons we find 
5am 2 



rf l = 

7 2V3e 



J x(t, x, y)n±(x, y)dtdxdy = 0.59275^ y^-rx™ = 1.072rx m . (4.43) 



If Xm > 1/10 one have to use the quantum formulas. For the energy losses one can use 
the approximate expression (the accuracy is better than 2% for any \) 



&/dt = 2/3amV [l + 4.8(1 + x)ln(l + 1.7x) + 2.44x 2 ] ^ . (4.44) 

Here x is t ne local value, so this expression for de/dt should be integrated over time and 
averaged over the transverse coordinates. For mean number of photons emitted by an 
electron there is the approximate expression (the accuracy is better than 1% for any xo) 



1.81xq7- 

[1 + 1.5(1 + X o)hi(l + 3xo) + 0.3x 2 ,] 1 / 6 ' 



Me) = h . ; ^ . r-TT77T^7-^ , r^rr^ (4-45) 



where the expression should be averaged over the transverse coordinates. For the project 
TESLA one gets for (Ae/e) c i = 4.3% according to Eq. 1)4.42)1 . while the correct result from 
Eq. (j4.44j) is Ae/e = 3.2%. For mean number of photons emitted by an electron we have 
correspondingly ~ 1.6 while the correct result is n 7 ~ 1.5. 

4.4 Mechanisms of e + e~- pair creation 

The probability of pair creation by a photon in the external field can be find from formulas 
(14.28)1 using the substitution rule (see e.g.JTT], [T0| ) : e — > —e, uj — > —u, uj 2 duj — > 
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—e 2 de. Performing these substitutions we obtain (see also Eq.(3.50) in [TP] ) 



dt "■ 2\/3 



7TCo> 



? + 7 j + / s = 3i77 K = Sflb (446) 

here e is the energy of the created electron, e' = u — e is the energy of created positron. 
There are different mechanisms of electron-positron pair creation 

1. Real photon radiation in the field and pair creation by this photon in the same field 
of the opposite beam. This mechanism dominates in the case x > 1- 

2. Direct electroproduction of electron-positron pair in the field through virtual photon. 
This mechanism is also essential in the case x ^ !• 

3. Mixed mechanism(l):photon is radiated in the bremsstrahlung process, i.e. inco- 
herently, and the pair is produced by photon in an external field. 

4. Mixed mechanism(2): photon is radiated via the beamstrahlung mechanism, and 
pair is produced in the interaction of this photon with individual particles of on- 
coming beam, i.e. in interaction with potential fluctuations. 

5. Incoherent electroproduction of pair. 

In the actual case x ^ 1 mixed and incoherent mechanisms mostly contribute. We 
start with the mixed mechanism (2). In the case Xm 1 the parameter k Eq. (J4.46|) 
containing the energy of emitted photon is also small and the incoherent cross section 
of pair creation by a photon is weakly dependent on the photon energy. To logarithmic 
accuracy 

a p = 28/9a 3 A 2 ln(<7„/A c ) (l + 396/t 2 /1225) . (4.47) 

If the term oc n 2 being neglected, the pair creation probability is factorized (we discuss 
coaxial beams of identical configuration). The number of pairs created by this mechanism 
(per one initial electron) is 

oo oo 

n p = J W 1 (g,t)dt J 2a p n ± (g,t')dt'd 2 g 

— oo t 

OO 

x J exp(— 3v 2 — 2u? 2 )[erf 2 (w) — erf 2 (iv) ^^dwdv 
o 

A 2 o 

= 0.1167a 3 Nrxm^- In ^. (4.48) 
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So for the project TESLA parameters the total number of produced pairs by this mech- 
anism is Nn™ ~ 1.4 ■ 10 4 per bunch in one collision. 

Now we turn over to discussion of incoherent electroproduction of pairs when both 
intermediate photons are virtual. To within the logarithmic accuracy for any \ and k one 
can use the method of equivalent photons 



duo 
where 



<J(T " n(uo)a p (uo), (4.49) 



n (u) = — In — , A = mil + k) 1/3 , 
ituo q m 

ijj f EY\ 1 / 3 1 

q m = m- 1 + + — =q F + q a . (4.50) 



e \ uj J o y 



Taking into account that the cross section of pair photoproduction is 

. . 28a 3 , m , m? , . 9/ o , , . . 

<T P (w) = -z — % ln — r~A ' u = — ( K < x ); a p K K ( K » !)> ( 4 - 51 ) 

»?7i qm{UJ ) UJ 

we obtain in the main logarithmic approximation for the cross section of the pair electro- 
production 

f n a \ ^ a f i m . m duo ,. 

a(2e^4e) = — — / hi— -In— — — , (4.52 
9m^ J q m {uJ) q m {v ) v 



u i ii 



where uj max = e/(l + x), ^ min = m 2 /uo max . If we put x = 0, a y = oo we obtain the 
standard Landau-Lifshitz cross section an (see e.g. [TT], [S]) 

28a 4 , o 9 . 

With regard for the bounded transverse dimensions of beam and influence of an external 
field the equivalent photon spectrum changes substantially. For x ~ 1 we have q{uo) = 
m(uj/e) 2 / 3 + l/a y and under condition r ) 2 ^ 3 \ c /a' y > 1 we find 

^ = ^1^7 4/3 -V n2 ?- ( 4 - 54 ) 
3m 2 \ a y J X c 

For TESLA parameters 7 2 ^ 3 A c /cr 2/ ~ 1 then 

^ = ^rV n V- 4.55 
81m 2 

This cross section is three times smaller than standard gll- F° r the project TESLA 
parameters the number of pairs produced by this mechanism is n v = La v /2 ~ 1.5 • 10 4 
in one direction per bunch in one collision, L = l/(47ra x a y ) is the geometrical luminosity 
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per bunch. So the both discussed mechanisms give nearly the same contribution for this 
project. It should be noted that the above analysis was performed under assumption 
that the configuration of beams doesn't changed during collision, although in the TESLA 
project the disruption parameter D y > 1. 

So we see that the beam-size effect and an external field effects are very essential for 
pair creation mechanism also. 
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A Appendix 



Basic equations 



The quasiclassical operator method developed by authors is adequate for consideration 
of the electromagnetic processes at high energy. The probability of the radiation process 
has a form (see jTUj, p. 63, Eq.(2.27)); the method is given also in jTT], for classical limit 

sec 



dw 



d 3 k 



(27T) 



10 



J dt 2 J dt 1 R*(t 2 )R(t 1 ) exp -j f (kx(t 2 ) - kx{t x )) 



(A-l) 



where k = (co>,k) is the 4-momentum of the radiated photon, k 2 = 0, x(t) = (t,r(t)), t 
is the time, and r(t) is the particle location on a classical trajectory, kx(t) = ut — kr(t), 
e' = e — u, let us remind that we employ units such that h = c = 1. The matrix element 
R(t) is defined by the structure of a current. For an electron (spin 1/2 particle) one has 



R(t) 
A(t) 



B(t) 



m 

; U 

e*p(t) 
2\fee' 

1 



Sf {p')e*u Si (p) = <f+ (A(t) + io-B(t)) ip s 



2Vee 



IE" 



m 



e + m 



+ 



IE" 



m 



e + m 



2£' V \'l , 



I E + m 
e' + m 



e* x p(t)) 



e' + m 



e* x (p(t) - hk)) 



(A.2) 



here p, p' = p — k is the initial and final momentum of an electron, e is the vector of the 
polarization of a photon (the Coulomb gauge is used), the four-component spinors u Sf , u Si 
and the two-component spinors (p Sf , tp Si describe the initial (sj) and final (sf) polarization 
of the electron, v=v(t) is the electron velocity, $ = v ~ l (v — n(nv)) ~ vj_, vj_ is the 
component of particle velocity perpendicular to the vector n = k/|k|, 7 = e/m is the 
Lorentz factor. The final expressions in (jA.2|) are given for radiation of ultrarelativistic 
electrons, they are written down with relativistic accuracy (terms ~ I/7 are neglected) 
and in the small angle approximation. If we are not interested in the initial and final 
particles polarization, then 



2s 



r> 



rl 



(A.3) 



where we have used the notation = i9(t 2 ), 1? = i?(£i). For polarized initial electron we 
have 



R*(t 2 )R(t 1 



2e'7 2 









IS* 













PP' + *j((p'-p) x V )C 



(A.4) 



82 



where £ is the vector describing the initial polarization of the electron (in its rest frame), 
p = 70(*i) and p' = T tf(t 2 ). 

If the particle moves along a definite trajectory, then substituting the classical values 
r(t) and v(t) in (jA.lJ) - (jA.3|) . we obtain the desired probability for this process. When 
there is scattering, Eq. fjA.lJ) must be averaged over all possible particle trajectories. This 
operation is performed with the aid of the distribution function, averaged over atomic 
positions in the scattering medium and satisfying the kinetic equation. We consider 
here general case with the external field (particle acceleration) taken into account. The 
emission probability per unit time is then (in this Appendix we follow the paper |9.) 

dW=(^\ = -y^—Re Jdrexp (-i-ur) [ d 3 vd 3 v'd 3 rd 3 r'£(ti' ',■&) 
\ dt I (27r) 2 uj J \ s J J 



x 



Fi(r, v, t)F f (r', v', r; r, v) exp (~ijM r ' ~ r )) ( A ' 5 ) 



The distribution function F in (jA.5|) satisfies the kinetic equation 

dF(r,v,t) dF(r,v,t) dF(r,v,t) 
dt dr dv 

= nj tr(v, v') [F(t, v', t) - F(r, v, t)] rfV, (A.6) 

where n is the number density of atoms in the medium, and er(v, v') is the scattering 
cross section. The normalization condition 

J d 3 r J d 3 vF(r,v,t) = 1 (A.7) 

should be also satisfied, as well as the initial condition for Ff 

F f (r', v', 0; r, v) = 5(r - r')5(v - v'). 

In Eq. (|A.5j) . we integrate over d 3 rd 3 (r—r'), taking advantage of the fact that Ff(r', v', r; r, v) 
can only depend on the coordinate difference r' — r: 



dW 



d 3 k 



(2tt) 2 uj 
where 

Fi(v,t) = J rf 3 rF,(r,v,t) 



ReJdrJ d 3 vd 3 v'£(& r ,&)Fi(v,t)F k (V,T;v), (A.8) 
o 



F k (V, r; v) = exp (-i—uT^j J d 3 r' exp f -z^k(r' - r)J F f (r', V, r; r, v) (A.9) 
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We can put 



F k (v',T;v) = U(#',#,T)5(\V\-\v\) 
C7(^', i?, 0) =5(tf'-tf), 



(A.10) 



if we take into account that in ultrarelativistic limit wv ~ 0(l/7 3 ), w ~ wj_, and the 
scattering cross section is 

a (v, v') = 5(\v'\ - |v|)ct« (A.ll) 
Making use of these results and Eq. ()A.6|) for Ff, we obtain the following equation for 

dU .us ( 1 Q/2 \ TT 8U 



n J d 2 d"o(ti\ •&") [U(0", , r) - [/« , r)] 



(A.12) 



where for scattering in the screened Coulomb potential the cross section in the Born 
approximation is 

AZ 2 a 2 

*) = 1^7 ^ , (A.13) 



1 



here $1 = — , a s is the screening radius (a s = 0.81asZ -1 / 3 , a# is the Bohr radius). 

A similar equation for the ^-dependence of U can be obtained by letting — > i? and 
w — > — w. 

If the final state of the charged particle is of no interest, the probability (jA.8|) must 
be integrated over d 3 v'. The resulting emission probability per unit time, normalized to 
a single particle moving at a speed v is then 



dW 



(2vr) 2 u 



■^Re J dr «p(-<^) [^>(tf,r) + (l + ^) *V(tf,r) 



(A.14) 



(jJTTl 

where a = — — , and = (Vo, V) satisfies the equation 

dV b 2 dV^ 

a7 + z /^- w W 



n J dV<r(0, V) [V^V, r) - V^(0, r)] , 



(A.15) 



here 6 = — . The initial conditions for V» are 



V O (0,O) = 1, V(0,O) = 0. 



(A.16) 
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If the scattering cross section cr(i?, "&) depends solely on the angle difference i? — 
as Eq. ()A.I3|) . then Eq. ljA.l 5|) is most conveniently solved by Fourier transforming with 
respect to the variable 

w( x > r ) = j^y J ^ ex P(^ x )^(#> r )> 

V^,r) = J c/ 2 xexp(-^x)^(x,r). (A.17) 

From Eq. (|A.15|) . with the initial conditions ()A.16|) . we have 
5^(x,r) ,b 

— i-A x Lp^(x, t) + zwx^ M (x, r) 

= (27r) 2 n [E(x) - £(0)] y? M (x, r), £(x) = J d 2 tfexp(ztfx)o-(tf); 

^o(x, 0) = (J(x), y>(x, 0) = -iV<?(x) (A.18) 

If the angular distribution of the radiation is of no interest, Eq. (|A.14|) must be integrated 
over the photon emission angles Bearing in mind that d 3 k = u 2 dujd 2q d, and that 

J 0V(tf,r)d 2 tf = -zVv>(x,t)| x=0 , (A.19) 

we obtain the following expression for the spectral distribution of emission probability per 
unit time 



00 

dW _ f , / ar\ 

atuRe / arexp I — 1- 



- 



du J V 2 / 





,2 / „2 



^ W (0,T)-i l + ^lV^r) 



£ ' 2 7 2 V ' e' 2 



(A.20) 



The equation like (|A.18jl in the classical limit and with w = was discussed in jlj, 
Appendix V. 

For the cross section (jA.13|) we have for £(x) ()A.18|) 

E(x) = ^^ Kl (x^), (A.21) 

where is the modified Bessel function. Bearing in mind that contributions to the 

cross section come from x ~ <C fii 1 , and expanding Ki(x'di) as a power series in x$\ 
we obtain the following equation for <^( x > T ) from Eq. ()A.18|) . to power-series accuracy: 

d(pn(x,T) .b 

— i-A x (p^(x, t) + 2\vxy? M (x, r) = it(x) 

u (x) = — x 2 f In -± + C - - J , (A.22) 
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where C = 0.577216... is the Euler constant. 

The "potential" u{x) is calculated in the Born approximation since Eq. (jA.13j) was used. 
The Coulomb correction to the potential u(x) were calculated in Appendix A of ^7]. The 
eikonal approximation was use which is valid for arbitrary value Za. Substituting this 
correction we have 



where the function f(Za) is defined in Eq. (l2.5j) . In absence of radiation and an external 
field (b=w=0) Eqs. fA~22|) and ifPTjl at the initial condition (p(x,0) = (2tt)" 2 give the 
solution of pure scattering problem (the initial condition for the distribution function 
F(0,t) is F(0,O) = (5(7?)): 



u(x) 



iZ 2 a 2 2 / xA „ 1 \ 
j 2 J ( In — + C - - + f(Za)\ , 



(A.23) 



<p(x,t) 



1 



u(x)t 



(27T) 2 




1 



(A.24) 



For another form of solution see Eq. (j3.25j) . 
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Figure captions 



• Fig.l The functions G(s) (curve 1) and 4>(s) (curve 2) in Eq. (j2.19j) . 

• Fig.2 The functions D^uq) (curve 1) and D 2 (u ) (curve 2) in Ea. (|2~2l?J) . 

• Fig. 3 The spectral intensity of radiation cu— — = x— — , x = — m gold m terms of 
3L mrf taken with the Coulomb corrections (see Eq. (j2.36j) ). 

— Curve BH is the Bethe-Maximon spectral intensity (see Eq. (j2.29)) ): 

— curve 1 is the logarithmic approximation udW c /du Eq. (|2.17jl . curve cl is the 
first correction to the spectral intensity udWi/du Eq. (j2.22j ) and curve Tl is 
the sum of the previous contributions for the electron energy e = 25 GeV; 

— curve 2 is the logarithmic approximation udW c /duj Eq. (|2.17jl . curve c2 is the 
first correction to the spectral intensity LodWi/du Eq. (j2.22|) and curve T2 is 
the sum of the previous contributions for the electron energy e = 250 GeV; 

— curves 3, c3, T3 are the same for the electron energy e = 2.5 TeV; 

— curves 4, c4, T4 are the same for the electron energy e = 25 TeV. 

• Fig. 4 The relative energy losses of electron per unit time in terms of the Bethe- 
Maximon radiation length L® ad : —L° rad in gold vs the initial energy of electron (curve 

1) and the total pair creation probability per unit time W£ (see Eq. (l2.49j) fin terms 
of the Bethe-Heitler total probability of pair creation W^ H (see Eq. ()2.47jl ) in gold 
vs the initial energy of photon (curve 2). 

• Fig. 5 The pair creation spectral probability p , x = — in gold in terms of the 

dx uj 

exact total Bethe-Maximon probability taken with the Coulomb corrections (see 
Eq.(E3ZD)- 

— Curve BH is the Bethe-Maximon spectral probability (see Eq. (J2.46|0 ; 

— curve Tl is the total contribution (the sum of the logarithmic approxima- 
tion dWp/de (J2.43|) and the first correction to the spectral probability dW^/de 
(j2.45j) ) for the photon energy uj = 2.5 TeV; 

— curve 2 is the logarithmic approximation dW^/de (|2.43|) . curve c2 is the first 
correction to the spectral probability dW^/de (|2.45j) ) and curve T2 is the sum 
of the previous contributions for the photon energy u = 25 TeV; 

— curves 3, c3, T3 are the same for the photon energy uj = 250 TeV; 

— curves 4, c4, T4 are the same for the photon energy uj = 2500 TeV; 
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Fig. 6 Anomalous magnetic moment (AMM) of electron in units a/2n (r in Eq. (j2.69j0 in 
gold vs electron energy in TeV. 

Fig. 7 The functions Re k+ + (curve 2) and Im k+ + (curve 1) versus the photon energy 
taken in units u e (because of this the curves are universal) for the case when the 
influence of a medium is taken into account only (Eq. (|2.91Jl ). The both curves are 
normalized to the asymptotics given by Eq. (|2.94J) in the limit ujp — > oo. 



Fig. 8(a) The amplitude M ++ of the coherent photon scattering in gold under in- 
fluence of the multiple scattering at the different momentum transfer to the photon 
A in terms of the amplitude ImM ++ (|2.1()fi|) calculated for the screened Coulomb 
potential. 

- Curve 1 is ImM ++ for A =0.4435 m. 

- Curve 2 is ReM ++ for A =0.4435 m. 

- Curve 3 is ImM ++ for A =0.0387 m. 

- Curve 4 is ReM ++ for A =0.0387 m. 

Fig. 8(b) The same as in Fig.l but for the amplitude M + _ in terms of the amplitude 
ImM + _ ()2.106|) calculated for the screened Coulomb potential. 

Fig. 9 The reduction factor Jbh for the energy of the initial electron e = 25 GeV. 
The curves 1,2 and 3 are respectively for j3 =0.01, 0.03 and 0.06. 

Fig. 10 The reduction factor Jlpm for energy of the initial electron £=25 GeV in 
tungsten (cu c =228 MeV). The curves 1, 2 and 3 are respectively for /?=0.01, 0.03 
and 0.06. 

Fig. 11 The energy losses spectrum — in units — , in the target consisting of two 

du 71 

gold plates with thickness l\ = 11.5 fim for the initial electrons energy e=25 GeV . 

— Curve 1 is for distance between plates l 2 = 2l\\ 

— Curve 2 is for distance between plates l 2 = 4Zi; 

— Curve 3 is for distance between plates l 2 = 6I1, 

— Curve 4 is for distance between plates l 2 = 8I1, 

- Curve 5 is for distance between plates l 2 = 101 1. 

de 

Fig. 12 The energy losses — in tungsten with thickness I = 0.088 mm in units 

du 

2a 

— , ((a) is for the initial electrons energy e=25 GeV and (b) is for e=8 GeV). The 

7T 

Coulomb corrections and the polarization of a medium are included. 

— Curve 1 is the contribution of the main term describing LPM effect; 
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— curve 2 is the correction term; 

— curve 3 is the sum of the previous contributions; 

— curve 4 is the contribution of the boundary photons; 

— curve 5 is the sum of the previous contributions; 

— curve T is the final theory prediction with regard for the reduction factor (the 
multiphoton effects). 

Experimental data from Fig. 9 of jH]. 

Fig. 13 The energy losses u— — in gold with thickness I = 0.023 mm in units — , 

dtu ' it 
((a) is for the initial electrons energy e=25 GeV and (b) is for e=8 GeV). The 

Coulomb corrections and the polarization of a medium are included. 

(2*1 (2) 

— Curve 1 is the contribution of the term Re = Re J3 ; 

(2) 

— curve 2 is the contribution of the term Re J% ; 

— curve 4 is the contribution of the term Re jf^ , all ()3.29|) ; 

— curve S is the sum of the previous contributions Re 

— curve 5 is the contribution of the boundary photons, see Eq.(2.20), [THj : 

— curve T is the total prediction for the radiation energy losses . 

Experimental data from Fig. 12 of 



Fig. 14 The energy losses u— — in gold with thickness I = 0.0038 mm in units 

dio 

2a 

— , ((a) is for the initial electrons energy £=25 GeV and (b) is for e=8 GeV). The 
Coulomb corrections and the polarization of a medium are included. 

(2) (2) 

— Curve 1 is the contribution of the term Re J{ = Re J3 ; 

(2) 

— curve 2 is the contribution of the term Re J2 ; 

— curve 4 is the contribution of the term Re J\ , all ()3.29|) : 

— curve S is the sum of the previous contributions Re J^; 

— curve 5 is the contribution of the boundary photons see Eq.(2.20), 

— curve T is the total prediction for the radiation energy losses . 

Experimental data from Fig. 13 of 



de 

Fig. 15 The energy loss — in iridium target with thickness / = 0.128 mm in units of 

inverse radiation length L® ad fEq. ()2.36|0 for the initial electrons energy e=287 GeV. 
The Coulomb corrections are included. 
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— curve 1 is the the Bethe-Maximon intensity spectrum; 

— curve 2 is the contribution of the main (Migdal) term describing LPM effect; 

— curve 3 is the correction term; 

— curve 4 is the sum of the previous contributions; 

— curve T is the final theory prediction with regard for the reduction factor (the 
multiphoton effects). 

Experimental data from Fig. 2 of [To] . 

Fig. 16 The same as in Fig. 15 but for the initial electrons energy e=207 GeV. 

Fig. 17 The same as in Fig. 15 but for the initial electrons energy e=149 GeV. 

Fig. 18 The reduction factor for iridium target with thickness I = 0.128 mm (4.36 
% L ra( i) versus x = u/e. The curves 1, 2, 3 are for energies 149 GeV, 207 GeV and 
287 GeV correspondingly. 

Fig. 19 The relative energy losses of electron per unit time including the contribution 
of the correction term Eq. ()2.17|) in terms of the Bethe-Maximon radiation length 

L ra d Eq. (j2.36|) : -L ra d in iridium (curve 1) and in lead (curve 2) vs the initial energy 
of electron. 

Fig. 20 The bremsstrahlung intensity spectrum ujda/duj in units 2ar\ versus the 
photon energy in units of initial electron energy (x = oj/s) for VEPP-4 experiment. 
The upper curve is the standard QED spectrum, the three close curves below are 
calculated for the different vertical dimensions of colliding beams (equal for two 
colliding beams a = a z = A z ):a = 20 fim (bottom), a = 24 \im (middle), o = 27 \xm 
(top). The data measured in [73] are presented as circles (the experiment in 1980) 
and as triangles (the experiment in 1981) with 6 % systematic error as obtained in 

Fig. 21 The bremsstrahlung intensity spectrum ujda/duj in units 2ar\ versus the 
photon energy in units of initial electron energy [x = uj/e) for the HERA experiment. 
The upper curve is the standard QED spectrum, the two close curves below are 
calculated with the beam-size effect taken into account: the bottom curve is for 
a z = A z = 50 fim, o y = A y = 250 \im\ while the top curve is for o z = A z = 
54 /xm, (x y = A y = 250 fim. The data taken from Fig.5c in [TH] 

Fig. 22 The bremsstrahlung intensity spectrum ujda/duj in units 2otq versus the 
photon energy in units of initial electron energy (x = oo/e) for linear collider with 
beam energy e = 250 GeV. The upper curve is the standard QED spectrum. The 
curve below is calculated with the beam-size effect taken into account for a x = 
553nm, a y = 5nm. 
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• Fig. 23 Spectral intensity of radiation of round beams in units am 2 a z for x r d=0.13 
calculated according to Eqs. (j4.30J) . (j4.33|) . 

• Fig. 24 The spectral radiation intensity dl/du of coherent radiation (fast falling 
with x = oj/e increase curve) and of incoherent radiation (the curve which is almost 
constant) in units Na 2 X c /a x for beams with dimensions a x = 553 nm, o y = 5 nm, 
for Xm=0.13. 
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TABLE 1 The characteristics of the LPM effect 
Listed are the charge of nucleus Z, the density g, the radiation length L rad , 
e e (see Eq. flZZfr ), ^e(see Eq. jjQHl) ), cj c Eq. (jOZj) ) , cu p Ea. (|ZlT5jl ). 
Z is the target thickness in percent of the radiation length L rad 
(all photon energies u; c and u; p are for the electron energy e = 25 GeV) 



Material 


Z 


g(g/cm 3 ) 


L rad (cm) 


^e(TeV) 


u; e (TeV) 


u c (MeV) 


u p (MeV) 


/ 


C 


6 


2.2 


18.8 


144 


580 


4.3 


1.5 


2.1, 6 


Al 


13 


2.7 


8.9 


68 


270 


9.2 


1.6 


3.5, 6 


Fe 


26 


7.87 


1.76 


14 


56 


46 


2.7 


2.8, 6.1 


Pb 


82 


11.35 


0.56 


4.3 


17.2 


150 


3.0 


2.7 


W 


74 


19.3 


0.35 


2.7 


10.8 


230 


3.9 


2.7, 6.4 


u 


92 


18.95 


0.32 


2.5 


10 


250 


3.8 


2.2, 4.2 


Au 


79 


19.32 


0.33 


2.5 


10 


240 


3.9 


0.1, 0.7 


Water 




1 


36.1 


277 


1100 


2.6 






Std.rock 


11 


2.65 


10.0 


77 


300 


8.1 







TABLE 2 Characteristic parameters of the radiation process in 
tungsten with the thickness I = 2.7%L rad 



e (GeV) 


u c (MeV) 


lu p (MeV) 


T 

A c 


ujt (MeV) 




uo m (MeV) 


25 


230 


3.9 


21.3 


1.6 


2.7 


2 


8 


23 


1.3 


21.3 


0.76 


5.7 


0.5 



TABLE 3 Characteristic parameters of the radiation process 
in gold with the thickness I = 0.7%L rad and I = 0.1%L rad 
(all photon energies u are in MeV) 



e (GeV) 




u p 


T c (0.7) 


T c (0.1) 


^(0.7) 


/U0.7) 


An(0.1) 




25 


240 


3.9 


5.82 


0.96 


1.6 


0.75 


0.12 


28 


8 


25 


1.3 


5.82 


0.96 


0.76 


1.6 


0.25 


3.0 
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(b) Tungsten, 8 GeV 
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